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APRESENTAÇÃO 


Um dos grandes desafios enfrentados atualmente pelos engenheiros nos mais 
diversos ramos do conhecimento, é de saber ser multidisciplinar, aliando conceitos de 
diversas áreas. Hoje exige-se que os profissionais saibam transitar entres os conceitos e 
práticas, tendo um viés humano e técnico. 

Neste sentido este livro traz capítulos ligados a teoria e prática em um caráter 
multidisciplinar, apresentando de maneira clara e lógica conceitos pertinentes aos 
profissionais das mais diversas áreas do saber. 

Apresenta temas relacionados a área de engenharia mecânica e materiais, dando 
um viés onde se faz necessária a melhoria continua em processos, projetos e na gestão 
geral no setor fabril. Destaca-se ainda a apresentação das áreas da engenharia e elétrica e 
eletrônica, com a busca da redução de custos e automação de processos. 

Da ênfase em alguns trabalhos voltados a realizar um levantamento econômico dos 
de processos e o estudo das áreas térmicas. 

De abordagem objetiva, a obra se mostra de grande relevância para graduandos, 
alunos de pós-graduação, docentes e profissionais, apresentando temáticas e metodologias 
diversificadas, em situações reais. 

Aos autores, agradeço pela confiança e espirito de parceria. 

Boa leitura 


Henrique Ajuz Holzmann 
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RESUMO: Fazendo uso do microcontrolador 
Arduino um shopping center será automatizado, 
desde setores mais simples até os mais 
complexos. Utilizando nove microcontroladores, 
os quais trocarão informações com um 
supervisório, nele será feito todo o monitoramento 
e controle do sistema, tendo em vista que todo o 
sistema funcionará de forma automática, o qual 
terá intervenção humana somente para ajustes 
e acionamento manual se necessário. Fazendo 
uso de sensores de luminosidade a iluminação 
do ambiente será monitorada e controlada. 
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SHOPPING CENTER 


Sendo esse um sistema de baixo custo, fácil 
acesso comercial, que atende as necessidades 
dos ambientes e proporciona bem-estar aos 
usuários. 
PALAVRAS-CHAVE: 
Arduino, Periférico. 


Automação, Sistema, 


STUDY AND IMPLEMENTATION OF 
MICROCONTROLLERS IN SHOPPING 
CENTER AUTOMATION 
ABSTRACT: Using the Arduino microcontroller, 
a shopping center will be automated, from the 
simplest to the most complex sectors. Using nine 
microcontrollers, which will exchange information 
with a supervisor, all monitoring and control of 
the system will be carried out there, considering 
that the entire system will work automatically, will 
have human intervention only for adjustments 
and manual activation if necessary. Using light 
sensors, the ambient lighting will be monitored 
and controlled. This being one low cost system, 
easy commercial access, which meets the needs 
of the environments and provides well-being to 

users. 
KEYWORDS: Automation, System, Arduino, 
Peripheral. 


11 INTRODUÇÃO 


Muito se ouve falar em automação e 
sistemas automatizados, ela está presente no 
dia a dia do ser humano, tanto na indústria, no 
comércio e nas residências. É uma tecnologia 


relacionada com a aplicação de sistemas 
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mecânicos, elétricos e eletrônicos, apoiados em meios computacionais, na operação e 
controle de sistemas diversos. Essa tendência tecnológica tem crescido e tem proporcionado 
melhorias nos processos e na economia (ROSÁRIO, 2009). 

As plataformas de microcontroladores de baixo custo, tipo Arduino e seus 
equivalentes, facilitam o acesso tecnológico da eletrônica e automação. Elas possibilitam a 
criação de sistemas automatizados de baixo custo, quando comparados aos equipamentos 
tradicionais utilizados atualmente. Arduino é um microcontrolador que engloba software 
e hardware livres (open source) que proporciona aos usuários uma plataforma para 
prototipação de projetos interativos, englobando a área da computação física, onde o 
software interage com sensores, atuadores e outros dispositivos eletrônicos (MCROBERTS, 
2015). 

Visando o conforto e bem-estar dos clientes surge a ideia de automatizar todo o 
sistema de iluminação de um shopping center. Dessa forma, o objetivo desse trabalho 
é desenvolver um sistema de automação fazendo uso do Arduino e seus periféricos de 
entrada e saída., que venha a atender todas as necessidades sem fazer uso de um CLP e 
seus periféricos, de forma a gerar economia para a empresa na implementação. 


2 | MATERIAIS E MÉTODOS 


Na implementação desse trabalho será usado o microcontrolador Arduino mega 
2560. O Arduino mega 2560 é baseado no microcontrolador ATmega2560, possui 16 
entradas analógicas, 54 pinos de entrada e saídas digitais. Dessas, 15 podem ser usadas 
como saídas PWM. Possui também 4 portas de comunicação serial (BLUM, 2016). 





Figura 1: Arduino Mega 2560. 


Fonte: Autoria própria. 


A programação do sistema é realizada no Arduino IDE. O Arduino IDE é um software 
Arduino totalmente gratuito que permite desenvolver e gravar códigos no microcontrolador. 
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Fazendo uso do mesmo, pode-se realizar o upload dos códigos para a placa. Caso ocorra 
algum erro no código ou na configuração, o software gera uma notificação informando o 
mesmo (MONK, 2017). 


e sketch mar20a | Arduino 1.0.3 Ed 
File Edit Sketch Tools Help 


sketch mar20a 


A 








Arduino Uno on COMA 











Figura 2: Arduino IDE. 


Fonte: Autoria própria. 


O código da programação é realizado na linguagem C++ (MONK, 2017). 

Fazendo uso dessas tecnologias, de suas programações e dos módulos sensor 
de luminosidade, os quais podem ser configurados para diferentes níveis de detecção. E 
também fazendo uso dos periféricos de entrada e saída, sabendo a função de cada um 
deles, será desenvolvido na plataforma a programação de acordo com a necessidade de 
cada ambiente e do próprio usuário. 

Na figura 3 é mostrado o sensor de luminosidade, no qual pode-se configurar a faixa 


de lúmens. 
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Figura 3: Luxímetro TSL2165. 


Fonte: Arduinomega. 


Na figura 4 e 5 temos os módulos relés que serão usados para acionar os contatores 
do sistema de iluminação. 





Figura 4: Modulo relé 2 canais. 


Fonte: Arduinomega. 
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Figura 5: Modulo relé 6 canais. 


Fonte: Arduinomega. 


Na figura 6 mostra o diagrama de comunicação entre os Arduinos e o supervisório. 


SUPERVISÓRIO 


Arduino 07 


Arduino 03 
Arduino 04 


Arduino 06 


Figura 6: Diagrama do sistema. 


Fonte: Autoria própria. 


Para supervisionar todo esse sistema, será implementado um supervisório para ter 
todo o controle. 
A seguir na figura 7 temos uma tela touch screen a qual ira mostra todo o sistema e 


seus respectivos dados 
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Figura 7: Tela para o supervisório. 


Fonte: Arduinomega. 


Além disso, esse sistema pode ser acessado remotamente e operado da mesma 
forma. Essa automatização pode ser usada para outras aplicações que vão bem além da 


iluminação. 


31 RESULTADOS E DISCUSSÃO 


Espera-se um sistema que venha a atender todas as necessidades voltadas a 
automação que envolve um shopping center, desde aplicações mais simples até as mais 
complexas. Um sistema que minimiza os custos, tendo em vista que, fazendo uso desses 
microcontroladores a implementação terá um custo baixo comparado ao uso de CLP. 

Sabedores também que ocorrerá uma otimização do tempo, pois não será mais 
necessário a intervenção humana para acionamento de setores, sendo que o sistema 
funcionará de forma automática. Utilizando o supervisório será feito todo o acompanhamento 
através de dados gerados pelo mesmo. 

Isso gerará uma satisfação para o shopping e para os seus clientes, pois se tornará 


um ambiente mais confortável e com uma iluminação adequada. 


41 CONCLUSÃO 


O resultado esperado desse projeto é conforto e comodidade ao cliente, junto a 
isso, uma melhoria no setor de operações, tendo em vista que, ocorrerá uma otimização 
do tempo nas atividades, pois, a maioria dos acionamentos é feito através de timer, sendo 
necessário muitas vezes o ajuste devido as mudanças de luminosidade do ambiente. 

É importante salientar que é um sistema desenvolvido com tecnologias de fácil 
acesso comercial, fácil utilização e baixo custo. Atendendo assim, a necessidade dos 


shopping centers. 
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RESUMO: Este artigo investiga o desempenho 
do inversor de três níveis NPC e do inversor 
HB. A implementação é baseada na modulação 
Level-Shift PWM. Uma análise dos inversores é 
feita com base no desempenho (i) de Qualidade 
de Energia, avaliando o THD de corrente; (li) de 
Eficiência, avaliando as perdas de condução e 
chaveamento e (iii) de Confiabilidade, avaliando 
o estresse térmico. A quantidade de dispositivos 
assim como o padrão de chaveamento 
influenciará o desempenho do inversor. 
PALAVRAS-CHAVE: Inversor Multinível. 
Eficiência. Qualidade de Energia. 


PERFORMANCE ANALYSIS OF NPC AND 
H BRIDGE THREE LEVEL INVERTERS 


ABSTRACT: This article investigates the 
performance of the NPC three-level inverter and 
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the HB inverter. The implementation is based on 
Level-Shift PWM modulation. An analysis of the 
inverters is made based on the Power Quality 
performance (i), evaluated in current THD; (ii) 
Efficiency, assessed as reduced conduction and 
switching, and (iii) Reliability, assessed or thermal 
stress. A number of devices as the switching 
pattern influences drive performance. 
KEYWORDS:  Multilevel Inverter. 
Power quality. 


Efficiency. 


11 INTRODUÇÃO 


A preocupação com o meio ambiente 
e desejo de libertação da dependência dos 
combustíveis fósseis está conduzindo cada 
vez mais o aprimoramento de técnicas e 
soluções alternativas para substituir as fontes 
convencionais e poluentes de energia. Carros 
elétricos, geração de energia através do vento 
(eólica) e através do sol (fotovoltaica) são as 
alternativas mais atraentes no contexto de 
sustentabilidade (ANDRADE et al., 2012). 
Todavia, para a utilização dessas alternativas 
há a necessidade de conversão de energia de 
CC para CA. Os inversores têm como objetivo 
justamente isso, de uma maneira geral os 
inversores são equipamentos que convertem 
energia contínua (CC) em energia alternada 
(CA) e vice-versa. Também há inversores que 
têm como função elevar ou diminuir a tensão, 
sem convertê-la para CC ou CA (FRANQUELO 
et al., 2008). 
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Os inversores são constituídos por semicondutores (na maioria dos casos MOSFETs) 
que fazem o chaveamento da energia que chega na carga, controlando o fluxo de tensão 
e corrente. Há várias formas de organizar os semicondutores nos inversores, cada forma 
de organização exige uma quantidade diferente de componentes e efeitos na potência de 
saída. Na literatura isso é chamado de topologia. Um desses efeitos é o nível de tensão que 
se obtém na saída do inversor. Com o chaveamento em execução, a tensão de saída varia 
em determinados níveis. No inversor de dois níveis, ora a tensão de saída será máxima, 
ora a tensão será mínima. No de três níveis a tensão será máxima (P), zero(O) ou mínima 


(N), como mostra a figura 01. 





Tempo (s) 


Figura 01 — Tensão na Carga do Inversor de Três Níveis 
Fonte: Autor (2019). 


O inversor de dois níveis é o mais comum e mais usado, devido ao fato de ser 
O primeiro inversor desenvolvido, utilizar poucos semicondutores e ser o mais básico 
construtivamente. Contudo, ele tem a desvantagem de utilizar semicondutores de alta 
potência, devido à tensão máxima reversa suportada pelo o mesmo quando se está em 
aberto, mesmo sendo usado para aplicações de baixa potência (DE OLIVEIRA., 2005). 

Os inversores multiníveis são os inversores que possuem mais de dois níveis 
na tensão de saída e vêm com objetivo de diminuir essas desvantagens que o de dois 
níveis possui. As vantagens que os multiníveis têm sobre os de dois níveis são focadas 
na melhoria da potência de saída e qualidade do sinal. Quanto maior o número de níveis 
do inversor, melhor a forma de onda da saída e menor a distorção harmônica total (THD). 
As pesquisas mostram também um equilíbrio maior da ondulação da corrente da saída 
(CELANOVIC; BOROYEVICH, 2001; SAEEDIFARD; IRAVANI; POU, 2007). 

Duas das topologias mais comuns dos inversores são Ponto Neutro Grampeado 
(NPC) e o Ponte H em Cascata (HB). O NPC foi apresentado pela primeira vez por A. Nabae, 
|. Takahashi e H. Akagi em 1981 (NABAE; TAKAHASHI; AKAGI, 1981). Essa topologia de 
inversor foi baseada no inversor de dois níveis, sendo acrescentado dois semicondutores 
em cada fase e sendo controlado por modulação de largura de pulso (PWM). A Figura 02 nos 


mostra um inversor NPC de três níveis e um inversor de dois níveis, ambos monofásicos. 
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O inversor HB (Figura 03) surgiu em 1992 por M. Marchesoni, M. Mazzucchelli 
e S. Tencon. Tais topologias têm características diferentes entre si, como número de 
componentes, modularidade PWM, eficiência, qualidade do sinal de saída, THD, perda 
por estresse térmico e tolerância a falhas. Dependendo da aplicação, a topologia dos 
inversores multiníveis pode ser definida levando em consideração esses fatores. 

Na literatura podemos encontrar trabalhos que analisam o desempenho dos 
inversores com base nas perdas de condução e chaveamento. No entanto, a análise com 
base no estresse térmico não está bem difundida. Neste artigo discutiremos os resultados 


e características obtidos por simulação do inversor de três níveis para as topologias NPC 
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Figura 02 — Inversor de três níveis, NPC e Figura 03 — Inversor multinível HB. 


inversor de dois níveis, respectivamente. 


21 REFERENCIAL TEÓRICO 


A comparação dos inversores se baseia nos níveis de tensão, consequentemente 
de corrente e potência, na qualidade de energia, eficiência e confiabilidade. Além dessas 
características o número de componentes é um aspecto que impacta em alguns dados de 


desempenho. Neste tópico será apresentado as topologias NPC e HB. 


2.1 O Inversor NPC 


O inversor NPC (ver figura 02) foi proposto por Nabae em 1981 (NABAE; TAKAHASHI; 
AKAGI, 1981). O inversor possui 4 chaves ativas e 2 diodos. O barramento CC é dividido 
com dois capacitores, cuja tensão de cada um é da metade da tensão total do barramento. 

Quando se utiliza uma única fonte, o equilíbrio da tensão nos capacitores do NPC é 
um problema relevante, especialmente nos inversores com maior número de níveis, porque 
o fluxo de potência ativa tende a desbalancear a tensão nos capacitores (ANDRADE et al., 
2012)(OLIVEIRA; CORREA, 2012). 

A tabela 01 apresenta os estados de chaveamento para o inversor NPC. Para a 
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modulação deste inversor utilizou-se a técnica de modulação Level-Shift PWM (LS-PWM) 
(FRANQUELO et al., 2008). 


Vo (V) 





Tabela 01 — Estado de Chaveamento do Inversor NPC 


A Figura 04 apresenta a estratégia PWM que se resume a comparação da tensão de 
referência com duas ondas triangulares. Considerando que Vdc = 100V, a comparação da 
senoide de referência (Vref = -50V) com a triangular superior (varia entre OV e 50V) define 
a modulação das chaves S1 e S3. Já a comparação com da tensão de referência com a 
triangular inferior (que varia entre -50 V e OV) define o chaveamento das chaves 32 e S4. 








Figura 04 — Modulação PWM para o Inversor NPC 
Fonte: Autor. (2019) 


2.2 O Inversor HB 


O inversor HB (ver figura 03), mais conhecido como ponte-H, possui apenas 4 
chaves ativas e uma única fonte. A tensão de três níveis é obtida por meio do chaveamento 
adequado de suas chaves (OLIVEIRA; CORREA, 2012). 

A tabela 02 apresenta os estados de chaveamento para o inversor HB. Para 
a modulação PWM do inversor foi utilizada duas tensões de referência e uma tensão 
triangular como portadora. 

Cada tensão de referência é utilizada para modular a tensão de um braço. Buscando 
um maior aproveitamento da tensão de referência, a tensão de referência do braço A é 
defasada de 180º da tensão de referência do braço B. A figura 05 apresenta o esquema da 
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modulação. 





Tabela 02 — Estado de Chaveamento do Inversor HB 
Fonte: Autor (2019). 
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Figura 05 — Modulação PWM para o Inversor HB 
Fonte: Autor. (2019) 


31 MÉTODO DA PESQUISA 


A simulação foi realizada com base na variação do índice de modulação em 
amplitude (0,1 < ma < 1,0) em dois casos de carga. A carga 01 com o fator de potência de 
0,5 e a Carga 2 com o fator de potência 0,75. Os valores da carga estão apresentados na 
tabela 08. 


L 


ART LS 
1,560 | 7,18mH 





Tabela 03 — Cargas e seus respectivos valores de resistência e indutância 
Fonte: Autor (2019). 


4 | RESULTADOS 


Os resultados obtidos foram avaliados com base em três aspectos, que são: 
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qualidade de energia; eficiência dos inversores e confiabilidade. 

A figura 06 apresenta a forma de onda da corrente na carga na topologia HB e NPC 
com o índice de modulação 1,0 e com o COS 0,75. 

A figura 07 apresenta a forma de onda da tensão na carga do HB e NPC com o 
índice de modulação 1,0 e com o COS 0,75. 


Corrente 


Corrente 





Figura 06 — Forma de onda da corrente do HB e NPC, ambos com COS 0,75. 
Fonte: Autor. (2019) 


4.1 Qualidade de energia 


Para avaliar o desempenho dos inversores com base na qualidade de energia, é 
utilizada a Taxa de Distorção Harmônica. A Taxa de Distorção Harmônica (do inglês, Total 
Harmonic Distortion - THD) é um critério muito utilizado para definir o conteúdo harmônico 
de um sinal alternado. A THD é definida como: 


100 [ao 
THD% = Lito V (1) 


A THD define o quanto que o sinal difere de uma onda senoidal pura. A presença de 


uma THD alta tem influência em vários pontos no estudo de inversores. 
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Figura 07 — Forma de onda da tensão do HB e NPC, ambos com COS 0,75 
Fonte: Autor. (2019) 


Primeiramente que a qualidade da energia diminui consideravelmente. Outro ponto 
em que a THD tem influência é no fator de potência, por isso a importância do cálculo da 
THD. 

A figura 08 apresenta os resultados da simulação para a THD de corrente para as 
duas topologias nas duas condições de carga. 

Observando a figura 09, percebe-se que a THD é maior com o fator de potência de 
0,75 nas duas topologias apresentadas. 

É interessante de perceber que durante o intervalo aproximado de 0,8 < índice de 
modulação < 1,0, o inversor HB possui maior THD. Quando o índice de modulação é igual a 
1 a diferença do THD é de 61,5%. No intervalo de 0,1 < índice de modulação < 0,8, o NPC 


apresenta maior THD com uma diferença de 58,7% com o índice de modulação igual a 0,1. 
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Figura 08 - THD de Corrente 
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Figura 09 - THD dos inversores com COS 0,75 


4.2 Eficiência dos inversores 


Para calcular a eficiência dos inversores, serão analisadas as perdas por condução 
e chaveamento. Para diferentes topologias de inversores ou diferentes tipos de técnicas 
de modulação proposta, torna-se muito importante determinar as perdas de potência nos 
interruptores. 

Entende-se por perdas totais o somatório das perdas por condução e perdas por 
chaveamento. Na figura 10 e 11 mostra o gráfico de perdas totais com COS 0,5 e COS 0,75 
para as topologias NPC e HB, respectivamente. 

Nas duas topologias percebemos a característica de que com um COS 0,5, há 
maiores perdas. Na figura 12 temos um comparativo com as duas topologias com o COS 
0,5. O NPC possui maiores perdas que o HB, contudo a diferença de perdas entre as duas 


topologias diminui com a diminuição do índice de modulação. 
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Figura 10 — Perdas totais do NPC 
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O aumento das perdas na topologia NPC se dá pelo fato de a topologia possuir 4 


chaves e dois diodos, enquanto a topologia HB possui apenas quatro chaves. 
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Figura 11 — Perdas totais do HB 


Outro ponto importante de análise são as perdas de condução e de chaveamento 


apresentadas de forma separadas. 
A figura 13 apresenta um comparativo das perdas de condução entre as topologias 


NPC e HB. 
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Figura 12 — Perdas totais dos inversores com COS 0,5 


Pode-se perceber que as perdas por condução de na topologia NPC são maiores. 
Isso se dá pelo fato de a corrente no inversor NPC ser maior, como pode ser vista pela 


figura 06. 
Em relação às perdas por chaveamento, pode-se perceber, pela figura 14, que as 
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maiores perdas estão na topologia HB. Isso se dá pelo fato de o chaveamento na chave 
desta topologia ser maior do que a topologia NPC, como pode ser visto pela figura 15. 
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Figura 15 — Tensão na chave S1 de cada topologia 


4.3 Confiabilidade dos inversores 


A confiabilidade é caracterizada pelo estresse térmico do inversor, quanto menor for 
o estresse que ele sofre maior a confiabilidade do inversor. Considerando que os inversores 
estão cada vez mais sendo utilizados em uma grande diversidade de aplicações, a garantia 
da confiabilidade desse dispositivo se torna algo imprescindível. 

Como indicador de confiabilidade, a tensão térmica é a variável mais adequada 
para avaliar Os inversores de potência. De acordo com as estatísticas realizadas por (NAMI 
et al., 2011), a proporção de vários estresses causa falhas nos componentes eletrônicos, 
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entre os motivos, destaca-se o estresse térmico que é responsável por 55% das falhas nos 
componentes. 

Nas figuras 16 e 17 vemos os resultados das simulações para o estresse térmico em 
cada topologia. A diferença de estresse térmico entre o COS 0,5 e COS 0,75 é pequena, 
com o estresse sendo maior para o COS 0,5, para todas as topologias. 

Na figura 18 temos o pior caso para todas as topologias, que é com o COS 0,5, 
podendo-se ver que o estresse térmico das topologias é quase igual, tendo maior diferença 
com o índice de modulação entre 0,4 e 0,5. 
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Figura 18 — estresse térmico dos inversores com COS 0,5 


5 | CONCLUSÃO 


Os inversores multiníveis são muito importantes para aplicação industrial. A 
depender de como a topologia é formada pode-se ter um inversor que tenha desempenho 
melhor na qualidade da energia, na eficiência e confiabilidade. 

As topologias NPC e HB demostraram resultados aproximados no que diz respeito 
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ao estresse térmico. No entanto, pelo fato da topologia NPC possuir mais dispositivos teve 
uma eficiência menor que a topologia HB. 
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RESUMO: A iluminação do setor industrial 
apresenta um diferencial, pois o consumo de 
energia é de quarenta por cento. Com a troca 
das luminárias convencionais por luminárias LED 
foi possível verificar uma redução de 59,4% do 
consumo de energia e 21,8% da quantidade de 
luminárias, além de um PayBack inferior a 5 anos 
sobre o investimento. Por meio disto a iluminação 
Led apresenta tanto a viabilidade de implantação 
como benefícios para o setor industrial. 
PALAVRAS-CHAVE: LED, Industrial, Iluminação, 
Viabilidade. 


FEASIBILITY OF REPLACING THE 
CONVENTIONAL LUMINAIRES BY 
LED LUMINAIRES IN THE INDUSTRIAL 
SECTOR 


ABSTRACT: The lighting of the Industrial sector 
has a differential since the energy consumption is 
of forty percent. With the exchange conventional 
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SETOR INDUSTRIAL 


lamps LED lamps, it was possible to verify a 
reduction of 59.4% of energy consumption and 
21.8% of the number of installations besides 
a PayBack lower than 5 years on investment. 
Therefore, the LED lighting provides both the 
installation feasibility and industrial sector 
benefits. 

KEYWORDS: LED, Industrial, Lighting, Viability. 


11 INTRODUÇÃO 

O princípio do surgimento da iluminação 
inicia no período neolítico, há 7 mil anos atrás 
quando o Homo Erectus descobriu que ao 
esfregar dois gravetos ou friccionar duas pedras 
conseguia-se faíscas que geravam fogo, que 
por sua vez é fonte de calor e luz. Com o tempo 
o homem começou a aprimorar novas formas 
de produzir iluminação sem a necessidade de 
gerar fogo(Tudo, 2018) [6]. 

O início de produzir luz se iniciou no 
ano de 1580 por meio da utilização de velas, 
que eram utilizadas para iluminar teatros 
colocando em castiçais. A partir do ano de 
1780 foi criado a primeira lâmpada a óleo, 
desenvolvida pelo químico suíço Aime Argant. 
Com o desenvolvimento da lâmpada a óleo foi 
desenvolvida a lâmpada a gás, que por meio 
desta foi criado uma mesa de controle do gás 
que controlava a luminosidade que a lâmpada 
emitia, em paralelo com a criação da lâmpada a 
gás, Thomas Drummond entre os anos de 1816 
a 1890 desenvolveu o Lime Light que a emissão 
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da luminosidade era proveniente do cálcio, que por sua vez começavam a dar lugar as 
lâmpadas de arcos voltaicos. Com o desenvolvimento das lâmpadas de arco voltaico, 
criou-se um gerador dínamo, para realizar as ligações das lâmpadas, apresentando um 
brilho mais intenso e uma eficiência melhor. As modernas lâmpadas elétricas começaram 
a ganhar atenção no início do século XIX, onde Thomas Edison aperfeiçoa a primeira 
lâmpada incandescente. Com este passo inicial, entre os anos de 1906 a 1933 houve 
uma grande evolução das lâmpadas elétricas, que deram início a construção dos primeiros 
refletores e lâmpadas de alta potência (Iluminação, 2018) [3]. 

Com todo o desenvolvimento em relação a iluminação, surgiu em meados de 1907 
uma descoberta pelo inglês Henry Joseph Round que matérias inorgânicas podem emitir 
luminosidade ao aplicar corrente elétrica, entretanto a descoberta foi esquecida devido o 
foco do desenvolvimento ser o novo sistema de orientação para o transporte marítimo dentre 
os anos de 1921 até 1942. A partir do ano de 2006 surge o primeiro diodo com capacidade 
de gerar 100 lumens por watts, contudo apresenta algumas restrições. Com a pesquisa 
voltada para a iluminação LED em alta foi produzido diodos com eficiência gigantesca de 
250 lumens por watts, com evolução da iluminação Led, as luminárias apresentam maior 
sustentabilidade, economia de consumo e vida útil maior que a convencional, sendo assim 
é a iluminação do futuro (Ledvance, 2018) [4]. 

As indústrias atualmente buscam gerar mais sustentabilidade, por meio desta visão, 
a implantação de luminárias LEDs e uma solução, portanto, para que as empresas possam 
investir é necessário um estudo de viabilidade minucioso, para demostrar a economia e os 
benefícios do Led na indústria. O objetivo do estudo é apresentar um comparativo entre a 
iluminação convencional pelo Led e apresentar um PayBack(retorno de investimento) de 
acordo com investimento de implantação das luminárias LEDs, além dos benefícios que o 
Led proporcionará a Indústria. 


21 FUNDAMENTOS PARA A INVESTIGAÇÃO 


Nesta seção, aspectos de dimensionamento de luminárias e o orçamento para 
investimento são considerados. Ademais, as fundamentações para cálculos, a estratégia 
de condução do trabalho, bem como a forma de análise dos resultados são apresentadas. 


2.1 Iluminação industrial 


A iluminação industrial é um dos fatores aos quais muitas indústrias acabam 
esquecendo como influenciadores para economia e sustentabilidade. A indústria 
com iluminação adequada apresenta um diferencial, devido a melhor visibilidade nos 
setores. Os benefícios que a iluminação ideal pode gerar são a economia e prevenção 
com acidentes por luminosidade insuficiente. Com o intuito de melhorar a eficiência, as 


luminárias dimensionadas de forma correta geram estes benefícios. Por meio disto podem 


Projeto, Análise e Otimização na Área das Engenharias Capítulo 3 é 


ser realizados estudos para apresentar a melhor estratégia. 


2.2 Estrutura para a investigação 


Para apresentarmos uma estrutura de investigação devemos levar em consideração 


algumas informações que ajudaram na investigação. 


2.2.1 Horas de consumo 


De acordo com a Lei nº 8.631 e o Decreto nº 774 foi criado um sistema de tarifas que 
define a Hora de Ponta(HP) com duração de três horas consecutivas dentro do período das 
17 às 22 horas, onde depende de cada concessionária, as HP foram estabelecidas devido 
o auto consumo de energia em horas específicos do dia, com o intuito de reduzir este 
consumo. A tarifa de energia e maior, podendo chegar 3 vezes mais que a HFP. No estudo 
apresentado os horários definidos são das 18 às 21 horas como HP. As horas restantes são 
consideradas como Horas Fora de Ponta(HFP) onde o consumo de energia é distribuído de 
forma a não sobrecarregar o Sistema Interligado Nacional(SIN).(Franco,1993) [2]. 


2.2.2 Contrato de venda de energia 


Com a cogeração das termoelétricas por biomassa do bagaço de cana de açúcar 
é necessário um contrato para venda de energia para a concessionária. De acordo 
com a Lei nº 13.299, de 21 de junho de 2016, dispõe valores do custo de acordo com a 
capacidade de geração da usina. O valor do contrato é utilizado para realização do cálculo 
de safra( Temer,2016) [5]. 


2.2.3 Safra e entressafra de uma indústria 


Em usinas termoelétricas por biomassa da cana de açúcar, apresenta dois períodos: 
Safra e Entressafra. Os períodos de Safras são onde ocorre todo o processo de fabricação 
de produtos e subprodutos além da cogeração de energia elétrica. A safra costuma 
corresponder a maior parte do tempo de funcionamento da termoelétrica, com isso durante 
esse período e utilizado para fins de cálculo, como o período que a empresa poderia vender 
a energia que está sendo gasta. Para os cálculos de Entressafra e considerado as HP e 
HFP, pois neste período a empresa não gera energia para suprir sua necessidade, fazendo 
com que tenha que comprar da concessionaria. Neste período e feito a manutenção das 
indústrias. 


2.2.4 Reatores das luminárias 


As luminárias convencionais necessitam de reatores para seu funcionamento, 
porem na maioria das vezes eles são desprezíveis nos cálculos devido sua perda se baixa, 
entretanto, os reatores com o passar do tempo geram perdas maiores que o fabricado. 
Por meio disto os reatores devem ser levados em consideração pois apresentam perdas 
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que com o tempo são significativas. Principalmente no setor industrial onde as luminárias 


possuem altas potências e geram perdas significativas. 


2.2.5 Luminárias leds 


As luminárias LEDs são constituídas por Diodos Emissores de Luz, que apresentam 


melhor eficiência, robustez e alta duração de funcionamento. 


2.3 Cálculos específicos para levantamento de dados 


Para que a investigação seja realizada algumas fórmulas podem ser utilizadas. 


2.3.1 Potência total 


Para iniciar os estudos e necessário realizar o cálculo da potência instalada e 


posteriormente da potência a ser instalada. Para isso utilizamos a seguinte formula. 
Pr — NpxP (1) 


Onde: 

P. — Potência Total 

N, — Quantidade de Luminárias 
P, — Potência das Luminárias 


2.3.2 Consumo de energia total 


O consumo de energia e calculado em função da potência instalada e das horas 
de funcionamento. Para os consumos de energia durante a entressafra e utilizado as 
seguintes fórmulas; 


Crçmp) = PrxHyp (2) 


Onde: 
Cr) > Consumo de Energia Horas de Ponta 
P. — Potência Total das Luminárias 


H,p —> Horas de Ponta de Mensal 


Cromrp) = PrxHgrp (3) 


Onde: 
Cure», > Consumo de Energia Horas Fora de Ponta 
P. — Potência Total das Luminárias 


H,» —> Horas Fora de Ponta Mensal 
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O consumo de energia durante a safra e realizado pela soma dos consumos em HP 
e HFP. 


Cr = CrçHp)XCr(HrP) (4) 


Onde: 
C. —> Consumo de Energia Total 
Cr) > Consumo de Energia Horas de Ponta 


C. ur», —> Consumo de Energia Horas Fora de Ponta 


2.3.8 Custo de energia(kWh) 


Para realizar o custo de energia e considerado a seguintes fórmulas para a 
entressafra e utilizado: 


Erçnp) = CrmpXPrmp (5) 


Onde: 
Emp) —> Custo de Energia Total em Horas de Ponta 


Cu») > Consumo de Energia em Horas de Ponta 


P, nr) Preço do kWh das Horas de Ponta 


Erçnrp) = Crçmrp)XPrçmep) (6) 


Onde: 
Eur» > Custo de Energia Total em Horas Fora de Ponta 


Cure) > Consumo de Energia em Horas Fora de Ponta 


P ed Preço do kWh das Horas Fora de Ponta 


r(HFP 
Para safra e utilizado o preço que foi definido em contrato com a concessionária. 


Er — CrxPre (7) 


Onde: 

E, — Custo de Energia Total 

C. —> Consumo de Energia Total 

Pr. —> Preço do kWh definido em contrato 


2.3.4 Cálculo de economia 


O cálculo de Economia e feito em função da diferença do levantamento convencional 


versus previsto LED. Para realizar o cálculo utilizamos as seguintes fórmulas: 
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Ecs = Cgr(so) — Ceres) (8) 


Onde: 

Ec; —> Economia Mensal Safra 

Cgrsc > Custo Geral Total Safra Convencional 
Cgs, > Custo Geral Total Safra LED 


Ecgs = Cgreesc) — Ceres) (9) 


Onde: 

Ec.s —> Economia Mensal Entressafra 

Cgr(esç > Custo Geral Total Entressafra Convencional 
Cgr(es > Custo Geral Total Entressafra LED 


A Economia anual utiliza a seguinte fórmula: 
Ecr = Mes * Elgg — Mg * ECs (10) 


Onde: 

Ec, —> Economia Anual 

Mes — Quantidade em Meses de Entressafra 
Ec.s —> Economia Mensal Entressafra 

Ms — Quantidade em Meses de Safra 

Ec; —> Economia Mensal Safra 


2.3.5 Cálculo de viabilidade 


O cálculo de viabilidade e realizado utilizando a seguinte formula: 
EE = Ir/Ecr (11) 


Onde: 
PB — Retorno de Investimento em Anos 
|. — Investimento Total 


EC. —> Economia Anual 


2.4 Análise das informações 


O método de análise dos dados para verificar a viabilidade se dá por meio do cálculo 
da fórmula 11, onde se o resultado obtido estiver abaixo de 5 o investimento e eficiente, 
entretanto, se o resultado for maior que 5 e menor que 10 o retorno apresentará viabilidade 


porém com um retorno menor que o esperado, e se o resultado for maior que 10 o projeto 
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não apresentará viabilidade. 


Como o intuito de seguir os passos segue um esquema que deve ser seguido. 


Levantamento da Projeto 
Iluminação ———————8" Luminotécnico 
convencional 
PayBack de Projeção de Gasto 
dam 4——————— na 
Viabilidade com Iluminação LED 


31 RESULTADOS 


Os dados coletados foram de uma indústria da região do triângulo mineiro, com 
apuração das informações de safra e entressafra 2017-2018 com duração de um ano. Os 
dados compilados são divididos em 4 etapas: 1º etapa — Levantamento de Informações 
e Despesas com Iluminação Convencional; 2º etapa — Projeto Luminotécnico; 3º etapa — 
Projeção de Despesas com Iluminação LED; 4º etapa — Viabilidade de PayBack. 


3.1 Etapa 1: Levantamento de informações e despesas com iluminação 
convencional 


Nesta parte, e apresentado a potência instalada no setor da indústria ao qual foi 
escolhido. Na Tabela 1 os dados são apresentados. 


Descrição o Peas Ca Potência Total (kW) 
Luminaria Vapor Metálico 250W 28 270 7,06 
Luminaria Fluorescente Tubular 20W 5 43,2 0,216 
Refletor Vapor Metálico 400W 8 432 3,456 
Luminaria Fluorescente Tubular 40W 16 82 19 12 
Luminarias Vapor de Sódio 400W 144 432 62,208 
Total 201 74,752 


Tabela 1- Potência instalada no setor. 


Com as informações apresentadas acima e possível calcular o consumo de energia 
mensal, considerando Horas de Ponta e Horas Fora de Ponta, para Entressafra. Na Tabela 
2 e Tabela 3 e possível verificar o proposto: 
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Consumo de 


Descrição Potência Total (kW) Hora de Ponta (horas) Eneraia (kWh 
Luminaria Vapor Metálico 250W 7,56 66 498,96 
Luminaria Fluorescente Tubular 20W 0,216 66 14,256 
Refletor Vapor Metálico 400W 3,456 66 228,096 
Luminaria Fluorescente Tubular 40W 1,312 66 86,592 
Luminarias Vapor de Sódio 400W 62,208 66 4105,/28 
Total 74,752 4933,632 


Tabela 2 - Consumo de energia elétrica mensal em Hora de Ponta. 


Consumo de 


Descrição Potência Total (KW) Hora Fora de Ponta (horas) Energia (kWih| 
Luminaria Vapor Metálico 250W 7,56 294 2222,64 
Luminaria Fluorescente Tubular 20W 0,216 294 63,504 
Refletor Vapor Metálico 400W 3,456 294 1016,064 
Luminaria Fluorescente Tubular 40W 1,312 294 385,/28 
Luminarias Vapor de Sódio 400W 62,208 654 40684,032 
Total 74,752 44371,968 


Tabela 3 - Consumo de energia elétrica mensal em Hora Fora de Ponta. 


Com os dados de Entressafra é realizado o cálculo de energia gasta durante a 
entressafra. Considerando o preço médio para HP e HFP como mostra a Tabela 4 a seguir. 


Consumo Total de Preço Médio Duração Gasto de Energia (R$) 


Energia (kW/h) do KWih (R$) (Meses) 


Horas de Ponta 4933,632 0,49606647 3,50 R$ 8.565,93 
Horas Fora de Ponta 44371,968 0,338317651 3,50 R$ 52.541,37 
Total 49305,6 R$ 61.107,30 


Tabela 4 - Gasto de energia elétrica no período de entressafra. 


Com base no consumo mensal(HP+HFP) e possível calcular o gasto de energia 
durante a safra. Considerando o preço do contrato de energia como mostra a Tabela 5 a 


seguir. 


Consumo Total de Preço do Duração Gasto de Energia (R$) 


Energia (MW/h) MWwih (R$) (Meses) 
Safra 49,3056 200 8,50 R$ 83.819,52 


Tabela 5 - Gasto de energia elétrica no período de safra. 
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A próxima compilação das informações são os gasto com manutenção, materiais e 
descarte de lâmpadas. Na Tabela 6 a seguir se os dados compilados durante a safra. 


Custo (R$) duração (Mensal) Gasto Total (R$) 


Manutenção 21,09 8,5 R$ 179,27 
Materiais 210,81 8,5 R$ 1.791,89 
Descarte de Lâmpadas 1,22 8,5 R$ 10,37 

Total R$ 1.981,52 


Tabela 6 - Gasto com manutenção, materiais e descarte de luminarias. 


3.2 Etapa 2: Projeto luminotécnico 


O levantamento do projeto luminotécnico para aplicação LED apresentou as 
seguintes luminárias para substituição. Na Tabela 7 a seguir estão os dados da iluminação 


considerados de acordo com o projeto do setor de estudo. 


o Descrição 
Luminaria Industrial LED Arandela 30º - 94W 

Luminaria Industrial LED Arandela 30º - 31W 

Projetor LED 207W 

Luminaria Industrial LED Plafonier - 31W 

Luminaria Industrial LED Suporte para poste 30º - 31W 
Luminaria Comercial LED 36W 

Projetor LED-EX 48W 

Projetor LED 270W 


Tabela 7 - Potência instalada no setor. 


Por meio destas luminárias que foram definidas de acordo com o projeto 


luminotécnico, podemos fazer uma projeção de custos com a implantação LED. 


3.3 Etapa 3: Projeção de despesas com iluminação led 


A etapa de realização e levantar uma projeção de gastos futuros com as luminárias 
LED implantadas. Por meio disto a Tabela 8 a seguir apresenta a potência a ser instalada. 
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Quantidade de Potência da Potência Total 


Pescnçae Luminarias | Luminaria (W) (kW) 

Luminaria Industrial LED Arandela 30º - 94W 17 94 1,598 
Luminaria Industrial LED Arandela 30º - 31W 9 31 0,279 
Projetor LED 207W 8 207 1,656 
Luminaria Industrial LED Plafonier - 31W 3 31 0,093 
Luminaria Industrial LED Suporte para poste 30º - 31W 1 31 0,031 
Luminaria Comercial LED 36W 13 36 0,468 

Projetor LED-EX 48W 10 48 0,48 
Projetor LED 270W 96 270 25,92 
Total 157 30,525 


Tabela 8 - Projeção de Potência LED instalada no setor. 


Com os cálculos de potência LED instalada podemos realizar os cálculos de 
consumo de energia. Na Tabela 9 e Tabela 10 está expresso o consumo de energia mensal 
em Horas de Ponta e Horas Fora de Ponta. 


ae Potência Total Hora de Ponta Consumo de 
Descrição 


O  W) (horas) Energia (kWih) 
Luminaria Industrial LED Arandela 30º - 94W 1,598 66 105,468 
Luminaria Industrial LED Arandela 30º - 31W 0,279 66 18,414 
Projetor LED 207W 1,656 66 109,296 
Luminaria Industrial LED Plafonier - 31W 0,093 66 6,138 
Luminaria Industrial LED Suporte para poste 30º - 31W 0,031 66 2,046 
Luminaria Comercial LED 36W 0,468 66 30,888 
Projetor LED-EX 48W 0,48 66 31,08 
Projetor LED 270W 25,92 66 1710,72 

Total 30,525 2014,65 


Tabela 9 - Consumo de energia elétrica mensal em Hora de Ponta LED. 


di Potência Total HoraFora de Consumo de 
Descrição 


kW Ponta (horas) Energia (kW/h 

Luminaria Industrial LED Arandela 30º - 94W 1,598 294 469,812 
Luminaria Industrial LED Arandela 30º - 31W 0,279 294 82,026 
Projetor LED 207W 1,656 294 486,864 
Luminaria Industrial LED Plafonier - 31W 0,093 294 27,342 
Luminaria Industrial LED Suporte para poste 30º - 31W 0,031 294 9,114 
Luminaria Comercial LED 36W 0,468 294 137,592 
Projetor LED-EX 48W 0,48 294 141,12 
Projetor LED 270W 25,92 654 16951,68 

Total 30,525 18305,55 


Tabela 10 - Consumo de energia elétrica mensal em Hora Fora de Ponta LED. 


Considerando os mesmos preços médios de HP e HFP para entressafra, na Tabela 
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11 e apresentado os gastos. 


Consumo Total de Preço Médio Duração Gasto de Energia (R$) 


Energia (kW/h) do KWih (R$) (Meses) 


Horas de Ponta 2014,65 0,49606647 3,50 R$ 3.497,90 
Horas Fora de Ponta 18305,55 0,338317651 3,50 R$ 21.675,82 
Total 20320,2 R$ 25.173,72 


Tabela 11 - Gasto de energia elétrica LED no período de entressafra. 


A seguir na Tabela 12 e apresentado os gastos com energia durante a safra 


considerando o prelo de contrato com a concessionaria. 


Consumo Total de Preço do Duração E 
Energia (MWih)  MWIh(R$) (Meses) Gasto de Energia (R$) 
Safra 20,3202 200 8,50 R$ 34.544,34 


Tabela 12 - Gasto de energia elétrica LED no período de safra. 


Os custos com manutenção, matérias e descarte para a implantação foram 
desconsiderados, devido a cotação para a realização do PayBack se tornou desprezível 
devido a durabilidade da luminária e benefícios oferecidos pela empresa a qual foi feito 


levantamento para apresentar viabilidade da implantação. 


3.4 Etapa 4: Viabilidade de payback 


Com todos os dadas as informações compiladas podemos realizar os cálculos de 
viabilidade de implantação. Na Tabela 13 a seguir e apresentado uma comparação entre a 


iluminação convencional e LED. 


Descrição luminação Iluminação 
Convencional LED 
Quantidade de luminarias 201 157 
Potência Instalada (kW) 14,752 30,925 
Consumo de Energia (kW/h) 49305,6 20320,2 


Tabela 13 - Comparação de iluminação Convencional X LED. 


Na Tabela 13 podemos comprovar que houve uma redução de aproximadamente 
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59,2% da potência instalada, além de uma redução de aproximadamente 21,9% da 
quantidade de luminárias instaladas. Na Tabela 14 e apresentado a comparação de custos 
entre iluminação convencional e LED. 


Iluminação Iluminação 
Convencional LED 
Entressafra 
Gasto de Energia Horas de Ponta R$ 8.565,93 R$ 3.497,90 
Gasto de Energia Horas Fora de Ponta R$ 52.541,37 R$21.675,82 
Safra 

Gasto de Energia por contrato R$ 83.819,52 R$ 34.544,34 
Gasto com Manutenção, Materiale Descarte R$ 1.981,52 R$ - 

Total R$ 146.908,34 R$ 59.718,06 


Tabela 14 - Comparação de custos de iluminação Convencional X LED. 


Na Tabela 14 e apresentado uma economia anual de aproximadamente 59,4%. 

Para apresentar a viabilidade do estudo em questão, foi levantando um orçamento 
de implantação, onde foram considerados as luminárias, os custos com maquinários e mão 
de obra. Na Tabela 15 está apresentado os custos da implantação. 


Custo (R$) 


Orçamento das Luminarias LED R$ 191.756,36 
Custo de Maquinarios R$ 30.000,00 
Mão de Obra R$ 40.000,00 

Total R$ 261.756,36 


Tabela 15 - Investimento para implantação LED. 


Por meio dos valores em questão e realizado os cálculos de PayBack. Na Gráfico 
1 será apresentado os cálculos, onde como comprovado na Tabela 14 a economia anual 
e de R$ 87.190,28 (oitenta e sete mil cento e noventa reais e vinte e oito centavos) será 
considerado como economia anual. 
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R$ 200.000,00 R$ 174.205,04 


R$ 150.000,00 
R$ 100.000,00 R$ 87.010,76 


R$ 50.000,00 
3 


-R$ 182,52 


R$ - 
“R$ 50.000,00 
“R$ 100.000,00 -R$ 87.374,80 


-R$ 150.000,00 





Gráfico 1 — PayBack em 5 anos. 
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Para a atender o estudo de viabilidade os métodos são realizados de forma simples, 
que atende a todos os requisitos necessários para apresentar dados concisos. de acordo 
com Calegari, 2017 apresenta metodologia simples para realizar os cálculos de viabilidade. 
Além de apresentar que o setor industrial apresenta uma porcentagem relevante do 
consumo de energia elétrica, considerando que a iluminação e dos fatores aos quais são 
deixados de lados, por ser simples e desprezíveis, porém podem apresentar um gasto 
expressivo devido a quantidade. 

Por meio dos métodos apresentados na tese os resultados podem ser validados, 
mostrando que seguindo todos os passos e realizando os cálculos corretamente. Podemos 
verificar que os cálculos apresentaram gastos razoáveis devido ser apenas da iluminação 
e de apenas um setor da industrial. Ao analisar os dados e possível verificar que com 
a implantação LED e houve uma redução de aproximadamente 59% do consumo de 
energia elétrica, além de apresentar uma redução aproximadamente 22% da quantidade 
de luminárias. Por meio disto houve uma economia anual de aproximadamente 87 mil reais 
que e significativa. 

A viabilidade foi possível devido os seguintes fatores. A empresa de levantamento 
do investimento apresentou benefícios onde abaixo de 5 anos os custos com troca de 
matérias e manutenção são depressíveis, além de produtos de alta durabilidade e 
qualidade. Podemos observar que com a implantação LED os benefícios gerados para a 
indústria são: Iluminação adequada e eficiente; redução de acidentes devido a iluminação 
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ineficiente; redução de mão de obra com manutenção; Redução com Gasto de Energia e 
um espelho significativo no aumento da produção, onde um lugar bem iluminado gerar uma 


percepção melhor. 


5 | CONCLUSÃO 


A viabilidade de implantação de iluminação LED em um setor industrial se mostrou 


viável com um retorno de investimento de aproximadamente 3 anos após a implantação. 


AGRADECIMENTOS 


O autor agradece o apoio recebido de todos os Docentes da Instituição UNIUBE e 


seus colegas de Graduação. 


REFERÊNCIAS 


[1] Calegari, Raphael; Como funciona o Sistemas de Iluminação na Indústria? ;Disponível em: 
<http://bolog.murrelektronik.com.br/sistemas-de-iluminacao-na-industria/ > Acesso em 16 de novembro 
de 2019. 


[2] Franco, Itamar; RESENDE, Eliseu; VASCNSELOS, Paulino Cícero de; CRUSIUS, Yeda Rorato; 
Disponível em: < http:/Awww.planalto.gov.br/ccivil 03/decreto/1990-1994/D0774.htm >. Acesso em 16 
de novembro de 2019. 


[3] Iluminação, Laboratório de; História da Iluminação: Primeiros Instrumentos; Disponível em: <http:// 
www.iar.unicamp.br/lab/luz/ld/slides/slides diciplina ar 210 b/primeiros instrumentos.pdf>. Acesso em 
07 de outubro de 2019. 


[4] Ledvance; A história do LED; Disponível em: <https://www.ledvance.pt/produtos/conhecimentos- 
sobre-o-produto/nocoes-basicas-do-led/historia-do-led/index.jsp>. Acesso em 06 de outubro de 2019. 


[5] Temer, Michel; Meirelles, Henrique; Filho, Fernando Coelho; Oliveira, Dyogo Henrique de; Osório, 
Fábio Medina ; Disponível em: < http://www.planalto.gov.br/ccivil 03/ Ato2015-2018/2016/Lei/L13299. 
htm >. Acesso em 16 de novembro de 2019. 


[6] Tudo, 10 em; O descobrimento do fogo; Disponível em: <https:/Mww.1Oemtudo.com.br/artigo/o- 
descobrimento-do-fogo/>. Acesso em 21 de novembro de 2018. 


Projeto, Análise e Otimização na Área das Engenharias Capítulo 3 EEN 


CAPÍTULO 4 


PROJETO PARA ELABORAÇÃO DE UMA PEN 


Data de aceite: 04/01/2021 


Rafael Ferreira da Silva 


Centro Universitário União das Américas — 
Uniamérica 

Foz do Iguaçu, Paraná 
http:/lattes.cnpqg.br/7898337334370742 


Welton Abreu Rosa 


Centro Universitário União das Américas — 
Uniamérica 

Foz do Iguaçu, Paraná 
http:/lattes.cnpq.br/6057660697254416 


Luciana Paro Scarin Freitas 


Centro Universitário Comunitário União das 
Américas — UniAmérica 

Foz do Iguaçu-PR 
http:/lattes.cnpq.br/2554057558293125 


Jorge Luis Ribeiro dos Santos Júnior 


Centro Universitário Comunitário União das 
Américas — UniAmérica 

Foz do Iguaçu-PR 
http:/lattes.cnpg.br/0033807609839899 


Luís Henrique Chouay Dall' Agnese 


Centro Universitário União das Américas — 
Uniamérica 

Foz do Iguaçu, Paraná 
http://lattes.cnpg.br/1502220513099067 


Grégori da Cruz Balestra 

Centro Universitário Comunitário União das 
Américas — UniAmérica 

Foz do Iguaçu-PR 
http:/lattes.cnpq.br/8659424533812354 


Projeto, Análise e Otimização na Área das Engenharias 


PLOTTER 


RESUMO: Quando pensamos em precisão e 
agilidade, lembramos de máquinas controladas 
por computadores, essas máquinas denominadas 
CNC podem executar trabalhos com alta 
precisão, a PenPlotter uma impressora CNC é 
capaz de escrever qual quer tipo de caligrafia 
e a plotagem de desenhos. Ela é movida por 2 
nemas instalados com uma polia dentada e um 
servo motor para movimento da caneta, através 
dos vetores gerados pelo software inkscap a 
impressora realiza os movimentos necessários 
para a plotagem. No mercado atual existem 
diversos tipos de PenPlotter com uma variedade 
de modelos, materiais e softwares, através 
de uma seleção escolhemos o que melhor 
se enquadra em nosso projeto. Esse projeto 
apresenta a construção e o funcionamento de 
uma PenPloiter. 

PALAVRAS-CHAVE: PenPlotter; CNC; Precisão. 


ABSTRACT: When we think of precision and 
agility, we remember machines controlled by 
computers, these machines called CNC can 
perform jobs with high precision, the PenPlotter 
is a CNC printer capable of writing any type of 
calligraphy and plotting drawings. Itis moved by 
2 nemas installed with a toothed pulley and a 
servo motor for movement of the pen, through the 
vectors generated by the inkscap software the 
printer performs the necessary movements for the 
plotting. In the current market there are several 
types of PenPlotter with a variety of models, 
materials and software, through a selection we 
choose the one that best fits our project. This 
project presents the construction and operation 
of a PenPlotiter. 
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11 INTRODUÇÃO 


No mundo atual globalizado, o requisito básico de qualquer indústria é produzir 
grandes quantidades de produtos de qualidade com baixa perca de material, baixo custo 
de instalação, com alto acabamento superficial e grande precisão dimensional. Isto pode 
ser conseguido através de máquinas operadas por computadores. Eles são basicamente 
conhecidos como máquinas CNC. Segundo o site Mecanicalndustrial “As máquinas CNC 
são máquinas automatizadas que funcionam por meio de sistemas CNC, ou seja, por 
Controle Numérico Computadorizado. Isso significa que as máquinas controladas por 
sistemas CNC são muito mais rápidas em processos de produção” 

Utilizando este tipo de máquina, os produtos são produzidos em uma taxa mais 
rápida com alta precisão e menos interferência humana. Os sistemas de CNCs mais 
utilizadas são os de dois eixos e de três eixos. 

Este sistema consiste em três partes importantes: Projeto mecânico, módulos de 


acionamento e softwares de controle. 


* | O projeto mecânico consiste no corpo estrutural; 
* | Os módulos de unidade consistem no microprocessador; 


* Fo Software do Sistema é usado para gerar o desenho na folha. 


A PenPlotter é a automação de máquinas que são operadas por comandos 
precisamente programados. O principal objetivo de uma Plotter CNC está voltado para a 
leitura de vários vetores de produtos. O princípio de funcionamento é muito semelhante à 
de uma máquina CNC. Nesse sistema, o operador humano não executa manualmente o 
processo de desenhar os vetores, a máquina, portando uma caneta e um servo motor para 
acionamento em seu eixo Z, através de um controle programado, produz um trabalho de 


alta qualidade e maior agilidade em comparação com o trabalho humano. 


2 | MATERIAIS E MÉTODOS 


Os dados foram coletados através da análise e estudo de artigos científicos 
referentes ao produto, assim como em referências dos fabricantes, onde contém 
informações detalhadas do projeto e datasheet dos motores. Com os dados dos artigos 
disponibilizados, foi possível realizar um comparativo entre vários destes, onde através de 
anotações de cada artigo, possibilitou adquirir informações necessárias para a produção de 


um produto com qualidade e funcionalidade. 
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2.1 Memorial de cálculos 


Para o correto funcionamento dos motores de passo, o software precisa das 


seguintes especificações dos motores: 


Passos que cada eixo realiza: 


Eixo X: $ 100=775.000 (x, step/mm) 
Eixo Y: $ 101=367.000 (y, step/mm) 


Eixo Z: $ 102=250.000 (z, step/mm) 


Aceleração de cada eixo: 


Eixo X: $120=10.000 (x accel, mm/sec"2) 
Eixo Y: $121=10.000 (y accel, mm/sec"2) 


Eixo Z: $122=10.000 (z accel, mm/sec"2) 


Limite máximo de cada eixo: 


Eixo X: $130=400.000 (x max travel,mm) 
Eixo Y: $131=300.000 (y max travel,mm) 


Eixo Z: $132=200.000 (z max travel,mm) 


2.2 Modelagem 


O seguinte projeto é composto por um conjunto de equipamentos mecânicos e 


eletrônicos: 
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A base estrutural desta máquina consiste no posicionamento paralelo de duas 
barras roscadas M10 (650mm de comprimento cada), com porca e contra por- 
ca, exigindo uma quantidade de oito porcas. 


A base de movimentação do projeto é composta por quatro barras lisas de 8 
mm (duas de 650 mm de comprimento e duas de 900 mm de comprimento), 
onde duas destas (900mm) estão posicionadas em sentido ao eixo X e as ou- 
tras duas (650mm), posicionadas em sentido ao eixo Y, suportadas pelos supor- 
tes dos Nemas 17 Steppers e suporte do servo motor. 


O carro de movimentação foi projetado e impresso em uma impressora 3D com 
filamento ABS, para comportar as travas das correias GT2, assim como oito 
LM8UU (rolamentos lineares) e servir como guia para as posições do eixo Z. 


Para o mecanismo de movimentação dos eixos X e Y, foram utilizados 2 nemas 
17 Steppers instalados com uma polia dentada, GT2 com 20 dentes, em cada 
eixo do motor, onde para a movimentação do carro, elas são interligadas por 
correias dentadas, fixadas no carro central. 


* | O sistema de impressão é constituído por dois tubos de 4mm, com 100mm de 
comprimento, fixados ao suporte para um micro servo SG90 (mais um extensor 
de cabo de 250 mm), onde ao eixo do micro servo está instalado com um su- 
porte ajustavel para canetas. 


* | O controle eletrônico escolhido para os componentes foi o Arduino UNO, onde 
o mesmo recebe o código G do computador, converte em informações para os 
motores e micro servos motores, utilizando uma CNCshield. Para os motores 
Nemas 17, utilizase dois Drivers A4988 conectados a CNCshield, e para o micro 
servo, utiliza-se uma porta padrão do Arduino. 


* | Aalimentação do circuito do projeto se baseia em um acionamento por fonte de 
alimentação 12v — 2A. 


31 RESULTADOS E DISCUSSÃO 


O sistema trabalha basicamente com dois motores de passo (um para o eixo X e um 
para o eixo Y) e micro-servo-controlador (para o eixo Z). Onde o circuito Arduino recebe 
a programação para a funcionalidade de dimensionamento da área de atuação do plotter 
e também o comando de atuação dos motores. Através do software inkscap é colocado 
um desenho ou até mesmo textos e é gerado o G-CODE que através do programa grbl 
o g-code que avia sido criado se transforma em coordenadas para os eixos e assim são 
plotados na folha que está colocado na mesa de desenho usando microcontrolador. Este 
sistema reduz o esforço humano e também reduz as chances de erro. 

Após uma detalhada análise de produtos similares existentes no mercado chegamos 
a conclusão que esse modelo seria a melhor opção para nosso projeto, por ter uma fácil 
construção e montagem de suas correias, como vemos na figura 1 abaixo: 





Figura 1 - Desenho da Pen Plotter no SolidWorks. 


Fonte: Autoria própria 
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Esse modelo se destaca por usar apenas 1 motor de passo no eixo X e um Motor no 
Eixo Y, Diferente do esboço mostrado abaixo Feito por Smokevicz, Fabio; Araújo, Ricardo. 
Desenvolvimento de uma Plotter a caneta de baixo custo para escrita e desenho em 
superfícies planas - Joinville: Instituto Federal de Santa Catarina, 2017. 





Figura 2 - Esboço da Plotter em SolidWorks 


Fonte: Autoria própria 


Abaixo podemos conferir nas Figuras 3 e 4 imagens reais do projeto após sua 
montagem. 








Figura 3 — Visão Lateral PenPlotter 


Fonte: Autoria própria 
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Figura 4 — Visão Geral PenPlotter 


Fonte: Autoria própria 


41 CONCLUSÃO 


Com a realização desse trabalho concluímos que com o objetivo de projetar e 
fabricar um projeto com funcionalidades precisas, exige o dimensionamento dos passos 
dos motores nemas, assim como alinhamento do código para funcionamento sincronizado 
de todo o sistema. Com a experiência do protótipo foi-se selecionado os programas mais 
adequados para a leitura do código g e acionamento dos motores, sendo assim, carregando 
uma quantidade de conhecimento teórico e prático ao projeto. 
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RESUMO: Este artigo tem o objetivo de 
apresentar uma avaliação das máquinas e 
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equipamentos de uma frota utilizando quatro 
métodos de depreciação, obtendo-se desta 
forma, o valor real de mercado, depreciação 
total da frota, identificação da quantidade 
de ativos sucateados, bem como prever os 
próximos ativos que atingirão o valor residual 
ou de sucata. A proposta foi motivada pelo 
conhecimento específico, em avaliação de 
máquinas e equipamentos, que o engenheiro de 
avaliação deve possuir e pela necessidade de se 
estabelecer o valor real de mercado, de modo 
a auxiliar na tomada de decisão em relação 
a destinação destes bens. Esta proposta se 
deteve nos problemas relacionados a pesquisas 
e métodos inadequados para o estabelecimento 
dos valores reais para se avaliar os bens. A 
comparação destes métodos de depreciação 
pode acarretar obtenção de valores mais 
fidedignos dos ativos. Os métodos empregados 
contemplam a vida útil, vida aparente e valor 
de aquisição dos ativos como parâmetro. Os 
resultados indicam que a comparação entre os 
métodos foi efetiva para o cálculo da depreciação 
e dos valores de mercado, pois se mostraram mais 
adequados para finalidades distintas. Assim, a 
comparação entre estes métodos permite que os 
resultados da avaliação tenham menor ou maior 
decrepitude do valor real de mercado, em que a 
contribuição se dá na obtenção de um resultado 
mais assertivo na tomada de decisão das 
máquinas e equipamentos, por parte de quem os 
detém ou os controla. Neste sentido, o método 
de Cole apresentou o maior valor para o custo 
total da frota podendo ser utilizado para a venda 
dos ativos, proporcionando um maior retorno ao 
proprietário. O método de Caires apresentou o 
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menor valor total da frota, o que seria melhor para avaliações, em que se vai pagar alguma 
tarifa, alíquota, taxa ou afins tomando como base a avaliação. 

PALAVRAS-CHAVE: Métodos de Avaliações. Máquinas e equipamentos. Depreciação. Valor 
residual. 


DEPRECIATION OF MACHINERY AND EQUIPMENT USING, LINE, COLE, 
CONSTANT PERCENTAGE AND CAIRES METHODS 


ABSTRACT: This article presents an evaluation of the machines and equipment of a 
fleet using four methods of depreciation, obtaining in this way the real market value, total 
depreciation of the fleet, identification of the amount of assets scrapped, as well as predicting 
the next assets that will reach the residual or scrap value. The proposal was motivated by the 
specific knowledge in evaluation of machines and equipment that the evaluation engineer 
must possess and by the need to establish the real market value, in order to assist in the 
decision making regarding the destination of these asseis. This proposal focused on problems 
related to research and inadequate methods for establishing real values to evaluate the 
assets. The comparison of these depreciation methods may lead to more reliable values of 
the assets. The methods used include the useful life, apparent life and asset acquisition value 
as a parameter. The results indicate that the comparison between the methods was effective 
for the calculation of depreciation and market values, since they proved to be more suitable 
for different purposes. Thus, the comparison between these methods allows the results 
of the evaluation to have a lower or greater decrepitude of the real market value, in which 
the contribution is made in obtaining a more assertive result in the decision making of the 
machines and equipment, by those who holds or controls them. In this sense, Cole method 
presented the highest value for the total cost of the fleet and could be used to sell the asseis, 
providing a higher return to the owner. The Caires method presented the lowest total value of 
the fleet, which would be better for evaluations, in which some charge, aliquot, fee or the like 
will be paid based on the valuation. 

KEYWORDS: Assessment Methods. Machines and equipment. Depreciation. Residual value. 


11 INTRODUÇÃO 


Os engenheiros de avaliação devem manter-se atualizados no que cerne os preços 
de mercado das máquinas e equipamentos, devido alterações na economia, variação nas 
demandas de mercado, globalização e principalmente pelo avanço tecnológico, que os 
impulsionam a tal atualização (BUSTAMANTE, 2000). 

Os ativos industriais assumem valores diferenciados, porque quem define esta 
valoração é a finalidade da avaliação, atrelado atualização anteriormente citada do 
engenheiro avaliador é exigida grande experiência, como em avaliações de complexos 
indústrias por exemplo (OLIVEIRA, 2007). 

Engenheiros de avaliação vêm se dedicando divulgação, pesquisa e estudo de 
técnicas de avaliação, posto que a necessidade de conhecimento do valor de um bem deve 
ser embasada em preceitos legais e fundamento científico (CAIRES, 1972). 
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A subjetividade dos métodos utilizados, na avaliação de máquinas e equipamentos, e 
a falta de parâmetros para realizar a depreciação em muitos casos não possibilitam um valor 
mais próximo do real com o valor de mercado. A comparação dos métodos proposto neste 
artigo, visa fornecer uma gama de resultados, de sorte a permitir um método compatível 
com a finalidade da avaliação, convergindo nestes meios, um valor mais assertivo. 

Nesta linha de raciocínio, pretendeu-se identificar a importância e quantificar o 
impacto da comparação dos métodos aplicados para obtenção do valor de mercado dos 
ativos. 

Ao final, foram calculados os valores de mercado e depreciação de cada ativo, bem 
como identificado quantos foram sucateados e os próximos que atingirão o valor residual, 
apresentando qual metodologia está mais fidedigna com o valor de mercado em relação a 
uma avaliação específica. 

Até que se atinja sua vida útil, ou valor de sucata, o método da linha reta demonstra 
uma depreciação continuamente decrescente. Nesta etapa, o termo “valor de sucata”, o 
qual é o valor de mercado reaproveitável de um bem, compete a ABNT 14653 a definição, 
na condição de desativação. Isto é, o valor que um ferro-velho pagaria por um ativo móvel 
(UBIRAJARA, 2013). 

Para Kuhnem (2001), a fórmula para o cálculo da depreciação linear da pela equação 
(1) abaixo. 


—C(Vn-Vr) 
e 


DL (1) 


Onde DL é a depreciação linear anual, Vn é a o valor do ativo novo, Vr o valor 
residual e n a vida útil. Da equação depreende-se que o valor da depreciação a cada ano é 
subtraído sempre a mesma magnitude, daí resultando numa linha reta decrescente. 

Conforme Assis (2011), o método da soma de dígitos ou de Cole considera que o 
valor de um bem ou serviço decresce a uma taxa decrescente, sendo o valor da depreciação 
no ano t. À equação (2) é usada para o cálculo da depreciação de Cole: 

(2.(N-t+1)) 


DESA CARNE (2) 


Onde DC é a depreciação Cole anual, Vn é a o valor do ativo novo, Vr o valor 
residual, N a vida útil e tvida aparente ou idade atual da máquina ou equipamento. 

Consoante Wammws (2013), diz que o método da porcentagem constante é também 
conhecido como método de Matheson ou Exponencial, ou ainda método da porcentagem 
fixa sobre o valor contábil. Este método necessita de incorporar o valor residual no cálculo da 
depreciação. O valor da depreciação, à medida que tende a obsolescência, vai diminuindo, 
sendo sempre maior no início da vida útil. Mostra-se eficiente em demonstrar o declínio do 
valor de mercado do ativo e determinação do valor residual. A equação (3) para obtenção 


do valor da depreciação é dada por: 
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N|Vr 
T= |— (3) 
Vn 
Onde N representa a vida útil, Vro valor residual e Vn o valor do ativo novo. 
Segundo Oliveira (2007), O método de depreciação utilizado por Hélio Roberto 
Caires, considera que a depreciação é uma função dependente de idade (t), práticas de 


manutenção (4), regime de trabalho (2), vida útil esperada (1) de acordo com a equação (4): 


A 


D L ' , — 
Ep (ou. Cc. (5) 


l+B.e 


(4) 


Nesta equação a maioria das incógnitas são dadas: A é igual a 1,347961431; B 
representa uma outra expressão cujo o termo é (A — 1); C é 3,579761431; “e” é o número 
de Euler, denominado em homenagem ao matemático suíço Leonhard Euler, é a base 
dos logaritmos naturais, e vale 2,7182; e finalmente a função de desgaste a qual está em 


função de Te |, que expressam carga e manutenção na equação (5) 


du, T) — 0,853081710 g9,067348748.7-0,041679277.4—0,001022860 TU (5) 


Com função de levar em consideração as condições de trabalho das máquinas e 
equipamentos, o fator de carga “T” é utilizado. Podendo representar se o ativo está sob 
condições de trabalho nulo, leve, normal, pesado ou extremo, assumindo os valores 0, 5 
10, 15 e 20 respectivamente. 

Um outro fator de notável importância é o fator de práticas de manutenção “u”, que 
estabelece o nível de manutenção que o ativo possui. Este fator classifica a prática de 
manutenção em inexistente, deficiente, normal, rigorosa ou perfeita, relacionando com os 


seguintes valores 0, 5 10, 15 e 20 respectivamente. 


2 | MATERIAIS E MÉTODOS 


O mercado de ativos usados é muito restrito ou não existe para alguns tipos 
específicos, de sorte que dificultam a realização da avaliação por meio de inferência 
estatística, então para os profissionais legalmente habilitados avaliar máquinas e 
equipamentos, na grande maioria das vezes se valem de métodos de depreciação, a fim de 
obter o valor de mercado destes. 

Seguindo esta concepção, na infinidade de métodos para o cálculo da depreciação 
de máquinas e equipamentos foram escolhidos alguns métodos, por serem consagrados, 
difundidos e bastantes utilizados na avaliação de bens móveis, conforme referencial teórico. 
No entanto, a execução na prática denota que estes métodos de depreciação possuem 
algumas defasagens e, em casos particulares, não retornam o real valor de mercado de 


um bem móvel. 
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Na tentativa de compensar o desfalque, foi implementada a utilização conjunta de 
quatro metodologias consagradas e frequentemente utilizadas na área de engenharia de 
avaliações, que são: o método de linha, método linear ou método de quotas constantes; 
método da soma de dígitos ou método e cole; método da percentagem constante, método 
do saldo decrescente, método de Matheson ou Exponencial, ainda método da porcentagem 
fixa sobre o valor contábil; e fisalmente o método de Hélio de Caires. 

A ideia de utilizar, para o cálculo das depreciações, os quatro métodos, é identificar 
qual o valor aparenta ser o mais adequado a ser escolhido como valor de mercado do ativo, 
fazendo um certo ajuste neste valor. 

Para realizar a avaliação neste artigo, foram coletadas informações de sessenta 
e três ativos do setor da construção civil no dia 14 de junho de 2016. Em seguida, foram 
consultados os manuais de operação, para que se encontrassem alguns parâmetros 
importantes requeridos pelos métodos escolhidos e características peculiares que viessem 
a ser Úteis no desenrolar dos cálculos. Dentre as informações coletadas acerca dos ativos 
destacam-se: ano de fabricação, vida útil, vida aparente, valor residual, valores dos bens 
novos, práticas de manutenção, fator trabalho, marcas, categorias, classes, status na obra 
de construção. Na impossibilidade de se encontrar algumas informações recorreu-se a 
internet, como por exemplo, para achar os valores de bens novos. 

Com base nestas informações foi possível utilizar as metodologias propostas neste 
trabalho, permitindo calcular, para as máquinas e equipamentos, o valor de mercado, a 
depreciação da frota, identificar os ativos que atingiram o valor residual e os próximos que 
atingiram este estado em 2018. E em seguida realizar comparações entre as metodologias. 

No trabalho são apresentadas tabelas contendo a listagem de máquinas e 
equipamentos selecionadas, os valores dos bens novos, também são apresentadas as 
tabelas de depreciação das máquinas e equipamentos das quatro metodologias, fluxograma 
das etapas dos procedimentos, gráficos de status dos ativos, gráficos das categorias dos 
ativos, gráfico das marcas dos ativos, gráficos do ano de fabricação dos ativos, e por fim 
são apresentadas gráficos com: os valores de mercado dos ativos, depreciação total da 
frota, quantidade de ativos sucateados e quantos atingiram o valor residual. 


Fluxograma 


Com intuito de tornar a análise do processo muito mais simples, foi esquematizado 
as etapas do trabalho no fluxograma abaixo: 
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Coleta de dados 


Definição dos 
parâmetros para 
utilização dos métodos 


Método da 
Método linear d ds Método de Cole porcentagem er sagens 
ecrescente a o Caires 
Resultados 


Comparação 





Figura 1- Fluxograma de Coleta de Dados. 
Fonte: Própria (2016). 


Primeiramente, efetuou-se a coleta de dados, verificando-se os manuais de 
operação e catálogos e na impossibilidade, ou seja, no extravio ou quando não havia, 
a coleta se deu por meio de inspeção visual, indo até o ativo e checando a placa de 
descrições técnicas contidas no chassi ou em parte específicas da lataria. Em seguida, 
foi montada uma base de dados contento a reunião de todas as informações necessárias 
para aplicação dos métodos e conhecimento da frota, definindo-se, assim, os parâmetros. 
Logo após, efetuou-se aplicação dos quatro métodos. Na sequência, os resultados são 
gerados pela implementação das equações de depreciação utilizando como ferramenta de 
cálculo o Excel, através das planilhas. Por fim, se comparou os resultados e se realizou as 
discussões. 


Características da frota 


Na coleta de dados, relativos aos sessenta e três ativos do setor da construção 
civil, foi possível inferir informações que permitiram uma visão panorâmica da frota, 
convergindo para um melhor entendimento das peculiaridades dos ativos. Perspectivas 
que estabelecerão bases solidas para se ter uma previsão das possíveis magnitudes de 
valores de depreciação, no âmbito de seu grau, observadas nos gráficos das Figuras de 2 
a 5, mostrando as características da frota. 

Na Figura 2, percebe-se que quarenta e quatro ativos estão em serviço (69,84%), 
sendo a maioria dos ativos, isto é, produzindo com qualidade, ao passo que seus desgastes 
se amplificam devido seus funcionamentos contínuos naturalmente, exemplificados pelo 


fator de carga, conforme estabelece o método de Hélio de Caires o qual possui vínculo 
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com fatores de carga, pois é diretamente proporcional a depreciação. Observa-se ainda, 
que existem quinze ativos (23,81%) que estão aptos para serem vendidos e quatro ativos 
(6,35%) estão aptos em um pátio aguardando liberação de ordem de serviço. 
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Figura 2 - Status das máquinas e equipamentos na obra. 
Fonte: Própria (2016). 


Denota-se da Figura 3 uma classificação referente às categorias dos ativos em 
que os pequenos equipamentos, veículos e industrias representam minoria, observa-se 
na categoria das máquinas e equipamentos, que a maioria são rodoviários representando 
47,63%, seguido de Máquinas com 33,33%, Pequenos Equipamentos com Manutenção 
em 9,52%, Veículos com 6,35% e industriais com 3,17%. Aproximadamente, 80,96% dos 
ativos são das categorias Rodoviários ou Máquinas, cujos valores em relação as demais 
categorias são maiores, mais utilizadas e remetem uma ideia de que a depreciação total 
anual da frota possa ser alta. 
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Figura 3 - Categoria das máquinas e equipamentos. 
Fonte: Própria (2016). 


Já a Figura 4, uma classificação referente às marcas das máquinas e equipamentos, 
de igual forma da classificação das categorias, procurou-se ordem na frota em relação às 
suas marcas, em que as marcas Volkswagen (31,75%), JCB (20,63%) e Scania (12,70%) 
se destacaram. Aproximadamente, 65,08% dos ativos são das marcas já citadas, e por 
se tratarem de marcas de alto valor, com exceção da marca Volkswagen, impulsiona- 
se a um raciocínio de que depreciações de ativos são exorbitantes. As demais marcas 
apresentam percentuais equivalentes as quantidades especificadas no gráfico, ou seja, 
para cada quatro marcas tem-se 6,35%, as que apresentam duas marcas, 3,17% e as com 
uma marca apenas, 1,59%. 


MARCAS DAS MÁQUINAS E EQUIPAMENTOS 


QUANTIDADE 





Figura 4 - Marcas das máquinas e equipamentos. 
Fonte: Própria (2016). 
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Concernente a Figura 5, umas das mais importantes informações e também 
parâmetros para os métodos são os anos de fabricação, desta forma procurou-se elaborar 
uma classificação pertinente. Conforme se observa, aproximadamente 68,25% dos ativos 
foram fabricados nos anos de 2010, 2009 e 2013, respectivamente, 19,05%, 17,46% e 
15,87%, mostrando, assim, que muitos ativos estão próximos de suas vidas úteis com 
respectivamente 7, 8 e 4 anos de vida aparente, com curtas esperanças de vida. Assim, 
observando-se ainda nos anos de 2011 e 2012, a partir de seis equipamentos, os percentuais 
foram mais baixos, na faixa de 9,52%, depois aqueles com quatro equipamentos (2006, 
2007 e 2014), 6,35%, seguidos de dois equipamentos (1981 e 2004), com 3,17% e os de 
um equipamento (2008), com 1,59%. 
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Figura 5 - Ano de fabricação das máquinas e equipamentos. 
Fonte: Própria (2016). 


Tabulação de dados para aplicação dos métodos 


Depois de ter analisado as características da frota fazendo quatro pequenas 
análises, submeteu-se toda a frota aos métodos de cálculos de depreciação. Para isso 
foi elaborado uma planilha, conforme a Tabela 1, com os cabeçalhos apelido, tipo, ano de 
fabricação, sado residual, vida útil, vida aparente e valor de aquisição. A coluna apelidos foi 
utilizada, devido ao controle de identificação dos ativos ser feita por apelidos, diferente de 
outras frotas as quais são identificadas por numeração, como por exemplo em empresas 
de coleta de resíduos sólidos. 

A coluna tipo informa o tipo do ativo, se é caminhão e suas variações, 
retroescavadeira, trator e suas variações, Ônibus, veículos, moto bombas, grupo gerador e 
grupo semirreboque. À coluna ano de fabricação descreve os todos os anos de fabricação 
dos ativos. A coluna, saldo residual é utilizada para definir o valor do ativo que sofreu 
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depreciação, ao final de sua vida útil, resultando em seu valor final. A coluna vida útil remete 
a provável tempo de uso. A coluna vida aparente representa a idade do ativo. E finalmente, 
a coluna valor de aquisição mostrando o custo de ativo como novo, ou seja, o valor de 


aquisição do bem. 
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Peq-Maquina Soldar Diesel 1.100,30 7.335,35 










BATATA [Cam Auto Bomba lanca | 2012] R$ 21382071] 8] sing 109740000 
CURUJA CamBasculante | 2007] R$ 6424351] 8]  olR$ 279.221,99 
CARA PRETA 2007| R$ 6424351] *  8| 10 279.221,99 
ACACIA [Cam Basculante | 2010/ R$ 5157157] Bs) q 206.075,50 
CRONOS [Cam Basculante | 2013) R$ 106277223] | a, 500.663,30 
JURUPOCA [Cam Basculante | 2013/R$ 10627723] | a] 489.000,00 
R10 [Cam Basculante | 2013] R$ 10627723] Bs] a 489.787,33 
TARDELU! [Cam Basculante | 2013/R$ 10627723] 8) a] 489.000,00 
MARLONE CamBasculante | 2014/R$ 10627723) *  &| 505.000,00 
EGIDIO [Cam Basculante | 2014] R$ 106277223] 8) a 505.000,00 
FABIAO Cam Basculante | 2014] R$ 10627723] Bs] 3 505.000,00 
HENRIQUES11DR |CamBasculante | 2014/ R$ 10627723] Bs] a 505.000,00 
BUCHUDO CamBau |] 2004) R$ 2927196] —  B| as 99.539,76 
VIVALDAO [Cam Betoneira | 2010/ R$ 6161507] Bs] q 291.153,00 
VIVALDINHO [CamBetoneira | 2o1]R$ 6161507) 8a) 7 291.153,00 
ARENA CamBetoneira | 2010] R$ 6161507] 8) [R$ 29115300 
MARINGA [Cam Carroceria | 2006/R$ 3100196] 8) ai 192.450,31 
MALARIA IGAM3  |CamCarroceria | 2006 R$ 3100196] | as] 192.450,31 
ARENAPOLIS [CamGuindauto | 2010/ R$ 4720136] Bs] q 231.559,27 
BREU PRETO [Cam Guindauto | 2007/ R$ 5772466] | o 280.755,07 
ABARE [CamGuindauto | 2013] R$ 7800550] | a 393.300,00 
ITAPECIRICA [Cam Guindauto | 2010] R$ 4720136] 8) q 231.559,27 
CORURIPE CamGuindauto | 2010)R$ 4720136) 8) 7 231.559,27 
AZEITE 2011] R$ 5659995] 8) 6 280.609,23 
AGUAPE Cam TanquePipa | 2012] R$ 4100093] | Ss 202.427,00 
GUARAVERMELHO [CamTanquePipa | 2010) R$ 5814867] 8) q 301.523,11 
ENCHENTE [CamTanquePipa | 2013) R$ 5814867] | a 276.000,00 
MANAUARA Cavalo Mecanico | 2010/ R$ 4046085] | q 197.558,52 
MANAOS Central ConcretoFixa | 2011] R$ 3263184] a) 6 434.821,00 
SUZANA Escavadeira | 2012]R$ s8sng9) a) Ss 426.881,88 
BELEM Ill Escavadeira | 2ois[R$ o289754] | a 489.000,00 
SANTOS Grua | 2or[R$ 45558776] 10] R$ 322251008 
ILUMINADO GrupoGerador | 201/ R$ 1032193] Bs] é 76.393,57 
INDAIA Motoniveladora | 2013) R$ 14657135] Bs] a 706.100,00 
MAURILIO Onibus |] 2o1n[R$ 3316007] | 244.798,14 
RAICA [R$ 174000] 8) a 11.600,00 
VIRA VIRA [R$ 174000] 8) o 11.600,00 

[Peq-Maquina Soldar Diesel | R$ 110030] 8] as 

; pr 

Peg-Motobomba |] [RS 353,34 | Ds 

Peg-Motobomba |] es 35334) | 36 

Retro Escavadeira |] R$ aris44s| 8] 

Retro Escavadeira |] oa 

Retro Escavadeira |] oa 

Retro Escavadeira |] E 

Retro Escavadeira |] E 

Retro Escavadeira TR 

Retro Escavadeira E 

Retro Escavadeira E 

Retro Escavadeira |] oo Ca 

Retro Escavadeira | [os] 

Rolo Compactador |] E 

Rolo Compactador |] E 

[os a 

Rolo Compactador |] Do 

Rolo Compactador |] TR: 

Semi Reboque [os 

rrator Esteiras] E E 

rrator Esteiras E 

Veiculo 









2BBB|B|B/6|B|8|6|8|8|6|8|8|6|8|8|6|8|8|2|8|8 | |S|B Ba a 8/88 jaja la jajala ja] a |a|a a |a|a 


| 2006 
NORMA eq-Motobomba | 1981 2.355,60 
AURENY Peq-Motobomba 2.355,60 
IGARAPE Peq-Motobomba 2.355,60 
BOCA DO ACRE Retro Escavadeira | 2009] R$ 41.184,49 169.950,00 
BARREIRINHA Retro Escavadeira | 2009] R$ 41.184,49 169.950,00 
PAUINI Retro Escavadeira | 2009] R$ 47.480,86 | 210.291,60 
NOVA OLINDA Retro Escavadeira | 2009] R$ 41.184,49 | 169.950,00 
BERURI Retro Escavadeira | 2009] R$ 4118449] 169.950,00 
ARINI Retro Escavadeira | 2009] R$ 4118449] 169.950,00 
CANACARI Retro Escavadeira | 2009] R$ 4118449] 169.950,00 
BEL Retro Escavadeira | 2009] R$ 4118449] 169.950,00 
BORBA Retro Escavadeira | 2009] R$ 41.184,49 | 169.950,00 
BALEIA Retro Escavadeira | 2009] R$ 41.184,49 | 169.950,00 
ZOVIS Rolo Compactador | 2009] R$ 5109354| 302.000,00 
INHOMIRIM Rolo Compactador | 2010] R$ 4308114 279.064,80 
COTUNDUBA Rolo Compactador 279.064,80 
SOCADOR Rolo Compactador 121.895,86 
PILAO S11DR Rolo Compactador | 2013] R$  18.587,04| 121.895,86 
ATUMAM Semi Reboque 2011] R$ 28.290,00 125.231,86 
PACU Trator Esteiras | 2006] R$ 105.283,65 | 605.000,00 
DONI Trator Esteiras 828.673,00 
MATRINXA | 2010/ R$ 1300500) 3) q 90.818,80 
GUASCA Veiculo | 20:2]R$ sa9876) 3) OS 49.695,15 
CARITE Veiculo | 20:2]R$ cao) 3) SS 44.434,89 
REDOMAO Veiculo | 2o]rs cs7700) a) OS 44.330,18 
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Tabela 1- Parâmetros a serem utilizados nos métodos. 
Fonte: Própria (2016 
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31 RESULTADOS E DISCUSSÃO 


Após utilização das metodologias obtêm-se os resultados das depreciações das 
máquinas e equipamentos, conforme Tabela 2. Neste momento é possível comparar as 
depreciações de cada ativo obtido através dos quatro métodos utilizados. A Tabela 2 conta 
com 4 colunas sendo a primeira exibindo os resultados das depreciações pelo método 
linear; a segunda pelo método de Cole ou soma de dígitos; a terceira pelo método da 
percentagem constante: e a última pelo método de Hélio de Caires. É notório que os 
valores, por cada método, são diferentes, uma vez que os gráficos representativos dos 
métodos são diferentes. Pelo fato do método linear depreciar o mesmo valor, seu gráfico 
resulta de uma linha continua e decrescente, quanto aos demais são gerados gráficos 
exponenciais, culminando num decaimento superior ao método de linha e diferentes entre 


si. Assim, destas diferenças gráficas implicam valores discrepantes. 
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MÉTODOS DE DEPRECIAÇÃO 


Método Linear Método P. Constante Método Hélio de Caires 


RS 552.237,06 R$ 11.112,70 
RS 6424351 R$ 64.243,51 
RS 6424351 R$ 64.243,51 
R$ 13519111 R$ 0,00 
RS 197.193,04 R$ 50.820,74 
RS 19136139 R$ 49.636,84 
RS 19175505 R$ 49.716,76 
RS 19136139 R$ 49.636,84 
RS 14952104 R$ 172.245,67 
RS 149.521,04 R$ 172.245,67 
R$ 14952104 R$ 172.245,67 
R$ 149.521,04 R$ 172.245,67 
RS 2927196 R$ 29.271,96 
RS 200.845,69 R$ 0,00 
RS 200.845,69 R$ 0,00 
R$ 200.845,69 R$ 0,00 
RS 3100196 R$ 31.001,96 
RS 3100196 R$ 31.001,96 
R$  161313,17 R$ 0,00 
RS 5772466 R$ 57.724,66 
RS 157.647,25 R$ 39.922,63 
RS  161313,17 R$ 0,00 
RS 16131317 R$ 0,00 
R$ 168.006,96 R$ 24,73 
R$ 100.891,29 R$ 2.049,85 
RS 21295264 R$ 0,00 
RS 108.925,67 R$ 28.015,88 
RS 137.460,46 R$ 0,00 
RS 301.641,87 R$ 38,32 
RS 21125631 R$ 4.322,77 
RS 19805133 R$ 49.636,84 
R$ 1.660.153,39 R$ 58.742,49 
RS 4955373 R$ 6,73 
RS 27976433 R$ 71.673,97 
RS 15872855 R$ 21,58 
RS 1.740,00 R$ 1.740,00 
R$ 1.740,00 R$ 1.740,00 
RS 1.100,30 R$ 1.100,30 
R$ 353,34 [RS 35334 /RS 35334 | R$ 353,34 
R$ 35334 [RS 35334 /RS 35334 /R$ 353,34 
R$ 353,34 [RS 35334 [RS 35334 [R$ 353,34 
RS 12876551 R$ 41.184,49 
RS 12876551 R$ 41.184,49 
RS 162.810,74 R$ 47.480,86 
RS 12876551 R$ 41.184,49 
RS 128.765,51 R$ 41.184,49 
RS 12876551 R$ 41.184,49 
RS 12876551 R$ 41.184,49 
RS 12876551 R$ 41.184,49 
RS 12876551 R$ 41.184,49 
RS 128.765,51 R$ 41.184,49 
RS 250.906,46 R$ 51.093,54 
RS 206.485,70 R$ 0,00 
RS 206.485,70 R$ 0,00 
RS 5165441 R$ 12.373,26 
RS 5165441 R$ 12.373,26 
RS | 7270640 R$ 11,04 
RS 105.283,65 R$ 105.283,65 
RS 130.687,50 R$ 130.687,50 
R$ | 13.005,00 R$ 13.005,00 
RS 8.198,76 R$ 8.198,76 
R$ 6.477,00 R$ 6.477,00 
R$ 6.477,00 R$ 6.477,00 


Tabela 2 - Resultados das depreciações por cada método. 
Fonte: Própria (2016). 


Seguindo a mesma ideia de apresentação dos resultados da Tabela 2, tem-se a 
Tabela 3 com a apresentação dos valores de mercado obtidos pelas metodologias propostas. 
Os valores de mercado de cada máquina e equipamento são resultantes da diferença entre 


os seus valores de aquisição e suas depreciações. Alguns valores resultaram em zero 
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pelo método e Hélio de Caires, isto se deve a característica matemática da equação, pois 


quando a depreciação tende a zero, à medida que vai se aproximando da vida útil, faz com 


que o denominador da equação (4) fique muito alto (em torno da ordem de dez a décima 


quinta potência), em virtude da equação (3) (equação de desgaste) também ficar elevada, 


por conta do expoente do número de Euler deixa-lo com valor alto. 


VALORES DE MERC ADO 


Método Linear Método P. Constante | Método Hélio de Caires 


R$ 
R$ 


R$ 
R$ 


õ [6 [ô [6 |d [6 |d [6 [6 |6 [6 |6 |6 |6 |6 | |6 | 6 |8 |6 |8 6 |8 8 |8 6 |8 8 [68 |6]8 [6] |6]8 |6]8 |6]8 |6]8 |6]8|6]8|6]8|6]8|6]8|6]8 


R$ 
R$ 
R$ 
R$ 


Tabela 3 - Resultados dos valores de mercado para cada método. 


545.162,94 R$ 
214.978,48 R$ 
214.978,48 R$ 
70.884,39 R$ 
303.470,27 R$ 
297.638,62 R$ 
298.032,28 R$ 
297.638,62 R$ 
355.478,96 R$ 
355.478,96 R$ 
355.478,96 R$ 
355.478,96 | R$ 272.411,72 R$ 
70.267,80 R$ 
90.307,31 R$ 
90.307,31 R$ 
90.307,31 R$ 
161.448,35 R$ 
161.448,35 R$ 
70.246,10 R$ 
223.030,41 R$ 
235.652,75 R$ 
70.246,10 R$ 
70.246,10 R$ 
112.602,27 R$ 
101.535,71 R$ 
88.570,48 |R$ 64.909,07 R$ 
167.074,34 R$ 
60.098,06 R$ 
133.179,13 R$ 
215.625,57 R$ 
290.948,67 R$ 


1.562.356,69 | RS 958.664,55 | R$ 2.226.139,39 | R$ 


26.839,84 R$ 
426.335,68 R$ 
86.069,59 R$ 

9.860,00 R$ 

9.860,00 R$ 

6.235,05 |R$ 623505 /RS 623505 |R$ 

2.002,26 R$ 

2.002,26 R$ 

2.002,26 R$ 
41.184,49 R$ 
41.184,49 R$ 
47.480,86 [RS 47.480,86 R$ 
41.184,49 R$ 
41.184,49 R$ 
41.184,49 R$ 
41.184,49 R$ 
41.184,49 R$ 
41.184,49 R$ 
41.184,49 R$ 
51.093,54 R$ 
72.579,10 R$ 
72.579,10 R$ 
70.241,45 R$ 
70.241,45 R$ 
52.525,47 R$ 
499.716,35 R$ 
697.985,50 R$ 
77.813,80 R$ 
41.496,39 R$ 
37.957,89 R$ 
37.853,18 R$ 


Fonte: Própria (2016). 
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1.086.287,30 
214.978,48 
214.978,48 
206.075,50 
449.842,56 
439.363,16 
440.070,57 
439.363,16 
332.754,33 
332.754,33 
332.754,33 
332.754,33 
70.267,80 
291.153,00 
291.153,00 
291.153,00 
161.448,35 
161.448,35 
231.559,27 
223.030,41 
353.377,37 
231.559,27 
231.559,27 
280.584,50 
200.377,15 
301.523,11 
247.984,12 
197.558,52 
434.782,68 
422.559,11 
439.363,16 


3.163.767,59 


76.386,84 
634.426,03 
244 776,56 

9.860,00 

9.860,00 

6.235,05 

2.002,26 

2.002,26 

2.002,26 
128.765,51 
128.765,51 
162.810,74 
128.765,51 
128.765,51 
128.765,51 
128.765,51 
128.765,51 
128.765,51 
128.765,51 
250.906,46 
279.064,80 
279.064,80 
109.522,60 
109.522,60 
125.220,82 
499.716,35 
697.985,50 

77.813,80 

41.496,39 

37.957,89 

37.853,18 
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Com o auxílio da Figura 6, fica claro que em quase todas as máquinas e equipamentos 
existe diferença ou distorção entre as metodologias de depreciação propostas, ou seja, a 
utilização da metodologia vai depender da finalidade da avalição o método a se escolher. 

Conforme se observa, o primeiro lugar com R$11.929.081,66 de depreciação da 
frota, é produto das variações de grandes depreciações resultando numa depreciação 
total maior em relação as demais. Em segundo lugar, tem-se a depreciação do método 
linha com R$9.597.903,41, pois não há variações de depreciações anuais, são quotas 
constantes, resultando numa depreciação menor do que a do método de cole. Em terceiro 
lugar, tem se a depreciação pelo método da percentagem constante com R$6.563.011,91, 
muito embora o valor de aquisição seja multiplicado por uma taxa constante, aquele 
sempre se atualiza no ano subsequente, resultando em variações depreciação diferentes 
no decorrer do tempo, porém, ainda assim, menor que as variações referentes ao método 
de cole. Em quarto e último lugar, tem-se a depreciação pelo método de Hélio de Caires 
com R$2.195.971,18, representando o menor valor entre os métodos, deve-se ao fato da 
elaboração da equação de Caires ser baseada numa amostra de 850 ativos móveis os 
quais podem diferir das características da frota escolhida para esta pesquisa, portanto 
retornado valores distorcidos, de forma a serem menores. 


Depreciação total da frota por método 


R$ 15.000.000,00 - 





R$ 11.929.081,66 


R$ 9.597.903,41 





g R$ 10.000.000,00 
a 


R$ 6.563.011,91 





R$ 5.000.000,00 
R$ 2.195.971,18 





Método de Cole Método Linear Método P. Método Hélio de 
Constante Caires 


MÉTODOS 








Figura 6 - Depreciação total da frota por métodos. 
Fonte: Própria (2016). 


Através da Figura 7 também é possível visualizar que dentre as máquinas e 
equipamento da amostra que existem alguns ativos cuja vida aparente (idade atual) já 
ultrapassara a vida útil, fazendo com que o ativo atinja o valor residual (valor de sucata), 
não podendo desta maneira, reduzir o seu valor abaixo desta referência. Ainda na Figura 
6, ilustra as quantidades de ativos sucateados e os que ainda não atingiram o tal valor, do 
total da amostra de 63 bem móveis, 18 já se encontram neste estado. 
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Situação das máquinas e 
equipamentos 


QUANTIDADE 


ATIVO SUCATA 
SITUAÇÃO DO EQUIPAMENTO 





Figura 7 - Situação das máquinas e equipamentos. 
Fonte: Própria (2016). 


Para finalizar o estudo é apresentada a Figura 8, de modo a identificar as 
próximas máquinas e equipamentos que atingirão o valor de residual no ano de 2018. 
Do total da amostra, onze ativos estarão neste estado de valor, sendo que dos 11, 10 são 


retroescavadeiras e um é rolo compactador. 


Projeção de Equipamentos 
Sucateados para o ano de 2018 


1 
q 
m 
= 
mi! 
— 
E 
É 


Retroescavadeira Rolo compactador 


Tipos de Equipamentos 





Figura 8 - Projeção de Equipamentos sucateados para o ano de 2018. 
Fonte: Própria (2016). 


41 CONCLUSÃO 


O que mais chamou atenção no desenvolvimento deste artigo, foram os cuidados 
em relação a depreciação a serem tomados, concernente aos métodos de depreciação 
a serem aplicados em cada um dos bens móveis. Cada citação do referencial teórico foi 
essencial para que se pudesse aprofundar no assunto e encontrar ainda mais informações. 
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Trabalhar com avaliação exige do engenheiro avaliador grande experiência acumulada, 
pois até o presente momento não se tem um método de uso geral, sendo composto de 
muito empirismo, como também pesquisas e bom senso. 

Em se tratando de máquinas e equipamentos, a avaliação se mostra complexa 
inerente ao grande leque de tipos ativos, entretanto pode ser sanado com os métodos 
aplicados. Por isso o engenheiro avaliador deve proceder com a coleta de dados de 
maneira minuciosa, verificando os valores reais de mercado, retornado as características 
dos bens móveis como práticas de manutenção, manuais, entre outros os quais servirão de 
parâmetros para o método apropriado. 

Durante a metodologia, diferentes métodos foram utilizados e conduziram 
a diferentes valores de depreciação remetendo-se a diferentes valores de mercado, 
oferecendo, desta forma, um auxílio para engenheiros mecânicos que lidam com estas 
avaliações de bens móveis. 

Através dos resultados obtidos ficou claro que os valores de depreciações pelos 
métodos de cole, percentagem constante e linear foram superiores ao de Caires. Desta 
forma, pode-se concluir que calcular depreciação por uma gama de métodos constitui uma 
ótima alternativa para uma aproximação ao valor de mercado do bem a ser avaliado. 

Outra grande contribuição da metodologia proposta, além de calcular os valores de 
mercado e a depreciação total da frota, foi de identificar a quantidade de equipamentos 
que atingiram o valor residual e os próximos irão atingir no ano subsequente (2018). De 
modo a melhorar ainda mais a qualidade de resultados oriundos da avaliação os quais 
permitirão, caso necessário, tomada de decisão em relação a futuro do ativo, se vende, 


aluga, descarta ou continua-se com bem em frente serviço. 
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RESUMO: Este artigo tem como objetivo 
analisar os impactos gerados pela falta de 
compatibilização de projetos no custo de 
uma creche tipo 1 padrão FNDE em etapa de 
execução localizada em Caruaru-PE. Para tanto, 
foi realizada a identificação e quantificação das 
incompatibilidades entre os projetos (básico e 
complementares) observadas durante a fase 
de execução de uma creche tipo 1 padrão 
FNDE. Com o intuito de levantar o impacto 
causado no custo da obra. Inicialmente, foram 
identificadas 72 incompatibilidades in loco 
durante a execução, sendo 74% entre os projetos 
estrutural e hidrossanitário, 22% entre os projetos 
arquitetônico e hidrossanitário e 4% entre os 
projetos arquitetônico e elétrico. Com base nas 
interferências levantadas, foi elaborada uma 
planilha de custo para os serviços que precisaram 
ser acrescidos, para as soluções adotadas 
durante a execução. Ao final, foi analisado o 
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FNDE EM CARUARU-PE 


custo acrescido devido as incompatibilidades 
entre os projetos, totalizando um custo de R$ 
7.591,08, (sete mil, quinhentos e noventa e 
um reais e oito centavos). Na discussão dos 
resultados indicou-se a utilização da ferramenta 
BIM para evitar problemas de incompatibilidades 
entre os projetos, visto que o mesmo é um projeto 
padrão replicado nacionalmente, ou seja, o custo 
acrescido é bem mais expressivo. 
PALAVRAS-CHAVE: BIM. Incompatibilidade. 
Custo. Obra pública. 


ANALYSIS OF THE IMPACTS 
GENERATED BY THE LACK OF 
COMPATIBILIZATION OF PROJECTS AT 
THE COST OF THE EXECUTION OF A 
TYPE 1 STANDARD FNDE NURSERY IN 
CARUARU-PE 


ABSTRACT: This article aims to analyze 
the impacis generated by the lack of project 
compatibility on the cost of a FNDE standard type 
1 daycare center in the execution stage located 
in Caruaru-PE. To this end, an identification and 
quantification of incompatibilities between the 
projects (basic and complementary) observed 
during a phase of execution of a type 1 daycare 
standard FNDE was carried out. In order to raise 
the impact impact on the cost of the work. Initially, 
72 incompatibilities were identified on the spot 
during execution, 74% between structural and 
hydrosanitary projects, 22% between architectural 
and hydrosanitary projecis and 4% between 
architectural and electrical projects. Based on 
the interferences raised, a cost spreadsheet 
was prepared for the services that needed to 
be added, for the solutions adopted during an 
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execution. In the end, the increased cost due to incompatibilities between the projects was 
analyzed, totaling a cost of R $ 7,591.08, (seven thousand, five hundred and ninety-one reais 
and eight cents). In the discussion of the results, the use of the BIM tool was indicated to avoid 
problems of incompatibilities between the projecis, since it is a standard project replicated 
nationally, that is, the added cost is much more expressive. 

KEYWORDS: BIM. Incompatibility. Cost. Public work. 


11 INTRODUÇÃO 


No Brasil, a construção civil enfrenta dificuldades derivadas do processo empírico 
de construções. A utilização de metodologias antigas tanto no gerenciamento quanto na 
execução pode resultar em atrasos e prejuízos para os empreendimentos. 

O gerenciamento dos custos é um conjunto de processos envolvidos em 
planejamento, estimativas, orçamentos, financiamentos, gerenciamento e controle dos 
custos, de modo que o projeto possa ser terminado dentro do orçamento aprovado. Pode ser 
utilizado algumas técnicas para prevenção de erros orçamentários como: análise de valor 
agregado, índice de estimativa de custos da obra, índice de estimativa de desempenho, 
previsão de custos futuros (em caso de atrasos) e tomadas de decisão (PMI, 2017). 

Segundo o TCU (2014), obra pública é toda construção do bem público, que pode 
ser executada pelo próprio órgão (forma direta) ou licitada para um terceiro construir (forma 
indireta). 

De acordo com Paiva (2016), compatibilizar é incorporar tudo que envolvem um 
projeto. Após a conclusão de todas as partes envolvidas inicia-se o processo, sanando 
conflitos e interferências através de sobreposição manual ou por meio de tecnologias 
capazes de acelerar a verificação. 

A incompatibilidade entre projetos é um dos fatores que mais gera desperdícios em 
uma obra. Quando solucionadas na fase de elaboração evita-se problemas no canteiro de 
obras, reduzindo de 5% a 8% os custos da construção. Com a implantação da ferramenta 
Building Information Modeling (BIM) é possível prever erros com mais facilidade (INBEC, 
2018). 

Segundo o Ministério da Educação - MEC, na resolução nº 6, de 24 de abril de 
2007, o Programa Nacional de Reestruturação e Aquisição de Equipamentos para a Rede 
Escolar Pública de Educação Infantil (Proinfância), é uma das medidas tomadas pelo 
Governo, visando garantir o acesso de crianças a creches e escolas, bem como a melhoria 
da infraestrutura física da rede de Educação Infantil (FNDE, 2017). 

O objetivo geral foi analisar os impactos gerados pela falta de compatibilização de 
projetos no custo de uma creche tipo 1 padrão FNDE em etapa de execução localizada em 
Caruaru-PE. 

Foram identificados e quantificados os erros de compatibilização através de dados 
coletados in loco, durante a etapa de execução. Onde foi a primeira creche desse tipo 
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empregada no município, que atualmente contempla outras quatro unidades, três concluídas 
e duas em fase de construção. 

Todas foram projetadas antes do decreto, que estabelece o uso da metodologia 
BIM em obras públicas, ou seja, foram projetadas sem a metodologia BIM, resultando em 
inconformidades entre os projetos e impactando no custo do empreendimento. 


21 METODOLOGIA 


Adotou-se uma metodologia de caráter exploratório baseado em um estudo de 
caso com dados de natureza quantitativa e qualitativa. Para esta metodologia foi realizada 
uma revisão bibliográfica para investigar, compreender, analisar e discutir os benefícios 
do sistema BIM se implementado em obras públicas no Brasil. Com o intuito de identificar 
e quantificar os problemas ocasionados pela falta de compatibilização dos projetos 
arquitetônico, estrutural, elétrico e hidrossanitário durante o processo de execução e 
analisar os impactos no custo de uma creche tipo 1 padrão FNDE situada no bairro de 
Nova Caruaru em Caruaru-PE. 

A pesquisa exploratória é quando a pesquisa se encontra na fase preliminar, ou 
seja, investiga e adiciona informações sobre a explicação do assunto afim de apresentar 
as evidências, além de descobrir novos fenômenos e formular novas ideias, através de 
pesquisa bibliográfica e estudo de caso. O estudo de caso é um tipo de pesquisa que busca 
identificar, analisar e explicar como e por que ocorre um fenômeno em particular, sendo 
assim o pesquisador aprofunda o estudo (PRODANOV; FREITAS, 2013). 

A pesquisa quantitativa trata-se de uma abordagem que emprega padrões 
sistemáticos com respostas pré-estabelecidas para facilitar a análise e comparação dos 
dados e obter resultados quantificáveis (BASTOS; KELLER, 2015). Por meio de dados 
quantificáveis foram reunidos os problemas gerados pela falta de compatibilização entre os 
projetos arquitetônico, estrutural, elétrico e hidrossanitário durante o processo de execução 
de uma creche tipo 1 padrão FNDE, que interferiram diretamente no custo. 

O embasamento teórico possibilitou amadurecer o conhecimento sobre o tema 
da pesquisa e forneceu o suporte necessário para o desenvolvimento de uma avaliação, 
da relação entre os problemas gerados pela incompatibilização dos projetos durante o 
processo de execução e os impactos no custo, através de registros fotográficos e planilha 
de itens acrescidos. 

A coleta de dados referentes ao objeto de estudo foi realizada com base em uma 
avaliação in loco em uma creche tipo 1 padrão FNDE situada no Bairro de Nova Caruaru 
em Caruaru-PE. O município de Caruaru-PE está a aproximadamente 119 km de distância 
da capital do estado de Pernambuco (Recife). 

Durante o processo de execução, foi realizado o levantamento a fim de identificar, 
registrar e quantificar as inconformidades derivadas das incompatibilizações entre o projeto 
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básico e os projetos complementares elaborados através de uma ferramenta 2D. Para 
analisar os impactos causados no custo, foram levantados os itens que precisaram ser 
acrescidos por causa das incompatibilidades e tomados como parâmetro de valores das 
composições apresentadas na planilha da própria creche. 

O município Caruaru-PE contempla 5 creches tipo 1 padrão FNDE em execução, 
nas quais estão apresentadas a localização na Figura 1, a creche abordada foi a creche 
sinalizada de vermelho. O terreno possui 40x60m, totalizando 2.400,00 m? e a área 
construída é de 1.311,97 m?2. Os projetos elaborados pelo FNDE para execução desta 
creche foram de 2015 e sua execução só foi iniciada no ano de 2017, sendo a primeira 


creche tipo 1 padrão FNDE do município. 


MAPA DAS CRECHES TIPO 1 PADRÃO FNDE Em gpa 
CARUARU - PE ER [E] CRECHE TPO 1 PADRÃO FNDE 
na à fo, e E CRECHE TPO 1 PADRÃO FNDE ESTUDADA, 


(ET) AE TA 
[uva i Fe A pur | a nan 





Figura 1 — Mapa das creches tipo 1 padrão FNDE em Caruaru-PE 
Fonte: Google Earth (2020). 


A Figura 2 apresenta um fluxograma da metodologia científica adotada para 
o trabalho, exemplificando as duas áreas de abordagem do trabalho, qualitativa e 
quantitativa. 

Para apresentar os resultados foram analisados e discutidos dados coletados em 
conjunto com o referencial teórico. Em seguida, planilhas sintéticas foram produzidas 
no Excel, composta pelos serviços acrescidos para solucionar as incompatibilidades e 
pelos valores unitários dos serviços com base na própria planilha da creche, onde seus 
valores são baseados em: SINAPI/PE 06/2017 desonerada, ORSE 06/2017, EMLURB 
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12/2014 e SEINFRA. O BDI adotado foi de 27,7%. E foram produzidos quadros com os 
erros identificados durante a execução nos projetos arquitetônico, estrutural, elétrico e 


hidrossanitário de uma creche tipo 1 padrão FNDE. 


VISITA IN LOCO 


ANÁLISE DOS IMPACTOS 
CAUSADOS PELA FALTA 
DE COMPATIBILIZAÇÃO ANÁLISE QUALITATIVA 
ENTRE OS PROJETOS 


IMPACTOS NO CUSTO QUANTIFICAÇÃO DAS 
INCOMPATIBILIDADE 
ANÁLISE DOS DADOS 











ANÁLISE QUANTITATIVA 


ERROS DE 


COMPATIBILIZAÇAO 


r 


ANALISE DOS DADOS 





SÍNTESE DAS ANÁLISES 


Figura 2 — Mapa das creches tipo 1 padrão FNDE em Caruaru - PE 
Fonte: Autor (2020). 


31 RESULTADOS E DISCUSSÕES 


Para obtenção dos resultados foi feito um levantamento e análise das interferências 
que ocorreram por falta de compatibilização entre os projetos, durante a execução da 
creche objeto de estudo. Após isso, foram analisados todos os problemas que ocorreram, 
as soluções adotadas e levantado o custo que levou para solucionar os problemas. 


3.1 Levantamento de interferências de projetos 


O processo de identificação das incompatibilidades entre os projetos foi realizado 
in loco durante o processo de execução. Após isso, foi levantada a quantidade de 
incompatibilidades identificadas entre os projetos arquitetônico, estrutural, elétrico e 
hidrossanitário. 

Em seguida, foram analisados os problemas individualmente e as respectivas 
soluções adotadas durante a etapa de execução, com a finalidade de levantar os impactos 
causados no custo pela falta de compatibilização. Os dados obtidos foram agrupados e 


sintetizados como apresenta o Tabela 1. 


Projeto, Análise e Otimização na Área das Engenharias Capítulo 6 EE 


DESCRIÇÃO QUANTIDADE PERCENTUAL (%) 


ARQUITETÔNICO / ELÉTRICO 3 4 
ARQUITETÔNICO / HIDROSSANITÁRIO 16 22 
ESTRUTURAL / HIDROSSANITÁRIO 53 74 
TOTAL 72 100 


Tabela 1 - Incompatibilidade entre projetos 


Fonte: Autor (2020). 


As interferências encontradas entre os projetos estrutural e hidrossanitário, 
subdividem-se em estrutural e água fria, estrutural e esgoto sanitário. Os problemas 
encontrados entre os projetos arquitetônico e hidrossanitário, resumem-se em arquitetônico 
e água pluvial. As incompatibilidades encontradas entre os projetos arquitetônico e elétrico, 


resumem-se em esquadria e elétrico, de acordo com a Tabela 2. 


DESCRIÇÃO QUANTIDADE PERCENTUAL (%) 
ESTRUTURAL / ÁGUA FRIA 14 19,44 
ESTRUTURAL / ESGOTO SANITÁRIO 39 54,16 
ARQUITETÔNICO / ÁGUA PLUVIAL 16 RRjop 
ESQUADRIA / ELÉTRICO 3 4,18 
TOTAL 72 100 


Tabela 2 - Relação de interferências 
Fonte: Autor (2020). 


Detalhadamente e respectivamente, os problemas identificados durante a execução 
da creche foram: os ramais de abastecimento de água fria coincidiam com as vigas 
superiores, as colunas de ventilação da rede sanitária coincidiam com as vigas superiores, 
os ramais de água pluvial não coincidiam com o projeto arquitetônico, os pontos de tomada 


coincidiam com as esquadrias, conforme apresenta o Figura 3. 
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* O ramal de abastecimento coincidia com a viga 
superior 


* À coluna de ventilação coincidia com a viga 
superior 


- O ramal de água pluvial não coincidia com o 
projeto arquitetônico 


* O ponto de tomada coincidia com o vão de janela 


* O ponto de tomada coincidia com o vão de porta 





Figura 3 - Problemas encontrados 
Fonte: Autor (2020). 


Para os problemas apresentados entre as vigas superiores e os ramais de 
abastecimento de água fria a solução adotada foi desviar os ramais de abastecimento de 
àgua fria da viga superior, deixando as tubulações aparentes (Figura 4). Para os problemas 
apresentados entre as vigas superiores e as colunas de ventilação da rede sanitária as 
soluções adotadas foram relocadas as colunas de ventilação da rede sanitária para shaft's 
e paredes, criar shaft's para relocar as colunas de ventilação da rede sanitária e desviar as 
colunas de ventilação da rede sanitária da viga superior, deixando as tubulações aparentes 
(Figura 5). 





” da 


Figura 4 e 5 - Incompatibilidade entre estrutural e água fria e Incompatibilidade entre estrutural 
e esgoto sanitário. 


Fonte: Autor (2020). 
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Para os problemas apresentados entre o projeto arquitetônico e os ramais de água 
pluvial a solução adotada foi criar shaft's para esconder os ramais de água pluvial (Figura 
6). Para os problemas apresentados entre os vãos de esquadria e os pontos de tomada a 
solução adotada foi relocar o ponto de tomada para lateral e para baixo (Figura 7). 





Figura 6 e 7- Incompatibilidade entre arquitetônico e água pluvial e Incompatibilidade entre 
esquadria e elétrico 


Fonte: Autor (2020). 


Observou-se as soluções adotadas para os problemas identificados durante a 
execução da creche (Figura 8). 
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o ponto de tomada foi relocado para baixo do vão de janela % 


o ponto de tomada foi relocado para lateral esquerda do vão de 
janela 


/o 


o ponto de tomada foi relocado para lateral esquerda do vão de 
porta 


Vo 


foi criado um shaft para esconder a decida de água pluvial 


a coluna de ventilação foi desviada da viga superior deixando a AG 
tubulação aparente ” 


foi criado um shaft para relocar a coluna de ventilação 
a coluna de ventilação foi relocada para uma parede 


a coluna de ventilação foi relocada para o shatf 


o ramal de abastecimento foi desviado da viga superior deixando 
a tubulação aparente 


0.00% 5.00% 10.00% 15.00% 20.00% 25.00% 


Figura 8 - Soluções adotadas 
Fonte: Autor (2020). 


A solução adotada para os problemas entre as vigas superiores e os ramais de 
abastecimento de água fria foi desviar os ramais de abastecimento de água fria de 950mm, 
()60mm e 975mm da viga superior, deixando as tubulações aparentes. Para os problemas 
entre as vigas superiores e as colunas de ventilação da rede sanitária foi relocar as colunas 
de ventilação da rede sanitária de 450mm e 975mm para shaft's, criar shaft's para as 
colunas de ventilação e desviar a coluna de ventilação da rede sanitária da viga superior. 

Os problemas entre o projeto arquitetônico e os ramais de água pluvial foi 
solucionado criando shaft's para esconder os ramais de água pluvial. A solução adotada 
para os problemas entre os vãos de esquadria e os pontos de tomada foi relocar o ponto de 
tomada para lateral esquerda e para baixo do vão de janela, respectivamente. A solução 
adotada para os problemas entre os vãos de esquadria e os pontos de tomada foi relocar o 
ponto de tomada para lateral esquerda do vão de porta. 


3.2 Avaliação das incompatibilidades no custo 


Após o levantamento das incompatibilidades entre os projetos observadas durante 
a execução da creche, foram verificadas as soluções adotadas para os problemas, 
para criar uma planilha sintética composta pelos serviços acrescidos para solucionar as 
incompatibilidades e pelos valores unitários dos serviços com base na própria planilha da 
creche. 

Os serviços acrescidos para solucionar os problemas totalizam R$ 7.591,08, (sete 
mil, quinhentos e noventa e um reais e oito centavos), considerando que se trata de um 
projeto padrão replicado nacionalmente, entende-se que o custo acrescido é bem mais 
expressivo. A Figura 9 apresenta o percentual de custo de cada serviço acrescido. 
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O serviço que exigiu maior aplicação de recursos foi o serviço de sistema de vedação 
vertical interno e externo (paredes) totalizou-se R$ 2.498,61 (dois mil, quatrocentos e 
noventa e oito reais e sessenta e um centavos) referente à 33%. No serviço de drenagem 
de águas pluviais totalizou-se R$ 1.266,69 (mil, duzentos e sessenta e seis reais e 
sessenta e nove centavos) referente à 17%. No serviço de instalação sanitária totalizou-se 
R$ 1.128,66 (mil, cento e vinte e oito reais e sessenta e seis centavos) referente à 15%. 
No serviço de instalação hidráulica totalizou-se R$ 1.093,60 (mil e noventa e três reais e 
sessenta centavos) referente à 14%. 


29% * sistema de vedação vertical interno e extemo 
E (paredes) 
* revestimentos internos e externos 
pintura 
- instalação hidráulica 
À - drenagem de águas pluviais 
0 
5% - instalação sanitária 
q - instalações elétricas - 220v 


Figura 9 - Percentual de custo por serviço 
Fonte: Autor (2020). 





Nos serviços de revestimentos internos e externos totalizou-se R$ 1.075,87 (mil e 
setenta e cinco reais e oitenta e sete centavos) referente à 14%. Nos serviços de pintura 
totalizou-se R$ 347,59 (trezentos e quarenta e sete reais e cinquenta e nove centavos) 
referente à 5%. E nos serviços de instalações elétricas (220V) totalizou-se R$ 180,06 
(cento e oitenta reais e seis centavos) referente à 2%. 

No sistema de vedação vertical interno e externo (paredes) aplicou-se R$ 2.498,61 
(dois mil, quatrocentos e noventa e oito reais e sessenta e um centavos). O item que exigiu 
maior aplicação de recursos foi o único item presente no serviço, ou seja, a alvenaria de 
vedação, representando 100% do valor total do serviço de sistema de vedação vertical 
interno e externo (paredes). Para levantamento do material que precisou ser acrescido 
foi analisado quais paredes precisaram ser acrescentados para criação de shaíts para 
solucionar os problemas de incompatibilidades e esconder as tubulações. 

Na drenagem de águas pluviais aplicou-se R$ 1.266,69 (mil, duzentos e sessenta e 
seis reais e sessenta e nove centavos). O item que exigiu maior aplicação de recursos foi o 
joelho 90 - 100mm, representando 85,54% do valor total do serviço de drenagem de águas 
pluviais. O tubo de PVC 9100mm representou 14,46%. 


Na instalação sanitária aplicou-se R$ 1.128,66 (mil, cento e vinte e oito reais e 
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sessenta e seis centavos). O item que exigiu maior aplicação de recursos foi o joelho PVC 
45º 75mm, representando 41,12% do valor total do serviço de instalações sanitárias. O 
tubo de PVC rígido 75mm, representou 32,37%. O tubo de PVC rígido 5Omm, representou 
23,17%. O joelho PVC 45º 5Omm, representou 1,42%. O joelho PVC 90º 50mm, representou 
1,32%. 

Na instalação hidráulica aplicou-se R$ 1.093,60 (mil e noventa e três reais e 
sessenta centavos). O item que exigiu maior aplicação de recursos foi o joelho 90 soldável 
75mm, representando 42,87% do valor total do serviço de instalações sanitárias. O joelho 
90 soldável 60mm, representou 41,87%. O joelho 90 soldável 50mm representou 5,12%. 
O tubo de PVC soldável de 60mm, representou 5,04%. O tubo de PVC soldável de 75mm, 
representou 3,84%. O tubo de PVC soldável de 50mm, representou 1,26%. 

Nos revestimentos internos e externos aplicou-se R$ 1.075,87 (mil e setenta e 
cinco reais e oitenta e sete centavos). O item que exigiu maior aplicação de recursos foi o 
emboço, representando 42,60% do valor total do serviço de revestimento interno e externo. 
O reboco, representou 37,62%. O chapisco, representou 19,78%. 

Na pintura aplicou-se R$ 347,59 (trezentos e quarenta e sete reais e cinquenta e 
nove centavos). O item que exigiu maior aplicação de recursos foi a pintura, representando 
61,98% do valor total do serviço de pintura. O emassamento, representou 38,02%. 

Nas instalações elétricas (220V) aplicou-se R$ 180,06 (cento e oitenta reais e seis 
centavos). O item que exigiu maior aplicação de recursos foi o condutor de cobre unipolar, 
representando 70,19% do valor total do serviço de instalações elétricas. O eletroduto PVC 
flexivel corrugado reforçado, representou 29,81%. 

O custo total dos serviços acrescidos foi de R$ 7.591,08 (sete mil, quinhentos e 
noventa e um reais e oito centavos), para solucionar os problemas de incompatibilidades 
entre os projetos encontrados durante o processo de execução. O item que maior 
representou no orçamento foi o serviço de sistema de vedação vertical interno e externo 


(paredes). O Quadro 1 apresenta todos os dados de forma sintética referente ao custo total. 
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OBRA: CRECHE PADRÃO FNDE - TIPO 1 - BAIRRO NOVA CARUARU - CARUARU-PE 


PLANILHA RESUMO DE CUSTO 


DATA BASE SINAPI/PE 06/2017 DESONERADA, ORSE 06/2017, EMLURB 12/2014, SEINFRA 


BDI: 27,7 % 


ITEM [CÓDIGO | FONTE DESCRIÇÃO DOS SERVIÇOS uno, a o TOTAL 





SISTEMA DE VEDAÇÃO VERTICAL INTERNO E EXTERNO (PAREDES) 
ALVENARIA DE VEDAÇÃO 


R$ 2.498,61 








SINAPI 





Alvenaria de vedação de 1/2 vez emtijolos cerâmicos (dimensões nominais: 
39x19x09); assentamento em argamassa no traço 1:2:8 (cimento, cal e 
areia) para parede interna 


m? 








2.498,61 








REVESTIMENTOS INTERNOS E EXTERNOS 


R$ 1.075,87 





SINAPI 


SINAPI 


Chapisco de aderência em paredes internas, externas, vigas, platibanda e 
calhas 

Emboço para paredes internas e externas traço 1:2:9 - preparo manual - 
espessura 2,0 cm 


212,82 


458,31 








SINAPI 





Reboco para paredes internas, externas, pórticos, vigas, traço 1:4,5 - 
AsSDO aa 








404,74 





PINTURA 


347,59 





SEINF 


DA 


Emassamento de paredes internas com massa acrílica - 02 demãos 


132,17 








SINAPI 





Pintura em latex acrílico 02 demãos sobre paredes internas, externas 








215,42 








89449 


SINAPI 


INSTALAÇÃO HIDRÁULICA 
TUBULAÇÕES E CONEXÕES DE PVC RÍGIDO 


Tubo PVC soldável 9 50 mm, fornecimento e instalação 


m 


13,10 


R$ 1.093,60 


13,76 





89450 


2 
24 
2a 
23 

3 
3.1 
3.2 
4.2 
4.3 


SINAPI 
SINAPI 


Tubo PVC soldável (9 60 mm, fornecimento e instalação 
Tubo PVC soldável 9 75mm, fornecimento e instalação 


20,05 


55,14 





4.4 


SINAPI 


Joelho 90 soldável - 5Omm, fornecimento e instalação 





4.5 


4.6 | 89521 


SINAPI 
SINAPI 


Joelho 90 soldável - 60mm, fornecimento e instalação 
Joelho 90 soldável - 75mm, fornecimento e instalação 








DRENAGEM DE ÁGUAS PLUVIAIS 
TUBULAÇÕES E CONEXÕES DE PVC 


R$ 1.266,69 





5.1 
5.2 


89848 
89744 


SINAPI 
SINAPI 


Tubo de PVC 9100mm, fornecimento e instalação 
Joelho 90 - 100mm, fornecimento e instalação 


m 9,60 19,08 183,17 


64,00 16,93 1.083,52 


un 








6 


INSTALAÇÃO SANITÁRIA 


R$ 1.128,66 





6.1 89712 


89511 


SINAPI 
SINAPI 


Tubo de PVC rígido 50mm, fornec. e instalação 
Tubo de PVC rígido 75mm, fornec. e instalação 


19,23 
23,88 


268,26 
365,36 





89739 


6.3 
6.4 


SINAPI 
SINAPI 


Joelho PVC 45º 75mm - fornecimento e instalação 
Joelho PVC 45º 50mm - fornecimento e instalação 


13,65 


464,10 





SINAPI 


Joelho PVC 90º 50mm - fornecimento e instalação 








INSTALAÇÕES ELÉTRICAS - 220V 
ELETRODUTOS E ACESSÓRIOS 


R$ 180,06 





SINAPI 


Eletroduto PVC flexível corrugado reforçado, 925mm (DN 1º), inclusive 
conexões 
CABOS E FIOS (CONDUTORES) 





Condutor de cobre unipolar, isolação em PVC/70ºC, camada de proteção 
em PVC, não propagador de chamas, classe de tensão 750V, 
encordoamento classe 5, flexível, com as seguintes seções nominais: 








SINAPI 





2,5 mm? 








126,38 





Quadro 1 — Planilha Sintética 
Fonte: Autor (2020). 
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TOTAL: 
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R$ 7.591,08 





41 CONSIDERAÇÕES FINAIS 


A creche tipo 1 padrão FNDE, é um projeto nacional para desenvolvimento da 
educação infantil no Brasil. No município de Caruaru contempla 5 (cinco) creche deste tipo, 
onde três delas já estão construídas e duas em etapa de construção. Como os projetos 
foram elaborados antes do decreto nº 9.983 de 22 de agosto de 2019 do Governo Federal, 
trata-se de uma elaboração sem a ferramenta BIM, onde deixa algumas incompatibilidades 
entre seus projetos. 

Após identificar e quantificar as incompatibilidades presentes nos projetos estrutural 
e hidrossanitário, arquitetônico e hidrossanitário, arquitetônico e elétrico, foi realizada uma 
análise das soluções adotadas durante a execução para as interferências encontradas 
entre os projetos. Devido a necessidade de acrescentar os serviços para solucionar os 
problemas ocorridos foi feito um levantamento dos itens para quantificar o impacto no 
custo, para atribuir os valores adotados na planilha dos itens acrescidos, foi utilizado como 
base os valores apresentados na planilha orçamentária da própria creche, que sua base é 
SINAPI/PE 06/2017 desonerada, ORSE 06/2017, EMLURB 12/2014 e SEINFRA, adotando 
um BDI de 27,7%. 

Foram identificadas 72 interferências durante a etapa de execução por falta de 
compatibilização entre os projetos, foi feito um levantamento dos itens que precisaram ser 
acrescidos para as soluções adotadas, o valor total desses serviços foi de R$ 7.591,08, 
(sete mil, quinhentos e noventa e um reais e oito centavos). Os projetos que apresentaram 
maiores interferências foram entre o estrutural e hidrossanitário, já o item que demonstrou 
um valor maior foi o serviço de sistema de vedação vertical interno e externo (paredes), 
com um valor de R$ 2.498,61 (dois mil, quatrocentos e noventa e oito reais e sessenta e 
um centavos). 

Para evitar problemas decorridos da falta de compatibilização entre os projetos, 
é indicado a utilização da ferramenta BIM, pois é possível observar com antecedência 
os problemas de incompatibilidade, sendo possível soluciona-los na fase de elaboração, 
potencializando o desempenho na fase de execução e otimizando o custo final do 
empreendimento. 

Com o incentivo do Governo Federal para a utilização da ferramenta BIM nas 
obras públicas, será possível fazer o uso das dimensões (3D ao 7D) da obra para obter 
melhores resultados e evitar problemas durante as etapas de construção e utilização do 
empreendimento. 

A ferramenta BIM vai além da modelagem de um produto ou compatibilização de 
projetos, nele é possível englobar todos os aspectos relativos ao empreendimento, desde 
os produtos e processos, até a documentação. 

Portanto, se os projetos da creche tipo 1 padrão FNDE fossem elaborados utilizando 
o BIM, os problemas ocorridos durante a execução seriam solucionados durante a etapa de 
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elaboração, onde evitaria impactos no custo nacionalmente. 

Propõe-se como sugestão para um trabalho futuro o desenvolvimento de uma 
análise dos impactos no tempo causados pelas incompatibilidades entre projetos da creche 
tipo 1 padrão FNDE. 
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RESUMO: Azeotropia é um problema recorrente 
na indústria para a obtenção de substâncias 
mais puras e/ou concentradas. A mistura líquido- 
líquido binária de água e 1-propanol apresenta 
um ponto de azeotropia quando a mistura contém 
cerca de 42% molar de 1-propanol à 1 atm. Uma 
maneira de contornar esse problema é a adição 
de um agente arrastador que altera as interações 
moleculares e muda o azeótropo. Este trabalho 
avalia os efeitos de 24 diferentes compostos 
arrastadores na mistura água e 1-propanol quanto 
ao rendimento da destilação. A adição de um 
arrastador na mistura possibilita uma separação 
mais eficiente, principalmente se formar um ponto 
de azeotropia heterogênea ternária. O efeito 
desse arrastador foi determinado pelo método de 
contribuição de grupos UNIFAC, utilizando-se da 
abordagem computacional do XSEOS. 
PALAVRAS-CHAVE:  Arrastador, Azeotropia, 
Equilíbrio Líquido-Líquido, Mistura água/1- 
propanol, XSEOS. 


ABSTRACT: Azeotropy is a recurring problem 
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XSEOS 


in the industry for obtaining purer and/or 
concentrated substances. The binary liquid- 
liquid mixture of water and 1-propanol shows an 
azeotropy point when the mixture contains about 
42 % molar of 1-propanol at 1 atm. One way 
around this problem is to add an entraining agent 
that alters molecular interactions and changes the 
azeotrope. This work evaluates the effects of 24 
different entraining compounds in the water and 
1-propanol mixture in terms of distillation yield. 
The addition of an entrainer in the mixture allows 
for a more efficient separation, especially if a 
heterogeneous ternary azeotropy point is formed. 
The effect of this entrainer was determined by the 
UNIFAC group contribution method, using the 
XSEOS computational approach. 

KEYWORDS: Entrainer, Azeotropy, Liquid-Liquid 
Equilibrium, Water/1-propanol mixture, XSEOS. 


INTRODUÇÃO 


Destilação  azeotrópica heterogênea 
é um processo largamente praticado para a 
desidratação de uma vasta gama de materiais, 
incluindo ácido acético, clorofórmio, etanol, e 
muitos álcoois superiores. Segundo Benedict 
e Rubin (1945), o princípio da destilação 
azeotrópica consiste em adicionar uma nova 
substância (arrastador) a uma mistura, de 
maneira a aumentar a volatilidade relativa das 
dois componentes chaves e tornar, assim, 
relativamente fácil a separação. A substância 
adicionada forma um azeótropo com um 


ou mais dos constituintes da mistura e, em 
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consequência, está presente na maioria dos pratos da coluna em concentração apreciável. 

Dadas as não idealidades, a divisão de fases, os limites de destilação presentes, 
e a possível existência de múltiplos estados estacionários em tal sistema, colunas como 
estas podem ser extremamente difíceis para simular e para operar. Uma abordagem 
computacional para o problema pode poupar muitas dores de cabeça (TAYLOR e KRISHNA, 
2000). 

Dentre das várias abordagens computacionais possíveis, destaca-se o XSEOS, 
devido ao seu baixo gasto computacional e sua fácil instalação. O XSEOS é um add-in 
gratuito para Excel para o cálculo de propriedades termodinâmicas. Com mais de 22.000 
linhas de comando, cobre várias áreas da Termodinâmica da Engenharia Química. A sua 
utilização e funcionamento foi descrito por Castier (2008). 

Dos vários modelos termodinâmicos que o XSEOS calcula, o UNIFAC é um dos 
métodos mais utilizados. O método UNIFAC supõe que a mistura líquida é uma solução 
das unidades estruturais (denominadas grupos e subgrupos) que constituem as moléculas 
em vez de considerá-las como uma mistura das próprias moléculas. Desta maneira, as 
propriedades de uma mistura podem ser representadas pela soma das contribuições 
individuais de cada um dos grupos que compõem a mesma (DONG et al., 2020). 


MATERIAIS E MÉTODOS 


Um ponto de azeotropia ocorre quando a fração molar da fase líquida (Xi) é igual 
a fração molar da fase vapor (Yi). O equilíbrio líquido-vapor da mistura 1-propanol e 
àgua, inclusive seu ponto de azeotropia, está representado na Figura 1. Um esquema 
de separação por coluna destilação e decantadores entre água e 1-propanol por meio do 
arrastador éter di-isopropílico está respresentado na Figura 2. 


w 
E | 
oo 


— UNIFAC 
UNIQUAC 
----NRTL 


Temperatura (K) 
Q  Q Q Q o [95] 195] w 
o o o o = = =! = 
Mo + o oo O ] + o 


w 
o 
Õ 





w 
q 
oo 


0.0 01 02 03 04 05 06 07 08 09 10 
Fração molar de 1-propanol (x1,y1) 


Figura 1. Equilíbrio líquido-vapor da mistura 1-propanol e água à 1 atm. 
Fonte: Adaptado de Devi et al. (2017). 
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Reposição de 
arrastador 












Leve2, arrastador reciclado 


Pesado2 








DECANTADOR 







COLUNA 


Água + AZEOTRÓPICA 


1-Propanol 
Pesado2 


(1-Propanol) 
Pesadol : 
Água 


Figura 2. Sequência geral de destilação azeotrópica da mistura água/1-propanol com 
separação com o auxílio do arrastador éter di-isopropílico. 


Fonte: Adaptado de Pla-Franco et al. (2019). 


No caso de um azeótropo heterogêneo, duas fases líquidas ocorrem. As frações 
molares das fases líquidas devem estar relacionadas com a fração molar geral da fase 
líquida de acordo com o seguinte balanço de conservação de mols: 


X o Mi=x00º motta qbeta. mjbeta (1) 
e 
ng] = mjeife + Mgj9ote (2) 


Onde: MI é o número de mols total das fases líquidas, ici é a fração molar de i na 
fase líquida alfa, M L “if? é o número de mols da fase alfa, x”: é a fração molar de i na fase 
beta e M Il”: é o número de mols da fase beta. 

No azeótropo heterogêneo, as fases líquidas estão em equilíbrio entre si e também 


com a fase vapor, logo as seguintes condições devem ser satisfeitas: 


— alfa cc qmfa nset — bate, q bata, 
Ra F. PÉ =x, 


XY, a Y; E (3) 
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Onde: P é a pressão total do sistema, Va é o coeficiente de atividade de i na 
fase líquida alfa, P*" é a pressão de saturação e i na temperatura do sistema e à di é 
o coeficiente de atividade de i na fase líquida beta. Os coeficientes de atividade foram 
calculados utilizando o modelo UNIFAC. As pressões de saturação foram calculadas de 


acordo com as equações presentes no software DIADEM. 


B 
Pa = exp(4A+>+ Cn? + DTS) 
F (4) 


Onde a temperatura (T) é dada em Kelvin (K) e Psat é dada em Pascal (Pa). As 
constantes A, B, C, De E variam para cada substância e foram retiradas da base de dados 
do DIADEM. 

Para que todas as condições citadas acima fossem satisfeitas, a minimização 
da seguinte função objetivo foi realizada no programa Excel para Windows, utilizando a 


ferramenta “Solver”: 


FO = Errol+ Erros + Erros (5) 
Sendo: 
NC 
Errol = > — y;º 
f=1 (0) 


NC 
Erro? = >. (x, «MI — E radio - MjUÍas ao preco) 
1 (7) 


NE 
ê 
re > (apre É pote = are l add 
t=1 (8) 


Onde FO é a função objetivo, Erro1 é o erro associado à condição de azeotropia, 
Erro2 é o erro associado a condição de conservação de número de moles das fases líquidas 
e Erro3 é o erro associado à condição de equilíbrio entre as duas fases líquidas. 

Durante o processo de minimização, as frações molares gerais da fase líquida, as 
frações molares de cada fase líquida, o número de mols da fase alfa, e a temperatura do 
sistema foram considerados as variáveis de decisão (ou seja, foram as variáveis ajustadas 
para minimizar a função objetivo). As frações molares da fase vapor e o número de mols da 
fase líquida beta foram calculados da seguinte forma: 


GA E pan 

E E O 
Yi p (9) 
ngjºste =1 — mjetfe (10) 


O valor de 1 na equação acima foi a base de cálculo escolhida para o número total 
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de moles nas fases líquidas. 

É necessário, primeiramente, o teste da metodologia em vários equilíbrios 
azeotropos ternários, os dados encontrados na literatura foram do equilíbrio água- 
propanol-arrastador. Na Tabela 1 temos os dados retirados de Dimian et al. (2004) para 11 
diferentes arrastadores. 


Azeótropo heterogêneo ternário Temperatura do 

Arrastador É st oC 
0,1900 0,0940 0,7160 59,23 
0,3810 0,2100 0,4090 75,30 
0,4850 0,2910 0,2240 81,40 
0,5655 0,3170 0,1175 85,40 


Tabela 1. Dados experimentais retidos da literatura das frações molares dos equilíbrios 
ternários Água+1-Propanol+Arrastador. 





Fonte: Baseado em Dimian et al. (2004). 


RESULTADOS E DISCUSSÕES 


Na Tabela 2 é apresentado os dados calculados utilizando-se do XSEOS para os 
mesmos 11 arrastadores da Tabela 1. 


Azeótropo heterogêneo ternário (% Erro) 
Temperatura do 


Arrastador azeótropo (ºC) 


Água 1-Propanol Arrastador 


RR 0,1873 0,1008 0,7119 59,92 
(1,42%) (7,23%) (0,57%) (1,16%) 


esti 0,3584 0,2015 0,4400 74,65 
p (5,93%) (4,05%) (7,58%) (0,86%) 
0,4805 0,2886 0,2308 82/11 


(0,93%) (0,82%) (3,04%) (0,87%) 


0,5449 0,3440 0,1111 85,59 
(3,64%) (8,52%) (5,45%) (0,22%) 


Octano 


Nonano 
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2,2 ,4-trimetilpentano 


1-hexeno 


1-hepteno 


1-octeno 


Ciclohexano 


Ciclohexeno 


1,3-Ciclo-hexadieno 


Média + DP 


0,3543 
(19,48%) 


0,1662 
(191,58%) 


0,3393 
(1,89%) 
0,4704 
(2,04%) 
0,2566 
(12,60%) 
0,2746 
(4,98%) 


0,2634 
(3,52%) 


0,2063 
(0,82%) 


0,0826 
(46,01%) 
0,1872 
(1,99%) 
0,2768 
(3,89%) 
0,1343 
(29,63%) 
0,1400 
(21,74%) 


0,1328 
(27,69%) 


0,4394 
(24,83%) 


0,7512 
(4,91%) 
0,4735 
(0,53%) 
0,2528 
(0,72%) 
0,6091 

(1,05%) 
0,5853 
(1,80%) 


0,6038 
(3,08%) 


22,00+56,34 13,85x15,12 4,87+7,01 


74,53 
(0,87%) 


57,00 
(0,12%) 
73,26 
(0,11%) 
81,45 
(0,15%) 
66,72 
(0,26%) 
68,15 
(0,51%) 


67,16 
(0,87%) 


0,55+0,39 





Tabela 2. Valores calculados pelo XSEOS das frações molares dos equilíbrios ternários 


Água+1-Propanol+Arrastador. 


Com base nos dados da Tabela 1, foi calculado o erro relativo de cada mistura 


ternária na Tabela 2, em que os valores dos erros estão entre parênteses. Além disso, 


calculou-se o desvio padrão (DP) dos erros. Apesar do 1-hexeno apresentar um alto erro 


relativo, ele apresenta um baixo erro absoluto. Além de que essa porcentagem mais elevada 


de erro se deve por o valor real se aproximar de zero, sendo que todo denominador que 


tende o valor para zero faz elevar o valor da divisão. Devido aos baixos valores dos erros, 


comprovou-se a eficiência da metodologia. 


Estimou-se os valores das composições globais e a temperatura do ponto azeotrópico 


para 24 diferentes arrastadores alcanos e alcenos, determinando se há ou não azeotropia e 


qual seria o tipo, se binária ou ternária. Os valores estimados estão na tabela 3. 


[ata (E que] 
Arrastador 
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Água 
0,0001 
0,0003 
0,0003 
0,0035 
0,0538 
0,0432 
0,3166 
0,3340 
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Temperatura 
do azeótropo 
(ºC) 


Tipo de 
Azeotropia 





2-metilhexano 
3-metilhexano 
Decano 
1-propeno 
1-buteno 
1-penteno 
1-noneno 
1-deceno 
Ciclopropano 
Ciclobutano 
Ciclopentano 
Cicloheptano 
Ciclooctano 
Ciclopenteno 
Ciclohepteno 


Cicloocteno 


Tabela 3. Valores estimados via XSEOS das frações molares dos equilíbrios ternários Água+1- 


0,4892 
0,4980 
0,7893 
0,0006 
0,0145 
0,1342 
0,7308 
0,9819 
0,0005 
0,0146 
0,0967 
0,0000 
0,0060 
0,0822 
0,4445 


0,5779 


0000 | oz | ox 
0003 | com | cos 


0,0003 0,9991 -32,66 


0,0003 0,8656 26,70 


ooo | os | 1244 
pos | os | as 


osso | oro | org 
posto | osso | azoa 
254 | os | zag 


0,4215 0,0005 88,24 


Propanol+Arrastador. 


0,6362 0,3638 92,09 


Binária 





Na Tabela 3 pode-se ver que dos 24 arrastadores: 6 não apresentam nenhuma 


azeotropia, 13 apresentam azeotropia binária e 5 apresentam azeotropia ternária. 


CONCLUSÕES 


A ferramenta XSEOS se mostra de fácil implementação e de baixa demanda 
computacional. Além disso, a metodologia aqui empregada se mostrou uma ferramenta útil 
na estimativa de equilíbrios azeotrópicos não disponíveis na literatura ou de difícil obtenção 
experimental, devido a ser uma ferramenta gratuita e os seus resultados serem próximos 


dos da literatura. 
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RESUMO: Muitas indústrias produzem elementos 
estruturais constituídos de dois ou mais materiais 
diferentes, buscando assim, a capacidade 
destes de suportar de modo mais eficaz cargas 
aplicadas. Esses elementos são chamados de 
materiais compostos e um exemplo destes são 
as barras compostas. Devido à importância 
dessas barras, este trabalho traz uma análise 
do comportamento mecânico destes materiais 
quando são flexionados devido a momentos 
fletores aplicados em planos dos eixos de 
simetria e fora destes. Para se obter uma análise 
mais eficiente do comportamento dessas barras 
tornou-se necessário o uso do software ANSYS, 
que utiliza o método dos elementos finitos 
(MEF) para a solução de diversos problemas 
em engenharia. Assim, com o uso do ANSYS, 
foi possível simular o comportamento mecânico, 
tensões e deformações atuantes, da barra com 
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boa precisão ao longo de todo o seu comprimento 
quando está fora flexionada devido a momentos 
fletores aplicados. 

PALAVRAS-CHAVE: Elementos finitos; ANSYS; 
Flexão; Barras. 


FINITE ELEMENTS DESIGN: BENDING IN 
BARS WITH DIFFERENT MATERIALS 


ABSTRACT: Many industries produce structural 
elements consisting of two or more different 
materials, thus seeking their ability to more 
effectively support applied loads. These elements 
are called composite materials and an example 
of these are the composite bars. Due to the 
importance of these bars, this work brings an 
analysis of the mechanical behavior of these 
materials when they are flexed due to bending 
moments applied in planes of the symmetry axes 
and outside of them. In order to obtain a more 
efficient analysis of the behavior of these bars, 
it was necessary to use the ANSYS software, 
which uses the finite element method (MEF) to 
solve various engineering problems. Thus, with 
the use of ANSY'S, it was possible to simulate the 
mechanical behavior, stresses and deformations 
acting, of the bar with good precision along its 
entire length when it is flexed out due to applied 
bending moments. 

KEYWORDS: Finite elements; ANSYS; Flexion; 
Bars. 


11 INTRODUÇÃO 


Com o intenso crescimento das cidades 


em todo o mundo aumenta, também, como 
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consequência, o número de estruturas de grande porte que são empregadas em edifícios, 
pontes, etc. Para isso, faz-se necessário a utilização de materiais que possam garantir uma 
integridade estrutural eficiente, a fim de que, possam ser evitados prejuízos econômicos 
e vidas não sejam perdidas. Pensando nisso, muitas indústrias produzem elementos 
estruturais constituídos de dois ou mais materiais diferentes, aumentando assim, a 
capacidade destes de suportar de modo mais eficaz cargas aplicadas. Esses elementos 
são chamados de materiais compostos e um exemplo destes são as barras compostas. 
Devido à importância dessas barras, este trabalho traz uma análise do comportamento 
mecânico destes materiais quando são flexionados devido a momentos fletores aplicados 
em planos dos eixos de simetria e fora destes. 

Para se obter uma análise mais eficiente do comportamento dessas barras tornou- 
se necessário o uso do software ANSYS, que utiliza o método dos elementos finitos (MEF) 
para a solução de diversos problemas em engenharia. O método dos elementos finitos é 
uma importante ferramenta computacional que permite realizar uma grande quantidade 
de operações matemáticas que se fossem feitos manualmente gastariam muito tempo e 
seriam muitos difíceis. De maneira geral, a ideia básica do MEF é dividir o corpo, domínio, 
em número finito de elementos, subdomínios, com determinadas formas geométricas, ou 
seja, este método divide um problema grande em outros menores que permite com que o 
computador possa realizar vários cálculos na solução dos problemas menores e somar tudo 
no final. Os subdomínios recebem o nome de elementos finitos que estão conectados entre 
si por um número finito de pontos nodais, nós, e, ao conjunto de todos esses itens, dá-se o 
nome de malha. A melhor precisão dos resultados depende, por sua vez, da quantidade de 
nós e elementos, do tamanho e dos tipos de elementos da malha. 

Portanto, com o uso do ANSYS, foi possível simular o comportamento mecânico, 
tensões e deformações atuantes, da barra com boa precisão ao longo de todo o seu 
comprimento quando está fora flexionada devido a momentos fletores aplicados. 


21 FUNDAMENTAÇÃO TEÓRICA 


Os elementos de viga são elementos lineares que apresentam uma das dimensões 
muito maior do que as outras duas (dimensões da seção transversal) e que está submetido 
a cargas perpendiculares ao seu eixo longitudinal. Nesses elementos, desenvolvem-se em 
suas seções transversais momento fletor e esforço cortante, sendo o primeiro responsável 
pela flexão e o segundo pelo cisalhamento da viga. O esforço cortante tem na maioria das 
vezes uma influência desprezível no comportamento da peça e podemos, com a finalidade 
acadêmica, desprezá-lo, estudando o efeito apenas da flexão isolada. Assim podemos 
classificar a flexão em: pura e simples. 

É necessário analisar também a posição do carregamento em relação à posição da 
seção transversal da peça. Convencionando por x e y os eixos principais centrais de inércia 
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da seção transversal da viga e aplicando carregamentos das seguintes formas: 





Figura 1. Carregamento Figura 2. Carregamento Figura 3. Carregamento 
no eixo y no eixo x fora dos eixos principais 


Analisando as situações, podemos classificar a flexão em: reta, a qual ocorre quando 
o carregamento está atuando em um dos eixos principais centrais de inércia da seção (x 
ou y). E oblíqua, a qual ocorre quando o carregamento é desviado em relação aos eixos 
principais centrais de inércia da seção. 


2.1 Barra sob Flexão Pura 


Considerando uma barra prismática AB possuindo um plano de simetria e submetida 
a conjugados iguais e opostos M e M' atuando naquele plano podemos observar que, se 
uma seção da barra AB for cortada em um ponto qualquer C, as condições de equilíbrio 
da parte AC da barra requerem que os esforços internos na seção sejam equivalentes ao 
conjugado M. O momento M daquele conjugado é chamado de momento fletor na seção. 





Figura 4. Barra sob conjugado M e M”. 


Chamando de o, a tensão normal em um ponto da seção transversale T, e T, as 
componentes da tensão de cisalhamento, expressamos que o sistema das forças internas 


elementares que atuam na seção é equivalente ao momento fletor M. 
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Figura 5. Representação do equivalente de forças. 


Da estática sabemos que a soma das componentes de forças em qualquer direção 
é igual a zero. Além disso, o momento fletor é o mesmo em relação a qualquer eixo 


perpendicular a seu plano, e é zero em relação a qualquer eixo contido naquele plano. 


SFx=0.:[foxdA=0 (Eq. 1) 
*My=0.:jfzoxdA=0 (Eq. 2) 
»*Mz=0.:J(-yoxdA) =0 (Eq. 3) 


2.2 Deformações sob Flexão Pura 


Analisando um trecho de viga submetida a flexão pura (figura 6) , o momento fletor é 
constante para qualquer seção transversal e a viga se flexiona de maneira uniforme. Assim, 
a linha AB na face superior da viga tem uma curvatura constante, o mesmo acontece com 


a linha A'B”, na face inferior da viga. 


Los 
LARA SER 
EEE A EO “xa “ 





Figura 6. Figura 7. Seções verticais e horizontais longitudinais. 
Representação da viga 
sob flexão pura. 


Supondo que a viga seja dividida em um grande número de cubos elementares com 
faces paralelas aos eixos de coordenadas. Todas as faces estão a 90º uma da outra. Assim 
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concluímos que na superfície 0,, 0,e T,. São zero. 

Logo, a Única componente de tensão que não se anula é a componente 0.. Desse 
modo em qualquer ponto de uma barra submetida à flexão pura, tem-se um estado uniaxial 
de tensões. 

Para M>0, observa-se que as linhas ABe A'B'diminuem e aumentam em comprimento 
respectivamente e nota-se que a deformação específica ex e a tensão ox são negativas 
na parte superior da viga (compressão) e positivas na parte inferior (tração). Apenas na 
chamada superfície neutra ex e ox são zero, pois essa superfície intercepta o plano de 
simetria, é paralelo às faces superior e inferior. 

A deformação normal específica é avaliada como a razão entre a variação de 


comprimento e o seu comprimento original, resultando em: 
EX = no (Eq. 4) 


O sinal negativo indica que a deformação é de compressão para um ponto y positivo 


(acima da linha neutra). 


2.3 Tensões e Deformações no Regime Elástico 


Consideramos que o momento fletor M seja de forma que as tensões normais na 
viga permanecem abaixo da tensão de escoamento do material o. Não haverá deformação 
permanente, e vale a lei de Hooke para tensão uniaxial. Considerando que o material seja 
homogêneo, e chamando de E seu módulo de elasticidade, temos na direção longitudinal 
x: 0, = E€.. 

Portanto a distribuição de tensões normais está diretamente associada à distribuição 
de deformações normais ex. Assim podemos dizer que ox varia linearmente com a distância 


y da linha neutra, ao longo da seção transversal. Logo: 


ni Es (Eq. 5) 
p 


Compressão 
i 





LN Y Linha Neutra 
(nem tração 
nem compressão) 
Tração 


Figura 8. Distribuição das tensões normais variando com a distância y. 


A partir das relações constitutivas, para um estado uniaxial em um material linear 


elástico e do equilíbrio estático chegamos em: 
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My (Eq. 6) 


O. z 
* I 
Onde M é o momento fletor, y é a distância do ponto de análise até a linha neutra e 
[o momento de inércia de área da seção transversal. 


2.4 Flexões em Barras com Diferentes Materiais 


As vigas são muitas vezes constituídas por mais que um material constituindo as 
chamadas vigas não homogêneas. Alguns exemplos de vigas deste tipo são vigas de 
madeira reforçadas a aço, viga constituída de dois materiais metálicos. A teoria da flexão 
de vigas sujeitas a momentos fletores pode ser facilmente adaptada ao estudo de vigas 
constituídas por dois ou mais materiais. 

Considerando uma barra formada por dois materiais diferentes unidos. Essa barra 
composta se deformará conforme já mencionado, pois sua seção transversal permanece 
a mesma em todo o comprimento. Assim, a deformação específica normal ox ainda varia 
linearmente com a distância y da linha neutra da seção. 


y (Eq. 7) 


p 





(a) (b) (c) 


Figura 9. Distribuição de deformação específica. 


No entanto, não podemos supor que a linha neutra passe pelo centróide da seção 
composta. Como os módulos de elasticidade E1 e E2 dos dois materiais são diferentes, 
as expressões obtidas para a tensão normal em cada material também serão diferentes. 
Considerando E2>E'1: 


Pr Bite É1y (Eq. 8) 
1— LbICx— 
p 
Eq. 9 
2 2€x p 


Obtemos uma curva de distribuição de tensões consistindo em dois segmentos de 
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reta. Conclui-se das equações acima que a forças dF que atuam no elementos de área dA 


da são: 

dr1 — ojdA — -*2 dA (Eq. 10) 
E (Eq. 11) 

dp> — o,dA — = dA i 


Chamando de n a relação E2/E1, temos: 


E (Eq. 12) 
dr; = ——— dA = — É (nda) 
p 
Nota-se que a mesma força dF, atuaria em um elemento de área ndA do primeiro 


material, assim: 


0/= 0, e 0, =n0, 


A análise para situações como esta deve ser da seguinte forma: 

Primeiro definimos um material de referência (de preferência, o de menor módulo 
de elasticidade) e modificamos as áreas dos outros materiais (em dimensão na direção 
paralela à linha neutra. A modificação é feita a partir da relação do módulo de elasticidade 
do material de cada parte com o do material de referência. 

Depois obtemos o momento de inércia | da área total transformada e calculamos 
as tensões atuantes. E finalmente multiplicam-se as tensões calculadas para a área 
transformada pelas razões dos módulos de elasticidade n para as áreas com materiais 
diferentes do material de referência. 

Exercício Teórico: A viga composta é feita de aço (A) e latão (B) e tem a seção 
transversal mostrada na figura. Se for submetida a um momento M = 5 kN.m, determine a 


tensão máxima no latão e no aço. 





Primeiro transformamos a barra em uma seção equivalente inteira de latão. 
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Eaço | 200GPa 
Elat — 105GPa 





= 1,904 assm: b'=nxbh=1,904x 150 mm = 285,6 mm 


285,6 mm 
——————————»+ 


50 mm 


50 mm 


150 mm 


Cálculo do centróide: 


= LyA o (150x50)x25 + (285,6x50)x75 ao 
VS SA” ” (150x50) + (285,6x50) 2! 62 mm 


Cálculo do momento de inércia da seção transformada: 
1= X(T+ Ad?) = =(5 bhº + Ad?) 
I= [5 (0,150x0,053) + (0,150x0,05)x(0,032782º) | + 
+ [> (0,02856x0,05º) + (0,2856x0,05)x(0,017218?) = 16,8308x10º m? 


Cálculo das tensões máximas: 


My | 5000x0,057782 


O máx lat — T = “16,8308x10-6 = 17,165 Mpa (trativa) 


My 5000x(0,1 — 0,057782) 


= n— = 1,904 = 23,879M É 
Omixaço = N T 904 x 16,8308x1075 3,879 Mpa (compressiva) 


31 MODELO COMPUTACIONAL 


3.1 Pré-processamento 

Para o modelo proposto foram utilizados parâmetros geométricos coerentes com 
a seção proposta na fundamentação teórica, bem como com a compatibilidade entre os 
carregamentos do modelo computacional e teórico. Foi obtida assim uma viga de seção 
retangular (150x100mm) e comprimento 980mm, podendo ser observada na figura 10. 
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Figura 10. Modelo 3D de uma viaga de seção retangular de dois materiais. 


Para os materias foram selecionados o Aço estrutural e o Latão predefinidos na 
biblioteca de materiais do Ansys, com as seguintes propiedades: 


Materiais Módulo de elasticidade[GPa]  Densidade[kg/m?] 
Aço estrutural 7850 
Latão 8300 





Tabela 1 


Primeiramente foi gerada uma malha de elementos quadriláteros sem qualquer 
com tamanho de elementos 50mm, contudo prosseguiu-se a simulação para caráter 
de comparação com os resultados teóricos e, também, confirmação da necessidade de 
refinamento. Verificou-se valores incoerrentes com aqueles esperrados pela teoria. Para 
corrigir o erro, diminuiu-se o tamanho dos elementos para 30mm, e para maior aproximação 
com o resultado teórico, foi reduzido ao tamanho final de 20mm. Através disso, foi possível 
obter novos resultados, com maior nível de confiabilidade. A malha final pode ser observada 
na figura 11. 
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0,00 200,00 (mm) o 
Tm "ss z X 


100,00 
Figura 11. Malha refinada final, composta de 13028 nós e 2352 elementos. 


As arestas frontal e traseira da viga foram fixadas (simulando apoios simples), e 
como na resolução analítica foi aplicado um carregamento tal que o momento fletor gerado 
fosse de 5 kN“*m. Para o cálculo do carregamento, admitiu-se uma carga P distribuída 
na face superior da viga. A mesma gera um momento interno cuja magnitude é dada por 
M = -P * L?, para o comprimento utilizado, e em conformidade com a fundamentação 
teórica, temos um carregamento de, aproximadamente, 41649N/m. O sistema está disposto 
na figura 12. 


ANSYS 


R18.0 
Academic 


0,00 200,00 (mm) sm 
To "3 z x 


100,00 


Figura 12. Esquema de aplicação da força e engaste 
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3.2 Solução e Pós-processamento 
As soluções levaram em média 18 segundos para serem obtidas, apresentando os 


seguintes resultados em ordem de tamanho de elemento: 


Tamanho de Tensão Máxima(Latão)[MPa] Tensão Máxima(Aço)[MPa] 
elemento 


Default(50mm) 11,544 -20,849 
30Omm 15,486 -22,617 
20mm NOS -24,481 


Tabela 2 





Para o resultado mais próximo do teórico obteve-se o seguinte mapa de tensões: 


ANSYS 


R18.0 
Academic 





“ 
0,00 400,00 (mm) ea 
Tm 


200,00 


Figura 13. Mapa de tensões normais ao longo da barra na direção z 


Para avaliação da deformação ocorrida na viga, é observado que a deflexão da 
mesma no apoio é O e máxima no centro, como prevê a teoria para essa configuração, e 


que pode ser verificada com o mapa deformação da figura 14. 
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Figura 14. Mapa de deformação total 


41 DISCUSSÕES 


4.1 Resultados analíticos e numéricos 


Tensão máxima(Latão)[MPal] 


17,165 





Tabela 3. Resultados analíticos 


Tamanho do elemento 


Tabela 4. Resultados numéricos 
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Gráfico 1. Representação das tensões analíticas e numéricas. 


4.2 Diferenças percentuais 


Para o latão as diferenças percentuais entre o modelo teórico e os resultados 
computacionais obtidos são: 


17,165 17,587 





Difrença [4%] 2,458 


Tabela 5. Diferença percentual entre tensão teórica e numéricas do latão. 


Para o aço as diferenças percentuais entre o modelo teórico e os resultados 
computacionais obtidos são: 


Teórico 


23,819 





Diferença [%] 


Tabela 6. Diferença percentual entre tensão teórica e numéricas do aço. 


O motivo dessas diferenças reside no fato de que o método de elementos finitos, 
o qual os programas computacionais trabalham, ser uma aproximação do método teórico. 
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O termo “discretização” é usado justamente porque passamos de um contínuo, 
modelo teórico, (a função originalu(x) está definida num domínio que é uma reunião infinita 
de pontos), para um conjunto discreto, o qual o domínio passa a ser uma reunião finita de 
intervalos, modelo computacional. Em cada um desses intervalos (retas), aproximamos a 
função original por um segmento de reta. Notadamente, quanto menor o comprimento dos 
subintervalos, ou seja, quanto menor a norma da partição, mais a função discretizadaud se 
aproximará da original u(x), é observável isso na figura abaixo. 

Estes subintervalos descritos acima são os segmentos de reta que tendem 
se aproximarem de u(x) são os tamanhos dos elementos descritos no tópico modelo 
computacional. A menor partição utilizada foi de 20mm, este que tem uma proximidade 


com a teoria em torno de 2,5%. 


u(x) 
Solução exata 


j Aproximação por 4 retas 


Aproximação por uma reta 





Figura 15. Aproximação por funções contínuas por partes 


5 | CONCLUSÃO 


No trabalho foi realizado o cálculo das tensões sofridas pelos diferentes 
materiaiscausada por um momento de 5 kN.m, aproximadamente, pelo método teórico e 
computacional, para percebermos as disparidades entre a teoria e o método de elementos 
finitos. E foi observável a propensão do valor teórico pelo método dos elementos finitos, no 
qual, quanto menor os pontos nodais menor é a diferença entre eles, ou seja, a geometria e 
o refino da malha é um fator determinante para a satisfatoriedade do resultado encontrado 
no ANSYS. 

Os pontos nodais utilizados foram por ordem de refino 50mm, 30mm e 20mm, para o 
primeiro ponto a diferença foi grande, principalmente para o latão chegando a 32,747% que 
sofreu tração modificando de forma desmedida os pontos nodais corpulentos, causando 
esta defasagem entre os valores. No entanto, para a menor partição de 20mm a variação 
entre o método computacional e o príncipio analítico foi em torno de 2,5% para ambos 
materiais, sendo este um valor aceitável para o tempo de processamento, em torno de 18 


segundos. 
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Na hipótese de diminuir de modo desmedido as partições, o tempo do procedimento 
iria aumentar demasiadamente, podendo diminuir a eficiência do trabalho, por uma maior 
precisão, não muito diferente da obtida, por um tempo excessivamente maior que até 
poderia danificar o computador. 

Para trabalhos futuros, seria interessante estabelecer a comparação dos resultados 
obtidos com outras metodologias de análise, como por exemplo cálculos analíticos e/ou em 
alguma viga real com dimensões e propriedades parecidas. Estas comparações serveriam 
como autenticação dos valores encontrados pelo software, viabilizando o estabelecimento 
de diferenças percentuais aceitáveis até para problemas mais complexos. 
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RESUMO: Partículas lignocelulósicas se 
destacam pelo baixo custo, disponibilidade, e 
menor toxidade. Espécies vegetais apresentam 
potencial para o reforço de compósitos 
estruturais. Diante disso, o objetivo deste 
estudo é a caracterização física e mecânica de 
compósitos de matriz de gesso reforçados com 
capim. Para tal propósito foram coletados, secos 
e peneirados, capins advindos de plantio na 
Universidade Federal de Lavras (UFLA). Esses 
compósitos foram confeccionados utilizando- 
se uma relação água / gesso de 0,6 e 5 % em 
massa de capim. Após o preparo e secagem, 
os painéis foram retirados do molde e cortados 
nas dimensões normatizadas para os ensaios 
físicos de densidade, umidade, inchamento e 
absorção de água e para os ensaios mecânicos 
de compressão e flexão estática de três pontos. 
Os resultados demonstraram a massa específica 
encontrada dos corpos de prova compósitos 
foi inferior à de literaturas comparativas, o que 
pode indicar uma boa classificação nos índices 
de mérito. Os valores de umidade, inchamento 
e a absorção de água obtidos também foram 
baixos, significando que houve diminuição da 
higroscopia, ou capacidade de absorver água 
do gesso, o que é desejado na maioria das 
aplicações de materiais multifásicos. Os dados 
obtidos dos ensaios mecânicos evidenciaram 
uma melhora nas propriedades de resistência 
do gesso a partir da adição de capim, um outro 
fator positivo. Portanto, os compósitos de matriz 
de gesso reforçados com capim representaram 
melhora tanto nas propriedades físicas quanto 
mecânicas, o que pode validar uma maior gama 
de aplicações para tais materiais. 
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EVALUATION OF PHYSICAL AND MECHANICAL PROPERTIES OF PLASTIC 
MATRIX REINFORCED WITH GRASS 


ABSTRACT: Lignocellulosic particles stand out for their low cost, availability, and lower 
toxicity. Plant species have the potential to reinforce structural composites. Therefore, the 
objective of this study is the physical and mechanical characterization of plaster matrix 
composites reinforced with grass. For this purpose, grasses collected from planting at the 
Federal University of Lavras (UFLA) were collected, dried and sieved. These composites 
were made using a water / gypsum ratio of 0.6 and 5% by weight of grass. After preparation 
and drying, the panels were removed from the mold and cut to standard dimensions for the 
physical tests of density, moisture, swelling and water absorption and for the mechanical 
tests of compression and static flexion of three points. The results showed that the specific 
mass found in composite specimens was lower than that of comparative literature, which 
may indicate a good classification in the merit indexes. The values of humidity, swelling and 
water absorption obtained were also low, meaning that there was a decrease in hygroscopy, 
or the ability to absorb water from the plaster, which is desired in most applications of 
multiphase materials. The data obtained from the mechanical tests showed an improvement 
in the resistance properties of the plaster from the addition of grass, another positive factor. 
Therefore, plaster matrix composites reinforced with grass represented improvement in both 
physical and mechanical properties, which can validate a wider range of applications for such 
materials. 

KEYWORDS: Vegetable fibers, Static flexion, Sustainability. 


INTRODUÇÃO 


Os sistemas construtivos evoluem a cada dia, principalmente, os que são 
ecologicamente corretos. As pesquisas científicas voltadas para a criação novos materiais 
de origem natural, para aplicação na construção civil, apresentam grande destaque do 
ponto de vista técnico, econômico e social. Tal fato é explicado pela insustentabilidade da 
construção civil, sendo esta responsável por gerar grande parte da produção de lixo no 
mundo. 

Um dos elementos de tal indústria é o gesso, produto cujo uso vem crescendo 
continuamente em todos os tipos de edificações, principalmente no sistema de vedações 
internas drywall (TENÓRIO, 2018). Entretanto, diante da grande produção de resíduos 
advindos do processamento do gesso, tornam-se necessários estudos para diminuição de 
uso, no mínimo parcial, deste material. 

A utilização de fibras vegetais como reforço em matrizes de gesso se apresenta, 
neste contexto, como uma solução sustentável, visto que são renováveis, ambientalmente 
corretas e podem melhorar as propriedades mecânicas, como a resistência à tração, 
a flexão e ao impacto desses materiais. Além do mais, possibilitam a alteração do seu 


comportamento desses compósitos após fissuração, diminuindo os efeitos de uma ruptura 
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brusca (SILVA et al, 2012). 

Diante do exposto, o objetivo desta pesquisa é a caracterização física e mecânica de 
compósitos de matriz de gesso reforçados com capim para comparação desses resultados 
com a literatura e validação da gama de aplicação de tais compostos. 


MATERIAL E MÉTODOS 


As fibras de capim utilizadas neste trabalho foram coletadas na Universidade Federal 
de Lavras; enquanto que o gesso foi adquirido no comércio de Lavras-MG. 

Para o tratamento do capim, inicialmente ocorreu a secagem (figura 01) e em 
seguida a cardagem manual para obtenção das fibras, permitindo a remoção das impurezas 
visíveis, decorrentes do processo de desfibramento. As fibras de capim passaram, ainda, 
por um processo de remoção das impurezas visíveis por meio da cardagem manual para 
proporcionar aderência à matriz quando da fabricação do compósito. 





Figura 01- Secagem do Capim 


A compactação dessas fibras foi realizada via Úmida. Neste processo foi adicionado 
água à medida em que as fibras ficassem umedecidas. Cada prensagem permaneceu 
durante o tempo de 23 minutos para verificar o seu comportamento de compactação. 

O gesso utilizado foi um pó branco, de elevada finura, equivalente ao cimento, 
elevada plasticidade da pasta, pega (aderência) e endurecimento rápido, pequeno poder 
de retração na secagem e estabilidade volumétrica que garantem desempenho satisfatório. 

Para o preparo dos compósitos de gesso e capim, foi utilizada a razão água/gesso 
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de 0,6, e capim 5 %. O gesso foi polvilhado sobre a água por um minuto, a suspensão do pó 
na água fi mantida em repouso por dois minutos e a pasta foi misturada durante um minuto. 

Na figura 2, mostra que após 2 dias de secagem, os painéis foram retirados do 
molde e em seguida cortados utilizando-se uma serra circular nas dimensões normatizada 


para os ensaios físicos e mecânicos. 





Figura 02- Corpos de prova Cortados. 


Para a avaliação da densidade aparente, mediu-se o comprimento, largura e 
espessura dos três corpos de prova. Assim foi obtido o volume dos corpos de prova. A 
densidade foi obtida pela razão entre massa, medida com uma balança de precisão, e o 
volume. 

Já a absorção de água e o inchamento foram encontrados mergulhando-se parte do 
corpo de prova em água e registrando-se suas massas e espessuras após 2 e 24 horas de 
submersão. 

Também foi aferida a propriedade física de umidade nesta caracterização, aferindo- 
se a massa dos materiais antes e após procedimento térmico de 24 horas em uma estufa. 

Os testes mecânicos, entretanto, foram feitos em máquina universal de ensaios 
do Departamento de Ciências Florestais (DCF) da UFLA, com escala de 600 Kg como é 
apresentado na figura 08. 
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Figura 03- Maquina universal de ensaio sendo teste de compressão e flexão. 


O ensaio de flexão estática de três pontos foi realizado conforme prescrições da norma 
ASTM D 790-17. Foram ensaiadas quatro amostras com dimensões de 20,/2x76,53x23,/2 
cm, sendo o vão do ensaio de 10 cm. Os corpos de provas foram posicionados de forma a 
permitir que a carga fosse aplicada no centro dos mesmos. 

Já o ensaio de compressão foi feito conforme as normas ASTM C297/ C297M-16. 
Foram ensaiadas duas amostras, sendo a carga aplicada em uma taxa de movimento 
constante até sofrer a ruptura. 


RESULTADOS E DISCUSSÃO 


Propriedades Físicas 


Os resultados das propriedades físicas aferidas neste projeto estão expostos na 
Tabela 1. 


U (%) AA 24h (%) | (%) p (g/em?) 
13,15 (+ 0,50) 36,00 (+ 1,50) 37,33 (+ 2,00) 1,25 (+ 0,08) 
U - Umidade na base seca; 


AA 24h - Absorção de água após 24 horas de submersão; 
| - Inchamento em espessura; 


p - Densidade aparente. 


Tabela 1. Resultados médios das propriedades físicas encontradas nesta pesquisa seguidos de seus 
respectivos desvios padrões. 
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Com relação à massa específica, o valor encontrado nesta pesquisa foi inferior ao 
encontrado por Oliveira et al. (2012), que em sua caracterização de compósitos de gesso 
com resíduos de etileno acetato de vinila e vermiculita encontraram uma massa específica 
média desse material de 2,60 g/cmº. 

Na maioria das aplicações de compósitos é requerido uma baixa massa específica 
associada à uma alta resistência, e o resultado mensurado pode indicar uma vantagem 
dentro desse contexto. Diferentes materiais são comumente classificados de acordo com 
essas propriedades nos índices de mérito. 

Já os valores de umidade, inchamento e a absorção de água obtidos também foram 
baixos em relação aos encontrados na literatura; Gencel et al. (2016), por exemplo, em 
seu projeto de caracterização de um compósito a base de gesso, diatomita e polipropileno, 
encontraram, dentre as amostras ensaiadas, o menor valor de absorção de água de 39,2 
% para um corpo de prova 100 % de gesso. Tal fato indica que os compósitos mensurados 
nesta pesquisa diminuíram a higroscopia, ou capacidade de absorver água, do gesso, o 


que é, majoritariamente, desejado. 


Propriedades Mecânicas 


Os resultados das propriedades mecânicas medidas nesta pesquisa se encontram 


na Tabela 2. 
ENSAIO DE COMPRESSÃO ENSAIO DE FLEXÃO ESTÁTICA 
MOR (MPa) MOE (MPa) MOR (MPa) 
1,48 (+ 0,35) 131,26 (+ 35,73) 2,31 (+ 0,69) 


MOE — Módulo de Elasticidade; 
MOR — Módulo de Ruptura. 


Tabela 2. Resultados médios das propriedades mecânicas encontradas neste estudo seguidos de seus 
respectivos desvios padrões. 


Os dados de MOE e MOR mensurados nos ensaios de flexão e compressão foram 
altos, sendo os valores de MOR, por exemplo, superiores aos dados obtidos por Gencel et 
al. (2016) para o ensaio do corpo de prova contendo apenas gesso, nos quais encontraram 
0,8 MPa para o MOR no teste de flexão estática de três pontos e 1,1 MPa para o MOR no 
ensaio de compressão; isso significa um ganho mecânico das amostras de gesso a partir 
do reforço de capim. 

Entretanto, os corpos de prova do mesmo autor citado acima, com adição de 
diatomita e polipropileno obtiveram maiores resistência, indicando melhores propriedades 


mecânicas que os materiais utilizados nesta análise científica. 
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CONCLUSÃO 


O objetivo desta pesquisa foi a caracterização física e mecânica de compósitos de 
matriz de gesso reforçados com capim para comparação desses resultados com a literatura 
e validação da gama de aplicação de tal composto. 

Perante os resultados obtidos, pode-se concluir que: 


* | À massa específica encontrada dos corpos de prova compósitos reforçados 
com capim foi inferior à de literaturas comparativas, o que pode indicar uma boa 
classificação nos índices de mérito e um aumento de possíveis usos desses 
compostos; 


* Os valores de umidade, inchamento e a absorção de água obtidos também 
foram baixos, o que significa que houve diminuição da higroscopia, ou capaci- 
dade de absorver água do gesso, o que é desejado na maioria das aplicações 
de materiais multifásicos; 


* | Os dados obtidos dos ensaios mecânicos evidenciaram uma melhora nas pro- 
priedades de resistência do gesso a partir da adição de capim, um outro fator 
positivo. 


Portanto, os compósitos de matriz de gesso reforçados com capim representaram 


melhora tanto nas propriedades físicas quanto mecânicas, o que pode validar uma maior 
gama de aplicações para tais materiais. 
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RESUMO: O transporte de cargas é uma 
atividade frequente e que exige certo esforço no 
meio da produção, seja para se ter a otimização 
do layout vigente ou para o alojamento de novos 
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LOCOMOÇÃO DE CARGAS 


maquinários no espaço de trabalho. Para isso 
é requerido aparatos mecânicos que venham 
a facilitar a locomoção desses equipamentos 
de médio e grande porte, ou seja, sistemas 
mecânicos capazes de mover equipamentos 
robustos. Nesse sentido, essa produção tem 
por finalidade desenvolver um sistema de 
transmissão para um guincho, semelhante 
ao de uma ponte rolante, que será instalado 
no laboratório de práticas mecânicas do IFPE 
Campus Caruaru. O projeto aborda desde a 
seleção do motor de indução que opera de modo 
satisfatório para os níveis de projeto, passando 
por etapas de redução por meio de elementos 
de máquinas flexíveis — correia e corrente — e 
pela caixa de engrenagens redutoras. Além 
disso, foi realizado o dimensionamento do tipo 
de cabo de aço e a seleção do material da alma 
do mesmo, esse processo foi realizado seguindo 
a NBR 4309, suas recomendações e cálculos 
necessários. Todas as etapas da metodologia 
são projetadas a fim de que as especificações do 
projeto sejam respeitadas e válidas a níveis de 
segurança como carga média para içamento dos 
maquinários, tempo de elevação e dimensões. O 
dimensionamento desses equipamentos também 
engloba: seleção de material dos componentes, 
tais como nos eixos e nas engrenagens e a 
análise de tempo de vida/ciclos de elementos 
de máquinas com nível de trabalho nominal. 
Engrenagens, de uma forma particular, serão 
projetadas baseadas na metodologia da American 
Gear Manufacturers Association (AGMA) que 
leva em consideração prioritária a resistência à 
flexão e ao desgaste. 

PALAVRAS-CHAVE: Sistema de Transmissão, 
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Projeto, Guincho. 


DESIGN OF MECHANICAL TRASMISSION SYSTEM OF A LIFTING WINCH FOR 
TRANSPORTATION OF EQUIPMENTS 

ABSTRACT: Transportation of loads it's a very commom activity and demands a certain 
effort in the environment of production, may it be for the layout optimization or only for the 
lodge of a new machine in the workspace. For this task, is required tools that make it easy 
the transportation of those equipment's of medium and high size, this means, mechanical 
systems there are capable of move toughness equipment's. In this sense, this production 
aims the development of a transmission system for a winch, like those of an overhead crane, 
that will be installed in Mechanical Practices Lab at the IFPE Caruaru. The project starts 
since the selection of the induction motor, whom work according to the levels of project fixed, 
passing through the reduction phases by means of flexible mechanical elements — belts and 
roller chains — and through the reduction gears. Besides that, it's was made the design of 
the wire rope, the selection of his material, following the NBR 4309 recommendations and 
calculations. Every phase of the methodology was design aiming that the specifications of 
project be respected and valid on the level of safety, like average load for rising of equipment, 
time of rising and dimensions. The design of these mechanical elements also includes 
the selection of material of shafts and gears, and the analysis of the lifetime and cycles 
by operating at normal conditions of power. Gears, in a way, was design according to the 
methodology of American Gear Manufactures Association (AGMA) and considers primarily 
the bending fatigue failure and pitting failure. 

KEYWORDS: Transmission System, Design, Winch. 


11 INTRODUÇÃO 


Os sistemas de transmissão de potência — mais conhecidos por redutores ou 
ampliadores de velocidade — são amplamente utilizados em equipamentos mecânicos 
devido à necessidade de adequação de torque e velocidade fornecidos por um motor, 
em relação ao demandado para a realização de trabalho (BIANCHI, 2009). Sistemas de 
transmissão mecânica têm um grande campo de aplicação no âmbito da mobilidade: como 
setor automotivo, naval e de mecanismos industriais. 

Conjuntos de transmissão são extremamente vantajosos porque são capazes de 
prover à carga um torque tantas vezes maior que as do motor quanto for a relação de 
redução (SANTOS JÚNIOR, 2002). 

Nesse sentido, o correto dimensionamento dos elementos de máquinas para 
solicitação mecânica de um sistema como esse não somente permite um aumento do 
tempo de vida dos componentes, como também assegura ao projetista níveis de segurança 
para o equipamento. De acordo com Niemann (2002), as transmissões por corrente são 
utilizadas para relações de multiplicação de até seis em eixos paralelos com um rendimento 
da ordem de 98%. 

As correias, segundo Melconian (2006), possuem por característica principal 
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a grande capacidade de multiplicação, de absorver impactos e sua perda associada, 
geralmente atrelada ao escorregamento e deslizamento, que representa apenas 5% de 
potência. 

Os sistemas de transmissão mecânicos encontram nas engrenagens seu principal 
elemento de máquina. Pelo fato de permitirem ajustar as características da fonte de rotação 
de acordo com a necessidade do sistema acionado, aumentando o momento, diminuindo 
a velocidade de rotação e vice-versa, com o mínimo de perdas de potência (FARIA, 2009). 
A transmissão por engrenagens apresenta um elevado rendimento em comparação com 
outros sistemas de transmissão mecânica. 


21 APRESENTAÇÃO DO SISTEMA 


O sistema de transmissão desse estudo é composto por um tambor de cabo de 
aço, a caixa de engrenagens redutoras, corrente e correia. Todas as etapas do sistema 
de transmissão estão anotadas com sub-índices numéricos, sendo assim as grandezas 
tais como velocidade e forças associadas a esses estágios estão expressas da forma: 1 — 
motor elétrico, 2 — polia do eixo AB, 3 — entrada do redutor, 4 — saída do redutor. 


c=80mm b = 90 mm a = 90mm 
Ae —————«———————+ 


e=150mm d=150mm 





| 
| 
CABO DE Aço 


! 
| 
| 
| 
| 
| 


CORRENTE 


Figura 1 — Componentes mecânicos do sistema de transmissão 


Fonte: Própria (2019) 


As especificações de projeto para esse sistema foram definidas com base nas 
dimensões do laboratório de práticas mecânicas do IFPE campus Caruaru. Segue, 
portanto as cotas e especificações de projeto: elevar cargas de até 713 kg (7 kN) com uma 
velocidade de 0,3 m/s até uma altura de 6 metros do solo. Motor elétrico de 4 polos (1800 
rom). As cotas a, b, c, de e são distâncias tomadas dos apoios dos mancais na base de 
apoio do sistema de transmissão e são, respectivamente: 90 mm, 90 mm, 80 mm, 150 mm 
e 150 mm. 
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31 PROJETO DOS COMPONENTES 


Partindo-se da velocidade angular do eixo do tambor, tem-se ideia da redução global 
que o sistema necessita, é feita uma ponderação do quanto será o fator de redução em 
cada um dos três elementos (correia, corrente e engrenagens). O organograma ilustra o 
processo de desenvolvimento do projeto. 









DIMENSIONAMENTO 
DO CABO DE AÇO 







ANÁLISE DE FALHA 
POR FLEXÃO E 
DESGASTE (AGMA) 







SELEÇÃO DO MOTOR 
ELÉTRICO 


ANÁLISE DA FADIGA 
NO EIXO DO TAMBOR 


Figura 2 — Organograma do projeto do sistema de transmissão 
Fonte: Própria (2019) 


3.1 Dimensionamento do cabo de aço e tambor 


O estudo do componente cabo de aço consiste em encontrar as cargas de trabalho 
a partir de um fator de segurança recomendado, segundo a NBR 4309, esse coeficiente em 
pontes rolantes varia de 6 a 8. O coeficiente de segurança é a razão da Carga de Ruptura 
Mínima CMR e a Carga de Trabalho CT, como explicita a equação (1): 


Ea (1) 


CT 


A carga de trabalho leva em consideração forças de atrito e inércia, Wé a carga a 
ser elevada em N, m o número de cabos de aço para a operação (m = 1), w é o peso por 
metro do cabo de aço e [é o comprimento suspenso do cabo (altura de elevação) que será 
de 6 metros. Da tabela 17-24 do Shigley (2008), tem-se que para içamento padrão, o cabo 
utilizado é o 6x19 que apresenta um w = 36, 18x 10º N/m,a é referente a aceleração que 
será de 0,015 m/s? e g= 9,81 m/s... 

Pelo fato do diâmetro do cabo ser um fator limitante, visto que não foi selecionado 
previamente, toma-se, portanto, uma média entre os valores do intervalo dado dos 
tamanhos padronizados da tabela 17-24 (Shigley, 2008). Obtendo d = 38mm. Aplicando 
esses valores na equação (2), tem-se: 
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W a 

CT = (— + col) (1 + -) (2) 
m g 
CT =7325N 


Tomando o coeficiente de segurança máximo, 8, se tem: CMR = 58,6 kN 
O manual técnico de Cabos de Aço (CIMAF, 2012) fornece os valores de CMR em tf, 
sendo assim CMR = 5,977 tf de um 6x19 alma de aço, o diâmetro da bitola do cabo (d = 


9,5 mm). Com isso, sabe-se que o diâmetro (D para o tambor se situa numa faixa de 


) 
tambor 
30 a 45 vezes o diâmetro do cabo, pondo uma margem de segurança, o diâmetro do tambor 
=351,5 mm. 


Tendo em mãos o diâmetro do tambor e a velocidade de elevação da carga, 


escolhido será da ordem de 37 vezes o diâmetro do cabo. D 


tambor 


determina-se a rotação do eixo pela expressão (3): 


V=oxr 


Resultando assim, w,=16,3 rad/s. A potência desenvolvida no eixo é tida pela 
expressão (4): 


Ha = (Trambor X CT) x w4 (4) 


Resultando em 2,975 cv, porém, ao considerar perdas nos estágios de transmissão, 
da correia, corrente e engrenagens de 5, 3 e 0%, respectivamente. Tem-se que a potência 
do motor deverá ser de: 


Ha 29755 


Hnotor = (1x0,97x0,95) 0,9215 >” º” 


Por questões de disponibilidade comercial e superdimensionamento de segurança 
para a aplicação, será selecionado um motor de potência 4,0 cv. Com rotação nominal de 
1800 rpm, 220/380 V e 4 pólos. 


3.2 Dimensionamento do eixo do tambor pelo critério de falha 


Para realizar a análise de fadiga, é necessário determinar o torque médio (T je o 
momento alternado (M ). As equações de equilíbrio por meio do diagrama de corpo livre 
(figura 3) fornecem que as reações nos mancais C e D: com as reações de ambos os 
pontos em z sendo nulas, C, = - 3459,645 Ne D, = 3519,645 NM. 
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Figura 3 — Diagrama do corpo livre do eixo CD 
Fonte: Própria (2019) 


O diagrama de Momento Fletor nos fornece o maior momento aplicado no eixo e 
resulta em M = 527,94 N.me T, é obtido pelo produto rx CT, T, =1230,25 N.m. Para um 
eixo rotacionando com flexão e torção constantes a tensão de flexão é completamente 
reversa e a torção estável, logo M = T =0. Pelo critério de falha de DE-Goodman, tomando 


um fator de segurança n = 2. Se terá por (5): 


nt yY 
a = (2/2 [a(o + E [BCs tm) |) 5) 
dep — 53,85 mm 


Alguns fatores e grandezas estão apresentadas na tabela 1: 


Limite de resistência (S,+) p/ E 149,625 
500 MPa Limite de Endurança (S,): 
SAE 1035-HR: MPa 
Fator concentrador de tensão 
Fator de superfície (k,): 0,665 o 1,49 
flexional (K,): 
Fator concentrador de tensão 
Fator de tamanho (k,): 1:13 
torcional (Krç): 


Tabela 1 — Fatores e tensões para análise de fadiga do eixo CD 
Fonte: Própria (2019) 










3.3 Etapas de redução do sistema 


É tomado como fator limitante de projeto que as taxas de redução devam ser menores 


que 10, para não super-requerer aspectos mecânicos dos equipamentos. A figura 4 mostra 
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esquematicamente a redução por estágios relacionado a cada elemento de máquina: 


1800 redução — 5,46 330 fredução = 4,0 825 fredução = 5,0 16.5 
4 — rpm o Wa — rpm Was — 2 rpm 4 — : 2 rpm 
MOTOR ERR | 6.14 DO Eixo As) RD | entrava Do | ND SAÍDA DO 
s REDUTOR di REDUTOR 
Figura 4 — Redução por estágios no sistema 
Fonte: Própria (2019) 





3.4 Dimensionamento e análise de vida das correias 


Na correia, há uma perda de 5% na transmissão da potência nominal do motor, 
assim, a potência desenvolvida na saída do motor 4,0 cv (aproximadamente 3 kW): H m= 
0,95 x H' mz 2,85 kW. 

A tabela 17-12 do Shigley (2008) relaciona as potências para diferentes velocidades 
e diferentes diâmetros mínimos, com base nessa, o diâmetro da polia motora: a, = 155mm. 
Com a redução estabelecida de 5,46, tem-se: a, = 5,46 x 155 = 846,30 mm. 

A velocidade periférica é tida a partir do diâmetro da polia motora, pela expressão 
(6), onde n é a rotação de saída do motor em rps: 


Wer = Tdp,N (6) 


A velocidade periférica de 14,6 m/s situa-se numa faixa de operação boa. 
Interpolando valores de velocidade periférica dado o diâmetro primitivo da roldana e o tipo 
de correia (tipo B) na tabela 17-12 do Shigley (2008), tem-se que a potência tabelada (H,.,) 
é de 2,89 kW, 

Da tabela 17-10, define-se o comprimento da circunferência interna e da tabela 17- 
11 do Shigley (2008), o comprimento a ser adicionado. São, respectivamente: L = 4500 
mme L.,=45 mm. 


Da equação (7), define-se a distância entre centros C: 


T Tm 2 A 
C = 0,25 [Lp — 2 (dp, + dpo)| + [Lp o 2 (dp; F dp,)| = 2(dp, o dp,) (7) 


Resultando C = 1445 mm. Esse dimensionamento respeita a desigualdade d, <C < 
2 
S(d, + d, ), portanto será mantida no projeto. 
2 1 


3.5 Dimensionamento das correntes 


Levando em consideração a capacidade de carga da corrente fazendo e considerando 
a perda de 3% de potência na seção, tem-se que a potência é então de 3,88 cv Uma 
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correção é feita a partir do fator de serviço (k.). Para condições moderadas de operação e 
motor elétrico, tem-se que o fator de serviço é 1,3, sendo assim, por (8): 


Hyroj — kçH 


A potência de projeto é de 5,04 cv A tabela 17-20 do Shigley (2008) nos fornece a 
capacidade em cavalo-vapor para várias velocidades da roda dentada. Foi necessário fazer 
uma interpolação para uma velocidade de 330 rev/min (velocidade da polia do eixo AB). 


BB 4 a | 50 | 60 | 80 | 100 | 120 | 10 | 160 | 180 | 200 | 240 | 





Tabela 2 - Capacidade em cavalos para corrente de filetes de fila única 


Fonte: Shigley (2008). Adaptado pelos autores 


Onde cada cor indica o tipo de lubrificação: tipo A (azul) — por gotejamento ou 
manual; B (verde) — por disco ou por banho; tipo C (amarelo) — por corrente de óleo; tipo 
C' (vermelho) — por corrente de óleo com esfolamento. O fator de correção do número de 
dentes (k) para uma corrente de filetes únicos dentada de 14 dentes é 0,78. Por (9): 


Hcorr = kHproj (9) 


A potência corrigida é de 3,93 cv. Segue abaixo os resultados na tabela 3: 

















aj al Om 6 ) Bo L0Q 120 14 160 180 ano 


| n 
Simples 01557 Qb62))  L5tl?  0,8658 306384) 5,2602] 122928 23,58 39,752 61,475 89154 133,16! Jl63604 253,26 


| Dupla 0,33282 1,1255774 2,67189 197186 5,5208528 8,9490178 H,897% 3999216 6754644 1045153 151,5618 209,3754 2781948 430552) 
Tripla  0,48945 1,65555 292405 21645 76596 1316055 0732 88)  993%3 153,699 222,885 307905 405,11 


Quadru 0,5646074 2,1685356 5,1866] 285714 1011067 1737193 40,56624 7763164 191,1196 2028827 2942082 4054345 5400252 





Tabela 3 - Relação dos tipos de correia para diferentes configurações 
Fonte: Própria (2019) 


Como a capacidade da corrente de nº 60, na configuração simples, foi a mais 
próxima e ligeiramente maior do que a Hd » ela será a aplicada no sistema. 

Ainda pode-se determinar alguns parâmetros da corrente para calcularmos limites 
de potência e fazer sua verificação. O limite de potência baseada na fadiga da placa de 


ligação é dado pela equação (10): 
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Hhrim,pl = 0,004 x Ne x na x p(-0.07P) (10) 


Onde: N, é o número de dentes da roda dentada menor (14), n, é a velocidade 
angular (330 rev/min), e p é o passo (0,75 in). Assim, Himpl = 5,5 cv. Fazendo o uso da 
equação (11), para o limite de potência baseado no rolete e na bucha e da equação (12), 


para o limite de potência baseado no desgaste excessivo. 


1000 x Kr x N,*º x pº8 


mM 1, 


(11) 


lim,RB — 


Nxn xp 
110,84 





Inn, 
) x (4,413 — 2,073 x p — 0,0274 x No) — ( 


Hhrim,Dc = ( 1000 


x(1,59xlogp+1873) (12) 


118 cv e 41,4 cv, respectivamente. Como a potência projetada não ultrapassa esses 
níveis, essa corrente será utilizada no sistema. Com isso, os aspectos geométricos e 
dimensionais da corrente seguem apresentados na tabela 5. 


Diâmetro do pinhão (D,): 


Diâmetro da coroa (D..): 13,37 in 
Número de passos (L): 





Tabela 5 - Especificações da corrente 
Fonte: Própria (2019) 


3.6 Dimensionamento das engrenagens redutoras 


Para o dimensionamento do par de engrenagens, tem-se que o fator de redução a 
título de projeto é 5, logo, por (13). Os sub-índices 3 e 4 são relativos ao pinhão e coroa, 
respectivamente. 


r,=5X7T3 (13) 


Iniciando com uma análise de interferência para engrenagens cilíndricas de dentes 


retos, tem-se que o número mínimo de dentes no pinhão é dado por (14): 


— 2k | 2 psiê 
Ns = ramo sa (mo + mg + (1 + 2mç)sin o) (14) 


Para uma razão de engrenamento igual a 5 e tomando o ângulo de pressão como 
20º e que k =1, o menor número de dentes do pinhão para que não ocorra interferência é 


16. Para a coroa, o maior número de dentes da engrenagem dada pela equação (15): 
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— Nósin?6-4 (15) 
*" 4-2Nssin? 6 

A expressão nos dá como resultado 100 dentes, logo, dado o intervalo obtido, foi 

projetado um pinhão com 18 dentes, e uma coroa com 90. Tomando como fator de projeto 


um módulo (m) de 6mm, tem-se, através dos devidos cálculos: 





Númerodedentes | 8 | 00 
0,2364 / 6,0 mm 
0.2955 / 7.5 mm 

13,5 

120 552 

528 


Tabela 6 — Dimensões do par engrenado 
Fonte: Própria (2019) 


















3.7 Análise de falha de flexão e desgaste das engrenagens 


Uma vez projetado os aspectos dimensionais, a análise das forças é imprescindível 
para o projeto contra falha de desgaste e flexão. Calculando a carga transmitida no pinhão 
por (16): 

my, = Pi, = SSO0ON 19 

Sendo H a potência que chega no redutor, passando pelos estágios de redução de 
correia e corrente, tem-se que: H, = H tor X Perdas = 4 x 0,97 x 0,92 = 3,97 cv = 2625,731 
W 

Com d, = 108 mme a celeridade n,= w,=82,5 rev/min. Portanto: a carga transmitida 
é de 5,628 kN Sabendo que o ângulo de contato q é 20º, pode-se determinar a componente 
radial e o módulo da força. A figura 5 abaixo mostra um diagrama de corpo livre das forças 


aplicadas no par de engrenamento. 


Fi, = 2,048 kN: Fs, = 5,989 kN 
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Figura 5 — Diagrama de corpo livre das forças aplicadas na engrenagem 


Fonte: Própria (2019) 


Realizando a análise da velocidade na linha primitiva, e sabendo que d é dado em 


polegadas e n em rpm, temos por (17): 


— mdans (17) 


V=—>º2 
12 


A análise de falha de desgaste e flexão da metodologia da AGMA é feita em unidades 
inglesas, portanto, a carga transmitida será W=1265,25 lbf As equações de tensão da 


AGMA, considerando todos os fatores, são duas a flexional (18): 


Pi KnK 
o; = WEKoKoKs (18) 


E a de resistência ao crateramento (19): 


| Km C 
0. =Cy W'KoKyKs E (19) 
p 


As equações do coeficiente de segurança AGMA S, para flexão e desgaste S, são 


expostas a seguir por (20) e (21), respectivamente: 








É Of adm KrKr ( 
Sc Zulu (21) 








Ro Ocadm KrKn 
O pinhão é de aço grau 1 com dureza de superfície e do dente de 240 HB (módulo 


de elasticidade 30 Mpsi), a coroa é de aço grau 1 e possui dureza de superfície e do dente 
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de 200 HB (módulo de elasticidade 2300 Mpsi). O pinhão irá rotacionar por 10º ciclos e 
a coroa por 2x10” ciclos. Largura da face do par engrenado 2,97in, possuem dentes não 
coroados e por consideração de projeto os mancais estão imediatamente adjacentes ao 
redutor e as engrenagens do redutor estão confinadas. 

O redutor também é comercialmente disponível, por fim, toma-se uma confiabilidade 


de 99%. Aplicando as devidas relações, tem-se todos os resultados tabelados abaixo (7): 


[esa [iss [Sea [rszows [1 fora 


Po [sescomesm [cha [100 [ka [100 + 


Tabela 7 - Fatores e coeficientes do projeto 
Fonte: Própria (2019) 





Calculando a falha por flexão no pinhão e na coroa por (18): 
(Of adm), = 8157,82 psi; (Of adm), = 6267 psi 
E consequentemente o coeficiente de segurança por (20): 
(Sr)3 = 3,75; (Sp)s = 4,53 
Calculando a falha por desgaste no pinhão e na coroa por (19): 
(Scadm), = 75724 psi; (Ocadm), = 76037 psi 
E, da mesma forma, obtendo os coeficientes de segurança por (21): 


(Sg)3 — 1,33; (Sg)4 — 1,178 


41 RESULTADOS E DISCUSSÃO 


Realizar esse trabalho permite não somente aplicar os conceitos de seleção de 


materiais e dimensionamentos de elementos de máquinas como também integrar nossa 
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visão no que diz respeito ao funcionamento de conjuntos mecânicos. Tomar a liberdade de 
projetar, fazer considerações a respeito de aspectos construtivos e aplicar metodologias de 
grande importância no âmbito de projetos mecânicos. 

É um trabalho multidisciplinar no âmbito da engenharia mecânica e nos permite 
deixar uma contribuição para o nosso campus, não só na pesquisa e produção acadêmica, 
como também para o laboratório de práticas mecânicas que possui um norte para a futura 
instalação do sistema de içamento. 
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RESUMO: Com constantes mudanças em 
nossos hábitos, surgem também novas 
necessidades. A partir daí surgiu a ideia de se 
criar uma solução para resfriar e manter bebidas 
frias através de um equipamento portátil e 
eficiente. Propondo-se a desenvolver, criar e 
validar um protótipo que atenda essa ideia este 
artigo busca apresentar informações e dados 
sobre o dimensionamento do equipamento, 
seleção de componentes, fabricação, validação 
e a apresentação dos resultados obtidos com o 
protótipo. 

PALAVRAS-CHAVE: Refrigerador Portátil, 
Peltier, Transferência de Calor, Prototipagem 3D. 
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DESIGN AND CONSTRUCTION OF A 
PORTABLE REFRIGERATOR BASED ON 
THE PELTIER EFFECT 


ABSTRACT: With constant changes in our habiis, 
new needs also appear. From there came the 
idea of creating a solution to cool and keep cold 
drinks using a portable and efficient equipment. 
To develop, create and validate a prototype 
that meeis this idea, this article aims to present 
information and data on equipment sizing, 
component selection, fabrication, validation, and 
the results obtained with the prototype. 
ABSTRACT: Portable Cooler, Peltier, Heat 
Transfer, 3D Prototyping. 


11 INTRODUÇÃO 


Atualmente existem no mercado 
vários tipos e modelos de refrigeradores, com 
grandes variações quanto às dimensões e 
potências, visando atender a todos os tipos de 
consumidores e suas necessidades em uso 
Podendo 


respectivamente manter alimentos e bebidas 


doméstico, comercial e industrial. 


resfriados e aumentar a conservação ideal de 
alguns alimentos. A capacidade de temperatura 
dos refrigeradores varia na faixa de -8ºC a 
-18ºC para o compartimento de congeladores 
e 2ºC a 18ºC no compartimento de resfriados 
(SILVA, 2019). bebidas 
esquentam quando expostas à temperatura 


Naturalmente, 


ambiente, sendo que as formas tradicionais de 
resfriamento de bebidas incluem uso de gelo ou 
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refrigeradores; métodos esses, que levam tempo para um consumo com paladar agradável 
e exigem grandes compressores mecânicos, caixas térmicas e um grande espaço físico. 
Devido às mudanças sociais surgiu a necessidade de desenvolver um equipamento de 
fácil acesso, portátil, compacto, com capacidade de refrigerar e conservar os líquidos em 
poucos minutos, ocupando um pequeno espaço físico, sem a utilização de gelo e sem 
necessidade de um compressor, o que diminui muito a emissão de ruídos. 

Para isso, é necessário o controle da climatização através de módulos de 
resfriamento termoelétrico, utilizando dispositivos eletrônicos como princípio de refrigeração 
(COSTA, 1982), obtendo vantagem sobre os métodos convencionais de refrigeração, pois 
são simples e não tem partes móveis a nível macroscópico. Além disso, são compactos, 
confiáveis, silenciosos e também não utilizam gás refrigerante que prejudicam a camada 
de ozônio (MORAN et al., 2013). 

Estes dispositivos são módulos (células ou pastilhas), baseados no efeito Peltier e 
são capazes de bombear calor para produzir um resfriamento ou aquecimento (MORAN 
et al., 2013). Importante destacar, que no núcleo de um refrigerador termoelétrico há 
dois materiais diferentes, neste caso semicondutores do tipo-n e do tipo-p. Para serem 
eficientes na refrigeração termoelétrica, os materiais devem ter baixa condutividade 
térmica e elevada condutividade elétrica (MORAN et al., 2013). O dissipador de calor é 
uma peça fundamental para o bom funcionamento de diferentes tipos de equipamentos e 
maquinários eletrônicos, pois possuem a função de facilitar a troca térmica do sistema com 
o meio ambiente, sendo a seleção ou desenvolvimento do dissipador de calor essencial na 
operação das pastilhas termoelétricas (CENGEL; GHAJAR, 2012), já que os módulos não 
retiram calor, e sim, transportam de uma face para outra. Por fim o objetivo desse trabalho 
é projetar e construir um refrigerador portátil para líquidos, além de validar sua eficiência 


após finalizado o protótipo. 


2 | MATERIAIS E MÉTODOS 


Para a fabricação e testes de funcionamento do equipamento que terá a capacidade 
de refrigerar e conservar líquidos, foi necessário o uso de pastilhas termoelétricas. 
Para a escolha do melhor módulo termoelétrico além do dimensionamento, foi efetuado 
testes práticos com os módulos TEC1-12706 e TEC1-12708, alimentados por uma fonte 
convencional de 12 V. 

Para o sistema de dissipação térmica, foram feitos testes práticos com quatro 
dissipadores de áreas de contato variando entre 54800 mm? a 295600 mm, juntamente com 
dois modelos de coolers de diferentes potências. Após definido o conjunto de dissipação, foi 
desenvolvido um modelo 3D para fabricação da carcaça e simulações aerodinâmicas, com 
o objetivo de melhorar o fluxo de ar e identificar pontos de melhoria na carcaça, produzida 
por meio de impressão 3D, utilizando filamento em PLA. 
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O sistema para a condução de calor da face fria da pastilha termoelétrica para o 
liquido apelidado de conjunto de condução (copo e base), será produzido por um prestador 
de serviços usinando um tarugo de alumínio de liga 6351-T6, utilizando um torno mecânico. 

Já com o protótipo finalizado, todo o conjunto será submetido a testes com uma 
câmera térmica FLIR embutida no Smartphone Caterpillar S60, e também comparações 
com o sistema de testes e com uma geladeira convencional modelo DFX41 da Electrolux. 


31 RESULTADOS E DISCUSSÃO 


Para estimar a carga térmica do sistema, é necessário fixar alguns parâmetros 
iniciais com o objetivo de realizar os cálculos para o dimensionamento do sistema. Como 
referência, foi considerado liquido (H,0) como objeto a ser refrigerado com volume 
de 220 ml, conhecendo a densidade da HO a 25º€, 0,9970 g/ml (LIDE, 2004) tem-se 
aproximadamente 220 g de massa; considerando a temperatura inicial do líquido e do 
ambiente a 25ºC e a temperatura final pretendida para resfriamento do líquido de 5ºC, 
em um tempo aproximado de 30 min. Sendo assim, é possível estimar a carga térmica 
utilizando a Equação (1) (MILHER; MILHER, 2014). 

Q=m-c:AT (1) 

Onde: Q - quantidade de calor sensível (cal ou J), m - massa (9), c - calor específico, 
constante em (cal/ºC) e AT - variação de temperatura (ºC). 

Portanto: Q =220-1-(25-5) = 4340 cal ou 18410 J 

A quantidade de calor que deve ser retirada do líquido para resfriamento até a 
temperatura desejada, é de 18410 J. Pode-se chegar à potência aproximada do equipamento 
através da Equação (2). 

P=QA (2) 

Onde: Q - Quantidade de calor sensível (cal ou J), P - potência (W) e t - tempo de 
operação (s). 

Portanto: P = 18410/1800 = 10,2 W 

A partir das condições operacionais estabelecidas para o projeto, chegou-se a uma 
potência mínima de 10,2 W, o que permite agora escolher a pastilha mais adequada. 

Com base na curva de performance fornecido pelo fabricante, foi cogitado para 
este projeto as pastilhas termoelétricas TEC1-12706 e TEC1-12708 (Figura 1), levado em 


consideração o custo, sua potência, suas características disponibilizadas pelo fabricante. 
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Figura 1. Ilustração pastilha termoelétrica (PELTIERMODULES, 2020) 


Analisando os dados fornecidos pelo fabricante (PELTIERMODULES, 2020), tem- 
se a corrente, tensão e potência máxima de operação do módulo. Para que a pastilha 
mantenha esses valores em operação, é necessário que se mantenha uma tensão e 
corrente constante. 

Com a principal fonte de consumo elétrico do projeto definido, levando em 
consideração outros elementos que irão ser agregados ao sistema, foi adquirido uma fonte 
de 12 V de 10 A, sendo uma fonte compacta quando comparada a outros modelos, e com 
potência similar, além de ser um modelo muito simples de ser usado. 

Para atender as necessidades de designer do protótipo, foi feito um projeto 3D para 
a criação do protótipo, que foi dividido em 2 partes (parte 1 - conjunto copo e base, e parte 
2 - conjunto de isolamento térmico). 

Conjunto de condução (Figura 2a), sendo constituído respectivamente pelo copo e 
base (Figuras 2b e 2c), é responsável pela troca de calor entre o líquido que se pretende 
resfriar e a pastilha termoelétrica. Para este item, foi utilizado um tarugo de alumínio que 
foi usinado em um torno. 

Parte 2 é responsável pelo isolamento térmico entre o conjunto de condução, o 
ambiente e o dissipador, fornecendo a sustentação para a fixação dos componentes (Figura 
2d). O item foi produzido utilizando uma impressora 3D (Ender 3) com parâmetros definidos 
em 60 mm/s, temperatura do extrusor em 205ºC, largura da camada é de 0,40 mm e 0,20 


mm de altura, utilizando o filamento cinza com diâmetro de 1,75 mm em PLA. 
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Figura 2. (a) - conjunto de condução em alumínio, (b) - copo em alumínio, (c) - base em 
alumínio, (d) - carcaça e (e) - conjunto montado, dimensões em mm 


Na escolha do dissipador foram feitas uma bateria de testes utilizando 4 modelos 
de dissipador com área de contato diferentes (Tabela 1). Todos foram montados em um 
conjunto com dois modelos de cooler. Sendo o primeiro cooler de modelo A9225-42RB- 
4AP-C1 da fabricante Cooler Master que possui 7,20 W, 64 CFM e ruído de 61 dB, e o 
segundo cooler de modelo RDL8025S da fabricante Ruilian Science, possuindo 0,98 W, 26 
CFM e ruído de 38 dB. 


Dissipador Área de contato aproximada (mm?) 
Aí 54800 
A2 127600 
A3 295600 
A4 173000 


Tabela 1. Dissipador com área de contato diferentes 


Os testes foram realizados utilizando a pastilha termoelétrica TEC1-12708 com 220 
mi de H,O na temperatura ambiente de 25ºC durante 30 min, monitoradas a cada 5 min. 
O termômetro utilizado possui uma escala que vária em 1ºC, e a cada troca do dissipador, 
o conjunto de condução era submergido em um balde com água durante 5 min para que 
ficasse em temperatura ambiente antes da próxima bateria de teste. 

Analisando a Figura (3a), o conjunto com cooler A9225-42RB-4AP-C1 e dissipador 
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A3, demostram melhor eficiência na dissipação térmica devido haver uma maior área de 
contato para troca térmica por convecção forçada. Comparando os demais conjuntos com 
dissipador AS, o dissipador A4 apresenta uma perda de eficiência aproximada de 12%, para 
o dissipador A2 a perda é de aproximadamente 24,50%, e para o dissipador Aí a perde 


chega a aproximadamente 29,50%. 


TEMPERATURA X TEMPO (a) TEMPERATURA X TEMPO (b) 
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Figura 3. (a) - cooler Cooler Master e (b) - cooler Ruilian Science 


Analisando a Figura (3b), o conjunto com cooler modelo RDL8025S e dissipador AS, 
demostra melhor eficiência na dissipação térmica devido a uma maior área de contato para 
troca térmica por convecção forçada. Comparando os demais conjuntos com dissipador 
A3, o dissipador A4 apresenta uma perda de eficiência aproximada de 12,50%, para o 
dissipador A2 a perda é de aproximadamente 25%, e para o dissipador A1 a perda chega 
a aproximadamente 37,50%. 

Após comparar os resultados (Figuras 3a e 3b), é possível notar que o dissipador 
A3, apresenta uma melhor eficiência independente do cooler utilizado. Seu tamanho 
inviabiliza o uso para este projeto, por não ser tão compacto e ter 700 g. Já o dissipador 
A1, tem um tamanho compacto excelente para o projeto, mas seu desempenho fica a 
desejar, por possuir uma menor área de contato para troca térmica por convecção forçada. 
Sendo assim, do ponto de vista técnico e comercial e a partir da análise das (Figuras 3a 
e 3b), nota-se que o conjunto com o cooler de modelo A9225-42RB-4AP-C1 e dissipador 
A4 (Figura 4a), é o que atende melhor às perspectivas do projeto refrigerador de líquido. 

Durante os testes para a escolha do dissipador e do cooler, foi notado uma leve 
condensação na parte central da base demonstrado na Figura (4b), a região compreendida 
dentro do círculo vermelho, sendo uma formação quase totalmente homogênea diferente 
das formações nas laterais indicadas pelas setas pretas. Após analisar, ficou claro que essa 
falha é proveniente da usinagem da base (Figura 4c). Por consequência desse mal contato, 
perde-se uma grande área para a troca de calor impactando diretamente no desempenho 


do equipamento. 
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Figura 4. (a) - conjunto com dissipador A4 e cooler A9225-42RB-4AP-C1, (b) - falha na 
usinagem e (c) - demonstração da falha (folga) na fabricação do conjunto condução 


Afim de quantificar essa falha, foram feitos novos testes adicionais utilizando a 
pastilha termoelétrica TEC1-12708 e cooler modelo A9225-42RB-4AP-C1 também com 
220 ml de H,O na temperatura ambiente de 25ºC durante 30 min monitoradas a cada 5 
min. Entretanto, dessa vez, diferente do teste feito anteriormente (A4 + Copo e Base), foi 
utilizado somente o Copo (A4 + Copo) em contato direto com a pastilha termoelétrica e os 
resultados podem ser vistos na Figura (5a). A diferença final é de 3ºC, perdendo nesse caso 
20% do seu desempenho uma quantidade muito relevante diante da falha apresentada. 

Logo em seguida foi feito um teste com a finalidade de quantificar a diferença no 
desempenho quando utilizado a pastilha TEC1-12706 comparando ao teste anterior (A4 + 
Copo). Os dados obtidos podem ser vistos na Figura (Db), é possível notar uma diferença 
final de 3ºC, equivalente a 20% menos desempenho. No entanto o seu consumo durante 
o teste, ficou em 5,40 A (64,80 W) contra 6,96 A (83,50 W) 28% a menos de consumo. A 
hipótese para essa diferença, é que com uma menor emissão de calor para ser dissipado, 
o dissipador e o cooler conseguiram ter um ganho na sua eficiência de 8%, uma vez que, 
o funcionamento da TEC1-12708 e TEC1-12706 é o mesmo. Sendo assim, optou-se por 
manter a utilização da TEC1-12708 no protótipo. 
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Figura 5. (a) - impacto da falha na fabricação do conjunto condução e (b) - diferença de 
desempenho entre as pastilhas 
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Após a definição do sistema de dissipação térmica e da pastilha termoelétrica, foi 
efetuado uma simulação aerodinâmica com objetivo de observar o impacto do fluxo de ar 
na dissipação de calor, conforme apresentado na Figuras (6a). Para a simulação, utilizou- 
se parâmetros do próprio software SolidWorks, onde é possível simular a interação das 
partículas de ar com a hélice do cooler dentro de uma zona delimitada, sendo também 
possível, introduzir a rotação do sistema, que para esse caso, foi fixado em 4200 rpm, 
gerando assim, interações bem próximas da realidade. 





Figura 6. (a) - direcionamento do fluxo de ar vista isométrica.e (b) - refrigerador portátil 


Finalizada a impressão da carcaça, inicia-se a montagem do protótipo. Foi incluído 
o circuito elétrico, que contempla a fiação, conector para energia, conector para o cooler e 
chave geral. Já na fixação da pastilha termoelétrica, foi utilizado pasta térmica nas áreas 
de contato com o dissipador e com o conjunto de condução facilitando o fluxo de calor. 
O dissipador foi fixado a carcaça por meio de parafusos, assim como a base do copo. O 
resultado final está apresentado na Figura (6b). 

Foi feito um teste com o protótipo finalizado sendo possível observar uma discreta 
melhora no seu desempenho quando comparado ao conjunto inicial A4, demostrando 
assim, que a carcaça cumpriu seu papel no isolamento do sistema (Figura 7a). 

Utilizando a câmera térmica fixada ao Smartphone Caterpillar S60, com uma 
precisão de 0,10ºC, uma resolução de imagem de 640 x 480 pixels e com uma temperatura 
ambiente de 21ºC, foi possível acompanhar a dissipação de calor no protótipo (Figura 18), 
sendo monitorada durante 30 min. 
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Figura 7. (a) - gráfico comparativo A4 x Protótipo final e (b) - imagem termográfica após 30 min 
em funcionamento 


Analisando as imagens, foi possível validar a simulação aerodinâmica. Comparando 
as imagens observadas pela câmera térmica, com a situação esperada (Figura 6a) e 
destacando a área compreendida dentro da marcação com linhas tracejadas (Figura7b), 
nota-se uma coloração mais avermelhada, sendo este o resultado da saída do ar quente da 
carcaça, que mesmo não incidindo diretamente, influencia na troca de calor do isolamento 
do copo com o ambiente, devido à sua proximidade. Uma situação prevista na simulação e 
marcada previamente como ponto de atenção. 

Neste teste (Figura 7b), o protótipo apresentou uma temperatura máxima de 28,1ºC 
em sua chave (devido a passagem da corrente elétrica), uma mínima de 8,1ºC na parte 
externa do copo (o liquido chegou a uma temperatura de 5ºC), sendo que na maior parte 
da carcaça a temperatura ficou por volta de 24,1ºC. 

Também foi possível notar o impacto da temperatura ambiente no desempenho do 
protótipo, atuando de forma constante em seu resultado, o que facilita aproximações para 
outras condições iniciais próximas as contempladas nos testes. 

Refletindo sobre todo o projeto, foi possível identificar pontos de melhorias que 
impactam diretamente no seu desempenho, como a correção da falha no conjunto de 
condução, a incidência de ar quente no isolamento do copo e uma possível tampa para o 
copo, uma vez que, com esses detalhes alterados, estima-se ser possível atingir em uma 
temperatura de 25ºC para 3ºC em 30 min. 

No entanto, levando em consideração todas essas mudanças e retomando os 
cálculos iniciais utilizando a Equação (1) para a carga térmica, tem-se: Q = 220-1-(25- 
3) = 4840 cal ou 20260 J. Aplicando a energia na Equação (2) da potência, tem-se: P = 
20260/1800 = 11,25 W. 

Levando em consideração que o sistema consome 83,50 W o rendimento do 
sistema como um todo, mesmo fazendo as correções possíveis, fica em torno de 14% uma 
eficiência energética baixa que em sua maior parte se dá ao baixo rendimento da pastilha 
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termoelétrica. 

No entanto, quando comparamos o desempenho do protótipo a uma Geladeira Frost 
Free modelo DFX41, ainda assim, tem-se um excelente desempenho. Para esse teste foi 
utilizado 220 ml de água a 25ºC e um copo padrão de vidro com espessura de sua parede 
lateral de 2,40 mm e com seu fundo de 5,25 mm em vidro com um volume de 225 ml 
monitorado a cada 5 min com duração total de 30 min. Os resultados podem ser vistos na 
Figura (8a). 
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Figura 8. (a) - gráfico comparativo com sistemas convencionais de refrigeraçãoe (b) - gráfico 
comparativo entre os equipamentos 


Importante destacar, que no momento dos testes a geladeira estava em 
funcionamento, soprando uma rajada de ar contínua no seu interior, promovendo uma 
melhor troca de calor. 

Em sua parte mais quente (Refrigerador), que no momento do teste a 2ºC, obteve 
uma temperatura de 16ºC. Já na sua parte mais fria (Freezer), que no momento do teste 
estava em -23ºC, obteve uma temperatura de 5ºC. O protótipo final partindo das mesmas 
condições conseguiu atingir 9ºC, entretanto após fazer as correções descritas anteriormente 
o protótipo poderá atingir 3ºC ou até menos nesse mesmo tempo. 

Adotando o protótipo como 100% de desempenho, tem-se que o refrigerador 
perdeu por uma diferença de 44% e o Freezer conseguiu ganhar por 25%. Já o protótipo 
quando forem aplicadas as melhorias descritas anteriormente, chegará a 38% a mais de 
desempenho, apresentando um resultado excelente, tendo em vista, suas proporções e 
limitações (Figura 8b). 


41 CONCLUSÃO 


Buscando criar um equipamento inovador, este trabalho atingiu com êxito seu 
objetivo de projetar, construir e testar um protótipo para o resfriamento de líquidos, 
utilizando o efeito Peltier; sendo possível baixar a temperatura de 220 ml de água a 25ºC 
para 9ºC em 30 min, além de colocar à prova as ferramentas, métodos e equipamentos 
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usados nesse projeto. 

Outro destaque importante, foi a comparação do protótipo final a um modelo de 
refrigerador convencional, e mesmo frente a este desafio, foi possível obter um excelente 
resultado, chegando a 44% a mais de desempenho que o Refrigerador, demonstrando 
que o protótipo, mesmo antes de aplicar correções e melhorias, pode sim, se tornar um 
excelente produto, proporcionando uma solução pequena, portátil, segura e eficiente para 
refrigeração. 
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RESUMO: O desenvolvimento de nanofluidos 
é extremamente importante para aprimorar a 
condutividade térmica dos fluidos base. Uma 
das maneiras de se melhorar a segurança das 
centrais nucleares está relacionada à melhoria 
da sua capacidade de transferência de calor. O 
estudo de novos fluidos que melhorem a taxa 
de remoção de calor é fundamental para se 
obter maior eficiência de sistemas energéticos. 
Dentre os diversos fatores que comprometem a 
eficiência dos sistemas energéticos, podem-se 
destacar as limitações termofísicas dos fluidos 
convencionais. Estas limitações inibem, de 
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forma bastante significativa, algumas aplicações 
industriais. Este artigo apresenta uma proposta 
de trabalho a ser realizada no Laboratório de 
Termo-hidráulica do CDTN, que tem como 
objetivo o estudo da melhoria das características 
de transferência de calor do refrigerante, utilizado 
no circuito primário dos reatores nucleares, 
através da adição de nanopartículas. Serão 
utilizados, principalmente, nanotubos de carbono 
em água leve, que é o fluido mais utilizado para 
refrigeração dos reatores nucleares. Como 
suporte ao trabalho será projetada, montada e 
calibrada uma bancada experimental, com intuito 
de estudar o comportamento termo-hidráulico 
destes componentes. 

PALAVRAS-CHAVE:  Nanofluidos; Eficiência 
Térmica; Termo-Hidráulica; Reatores Nucleares. 


11 INTRODUÇÃO 


Nanofluidos são conhecidos | por 
apresentar possibilidade de transferência de 
calor superior e são, portanto, investigados 
E dj O) 


estudo de novos fluidos que melhorem a taxa 


para aplicações em engenharia [1, 


de remoção de calor, é fundamental para se 
obter maior eficiência de sistemas energéticos. 
Dentre os diversos fatores que comprometem 
a eficiência dos sistemas energéticos, pode- 
se destacar as limitações termofísicas dos 


fluidos convencionais, inibindo de forma 
bastante significativa algumas aplicações 
industriais. Neste trabalho, pretende-se 


estudar as características de transferência de 
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calor de fluidos comumente utilizados, através da adição de nanopartículas, constituídas 
de nanotubos de carbono, em água leve. Água é o fluido mais utilizado, atualmente, na 
refrigeração de reatores nucleares. 

A melhoria das propriedades térmicas é limitada pela condutividade térmica do 
fluido de trabalho. Com objetivo de superar tais limitações, diversos estudos de fluidos com 
partículas sólidas em suspensões vêm sendo realizados, desde a primeira publicação de 
Maxwell no final do século XIX [4]. Segundo [5], entretanto, tais estudos estão limitados à 
utilização de partículas com escala micro, que apresentam algumas desvantagens como: 


* | Sedimentação — as partículas, após algum tempo, sedimentam-se, criando uma 
camada junto às superfícies reduzindo a capacidade térmica. 


*  Desgaste — é possível reduzir a sedimentação, através do aumento da veloci- 
dade de circulação do fluido, mas com aumento do desgaste dos equipamentos 
e tubulações, etc. 


*  Obstrução — devido às dimensões das partículas os canais tendem a ficar obs- 
truídos, particularmente nos estreitos canais de refrigeração. 


* | Perdas de carga — aumentam de forma bastante significativa. 


* | Condutividade térmica — é proporcional à concentração de partículas, mas po- 
tenciando problemas anteriores. 


21 NANOFLUIDOS 


Nanofluidos são constituídos por dispersões coloidais com um refrigerante 
tradicional como uma base, na qual as nanopartículas são suspensas. Em 1995, Choi [6], 
foi o primeiro a utilizar o termo “nanofluidos”, e publicou os resultados de sua pesquisa 
teórica. Desenvolvimentos subsequentes em engenharia de nanofluidos contribuíram para 
o rápido crescimento da nanotecnologia e tecnologias de superfície ao longo da última 
década [7]. 

As suspensões coloidais mostraram substancialmente intrigantes do ponto de vista 
dos desempenhos térmicos, sobre quatro pontos [6]: 


1. Aumento da condutividade térmica (aprox. 150%); 
2. Aumento monofásico do coeficiente de transferência de calor (aprox. 60%); 


3. Aumento do fluxo crítico de calor com regime de ebulição nucleada estendida 
(aprox. 200%); 


4. Aumento da eficiência de resfriamento. 


Verificou-se que apresentam propriedades térmicas melhoradas, por exemplo, 
quando atuando como líquido de arrefecimento, elevam o ponto de fluxo crítico de calor e 
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molhabilidade da superfície em pequenas concentrações, o que é uma característica útil 
nos reatores nucleares. 

Em trabalho recente [8], avaliaram experimentalmente nanofluidos de nanotubos de 
carbono de parede múltipla (MWCNTs- Multi-walled Carbon Nanotubes) a base de água, 
com concentrações em volume de O < P < 0,24%, para vazões mássicas de 20 g/s até 90 
g/s. O calor fornecido pela parede do tubo através da resistência na seção de teste variou 
entre 500 W até 900 W. A avaliação do desempenho termo-hidráulico das três amostras 
de nanofluidos testadas foi realizada especificando o coeficiente relativo de transferência 
de calor por convecção, em função da potência de bombeamento relativa obtida nas 
potências de bombeamento na seção de teste. As três amostras estudadas, apresentaram 
um incremento na condutividade térmica, em média, de 6,3%, 8,8% e 17% superiores 
a do fluido base (água destilada). Demostrando o efeito da concentração no aumento 
desta propriedade. No entanto as três amostras também apresentaram um incremento na 
viscosidade em média de 8,8%, 17,8% e 11,6%. 

Muito se tem estudado a respeito do desempenho da transferência de calor em 
diversos tipos de escoamento utilizando nanofluidos, constituídos de nanopartículas de 
diferentes materiais e diversos fluidos base, conforme mostrado na Tabela 1. 


, . ) Incremento do coeficiente 


RS Turbulento ALO, em água +45 


; -7 para 3.5vol% 


, +26 para 1vol% 


+15%, mesmo número de 
Reynolds. 

-10%, mesma potência de 
bombeamento. 


[14] Turbulento TiO, e ALO, em água 





Tab. 1. Comparação dos trabalhos realizados que analisaram a transferência de calor por 
convecção em nanofluidos. 


Fonte: [8] 


Ainda que muitos dos estudos relatados mostrem que houve aumento da transferência 
de calor por convecção forçada utilizando nanofluidos, escoando em regime laminar e 
turbulento, outras pesquisas mostram inconsistências nos resultados apresentados, ou 
seja, sob certas condições relatam incremento no coeficiente de transferência de calor. 
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Há situações onde não é verificado acréscimo tampouco diminuição no coeficiente de 
transferência de calor por convecção quando nanopartículas são adicionadas aos fluidos 
base. Contudo, nota-se em alguns casos a redução do valor do coeficiente de transferência 
de calor por convecção. Portanto, ainda são necessários mais estudos que validem e 
verifiguem o desempenho termo-hidráulico destes novos fluidos de arrefecimento [8]. 

Foi realizado um estudo sobre a utilização de nanofluidos para sistemas de 
resfriamento de emergência do núcleo, que é uma característica de segurança das centrais 
nucleares [3]. Tais estudos mostraram que houve um aumento do fluxo de calor crítico, 
através da injeção de nanofluido em caso de perda da refrigeração por acidente. 

Foram obtidas experimentalmente as curvas de inibição (temperatura versus tempo) 
em pequenas esferas (-1 cm) metálicas expostas a água pura e nanofluidos à base de 
àgua com alumina, sílica e nanopartículas de diamante em baixas concentrações (<= 0,1% 
Vol) [16]. As esferas utilizadas eram feitas de aço inoxidável e zircaloy, e foram resfriadas 
a partir de uma temperatura inicial de, aproximadamente, 1000ºC. Os resultados mostram 
que o comportamento de têmpera em nanoíluidos é quase idêntica à da água pura. No 
entanto, verificou-se que algumas nanopartículas acumular na superfície da esfera, o que 
resulta na desestabilização da película de vapor em testes subsequentes com a mesma 
esfera, acelerando, assim, muito o processo de resfriamento. 

Utilizou-se nanofluidos constituídos de dispersões diluídas de alumina, Zzircônia, e 
nanopartículas de sílica em água. Vários parâmetros que afetam a transferência de calor 
(ponto de ebulição, viscosidade, condutividade térmica, e a tensão superficial de ebulição) 
foram medidos e de forma consistente com outros estudos nanofluido [17], verificou-se 
ser semelhante aos de água pura. Entretanto para pool boiling experimentos mostraram 
melhorias significativas do fluxo crítico de calor nos nanofluidos (até 200%). 


31 METODOLOGIA 


O principal objetivo desta pesquisa é o estudo e melhoria da condutividade térmica 
em sistemas de nanofluidos com adição de nanotubos de carbono. Para isto existem vários 
métodos, que se pretende utilizar, para a medição da condutividade térmica, ou sejam: 
oscilações de temperatura, placas paralelas em regime permanente (steady-state parallel 
plate) e a técnica do fio quente transiente (transient hot wire), sendo esta última a mais 
utilizada [18]. 

A técnica do fio quente funciona medindo a resposta temperatura/tempo do foi a um 
pulso elétrico abrupto. Nesta técnica o fio funciona como aquecedor e termômetro, sendo a 
condutividade térmica calculada através de uma derivação da Lei de Fourier. 

Sistemas com suspensão de nanotubos de carbono têm apresentados os valores 
mais elevados para a condutividade térmica. Um estudo realizado por [19] sobre nanofluidos 
constituídos de nanotubos de parede múltipla adicionados à óleo, com objetivo de medir a 
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condutividade térmica efetiva. Os resultados experimentais cnamaram atenção por serem 
muito superiores aos indicados pelas previsões teóricas. 

Em trabalho recente [9], estudaram nanofluidos à base de água de ALO, e ZrO, 
foram caracterizados quanto à sua utilização promissora em aplicações de transferência de 
calor. Prepararam-se três concentrações diferentes de soluções dispersas de nanofluidos 
citados (0,01% vol., 0,05% vol., e 0,1% vol.) de nanofluidos comerciais. As medições 
experimentais foram realizadas a diferentes temperaturas. Condutividade térmica, 
viscosidade e densidade dos nanofluidos foram medidos. Concluiu-se que a concentração 


em volume, tamanho de partícula/forma e temperatura são variáveis importantes. 


4 | RESULTADOS E CONCLUSÃO 


Esta pesquisa iniciou-se em março deste ano (2016). No momento encontra-se em 
andamento a revisão bibliográfica, incluindo pesquisa dos nanofluidos mais adequados a 
serem adquirido para testes na bancada experimental. 

Devido ao potencial da utilização de nanofluidos em sistemas de refrigeração de 
reatores nucleares, mesmo com a pesquisa ainda em andamento, foi verificada a viabilidade 
dessa aplicação. Constatando diversos benefícios nessa possibilidade de uso, com isso 
pode-se melhorar as características da transferência de calor de fluidos habitualmente 
utilizados como refrigerante, adicionando nanopartículas constituídas de nanotubos de 
carbono. E atuando como líquido de refrigeração primária em sistemas de segurança, 
contribuindo para a atenuação de acidentes graves. 
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RESUMO: Atualmente, a eletricidade em todo 
o mundo é obtida de fontes não renováveis. 
No entanto, fontes limpas e renováveis têm 
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se mostrado mais atraentes para a geração 
de energia, porque representam uma solução 
potencial para o problema do aquecimento global. 
Portanto, há um aumento constante de usos 
de fontes alternativas de geração de energia, 
como o vento. No Brasil, a geração eólica ocorre 
principalmente com o uso de turbinas de grande 
porte, em regiões onde a velocidade média é 
superior a 7 m/s, sendo a geração eólica por meio 
de microgeradores ainda incipiente. Técnicas 
para a geração de energia eólica em escala micro 
(principalmente em locais de baixas velocidades 
de vento) precisam ser aprimoradas de modo a 
se obter uma maior eficiência de conversão e 
favorecer a viabilidade deste aproveitamento. 
Este estudo foi desenvolvido para avaliar o 
comportamento do escoamento do ar em difusor 
para aplicação na geração de energia eólica. 
Os modelos D2, D8, D16 e D13 foram os que 
promoveram maior aumento de velocidade do 
vento dentro do bocal, porém com desvantagem 
em relação ao modelo D18 devido a possuírem 
maior relação de comprimento L/D o que na 
prática pode contribuir para um maior custo de 
aplicação destes modelos. O modelo D16 se 
destacou por apresentar uma razão de potência 
de 8,86, e uma relação entre a velocidade máxima 
e a velocidade de corrente livre igual a 2,07. Para 
futuros estudos, pretende-se ainda trabalhar com 
melhoramento das geometrias, suavizando as 
curvas do bocal criando-se um bordo de ataque 
que permita o alinhamento do escoamento para 
minimizar os efeitos de separação do fluido, 
encurtando-se também o comprimento do 
bocal utilizando espécies de flanges na saída 
do dispositivo para melhorar aproveitamento de 
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area de seção de escoamento e evitar a recirculação do ar para dentro do bocal. 
PALAVRAS-CHAVE: Difusores aumentados, Dinâmica dos Fluidos Computacional, Potência 
Eólica. 


COMPUTATIONAL SIMULATION OF AIR FLOW IN DIFFUSER FOR 
APPLICATION IN WIND POWER GENERATION 


ABSTRACT: Currently, electricity worldwide is obtained from non-renewable sources. 
However, clean and renewable sources have been shown to be more attractive for power 
generation, because they represent a potential solution to the problem of global warming. 
Therefore, there is a constant increase in the use of alternative sources of energy generation, 
such as the wind. In Brazil, wind generation occurs mainly with the use of large turbines, 
in regions where the average speed is above 7 m/s, with wind generation through micro- 
generators still incipient. Techniques for generating wind energy on a micro scale (mainly 
in places with low wind speeds) need to be improved in order to obtain greater conversion 
efficiency and favor the viability of this use. This study was developed to evaluate the 
behavior of air flow in a diffuser for application in the generation of wind energy. The models 
D2, D8, D16 and D13 were the ones that promoted a greater increase in wind speed inside 
the nozzle, however with disadvantage in relation to the model D18 due to their greater L/D 
length ratio, which in practice can contribute to a greater cost of applying these models. The 
D16 model stood out for having a power ratio of 8.86, and a ratio between maximum speed 
and free current speed equal to 2.07. For future studies, it is also intended to work with the 
improvement of geometries, smoothing the nozzle curves creating a leading edge that allows 
the flow alignment to minimize the effects of fluid separation, also shortening the nozzle length 
using types of flanges at the outlet of the device to improve the use of the flow section area 
and avoid the recirculation of air into the nozzle. 

KEYWORDS: Diffuser augmented, Computational Fluid Dynamic, Wind Power. 


11 INTRODUCÃO 


O vento tem sido uma das principais fontes de geração de energia renovável no 
Brasil nos últimos anos. Regiões com velocidades de vento mais baixas também poderiam 
se beneficiar da geração de energia eólica em pequena escala e geração descentralizada 
usando micro turbinas com rotores adequados para partidas de baixa velocidade. Além 
disso, o uso de dispositivos auxiliares, como difusores, pode contribuir para o aumento da 
velocidade do vento que atravessa o rotor e, consequentemente, causar um aumento da 
energia eólica gerada. 

A exploração da energia eólica ocorre através da conversão da energia cinética 
contida nos ventos em energia mecânica e, posteriormente em energia elétrica, por meio 
de um gerador. A quantidade de energia cinética contida em um fluxo da massa de ar é 
função da velocidade do vento, da área varrida do rotor da turbina e da massa específica do 


vento. Na Figura 1 está representado um volume de controle contendo um rotor de turbina 


Projeto, Análise e Otimização na Área das Engenharias Capítulo 13 EE 


eólica localizado entre as seções 2 e 3 de um difusor. 


Jo 





Figure 1. Volume de controle para análise do escoamento do vento em um difusor. 


Considerando o vento com velocidade de corrente livre Voo, velocidade do vento na 
entrada do volume de controle que inclui a turbina U2, pressão P1 igual a P4 e igual a Paim 
local, e desprezando diferença de pressão entre as seções 2 e 3, a potência eólica P do 


vento pode ser calculada pela Equação 1. 


P=-pAU) (1) 


Onde p é a massa específica de ar (kg/m?), A é a área varrida pelo rotor da turbina 
(m?) e U, é a velocidade do vento que entra no rotor (m/s). 

Em busca da geração alternativa de energia eólica baseada nas regiões com 
menores velocidades do vento, vários autores colocaram em prática suas pesquisas para 
acelerar o vento antes que ele passasse pelas pás da turbina usando um dispositivo tipo 
difusor. A turbina contendo este dispositivo é do tipo DAWT (turbina eólica com difusor 
aumentado). Após o primeiro modelo proposto por Lilley e Rainbird (1956), vários estudos 
experimentais e computacionais foram realizados (Oman et al., 1975), (Foreman et al., 
1978), (Phillips et al., 2002). Estes estudos mostraram resultados experimentais que 
confirmam que as turbinas eólicas com difusor podem extrair muito mais energia do vento 
do que qualquer turbina sem difusor em condições de vento similar com o mesmo diâmetro 
do rotor. Estudos teóricos e numéricos foram realizados para investigar e melhorar o 
desempenho de DAWTs (Abe et al., 2005); (Ohya et al., 2008); (Jafari e Kosasih, 2014); 
(Wang et al., 2015). A maioria dessas pesquisas apresentou conclusões semelhantes com 
os resultados experimentais iniciais. 

Ao longo de anos, este dispositivo foi sendo modificado para se conseguir o máximo 
de aproveitamento do vento com o mínimo de custo, surgindo com isso um termo conhecido 


como “Wind Lens Turbine” para o sistema de geração eólica que inclui a turbina e um 
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difusor de comprimento bastante reduzido com um grande flange conectado na sua saída. 
Devido a formação de vórtice na parte posterior do flange, uma enorme quantidade de 
fluxo de massa pode ser escoada pela turbina. Trabalhos experimentais de vários autores 
demonstraram que o limite de Betz pode ser excedido quando se utiliza uma turbina com 
este tipo de dispositivo difusor flangeado (Oka et al. 2016). O limite de Betz indica que, 
mesmo para as turbinas mais eficientes (turbinas de 2 ou 3 pás de eixo horizontal), é 
possível recuperar apenas um máximo de 59% da energia do vento, o que significa que o 
coeficiente teórico máximo de potência é, aproximadamente, 0,59. 

Govindharajan et al. (2013), realizaram uma investigação numérica para a otimização 
no projeto de um difusor flangeado do tipo Wind Lens, que trata do efeito das regiões de 
baixa pressão geradas pelo dispositivo sobre a vazão mássica através da turbina eólica e 
aumento de torque. 

Kunjir et al. (2018) realizaram uma análise da performance de turbinas eólicas 
usando tecnologia de Wind Lens em que eles analisaram que a estrutura de anel com 
abas simples cria turbulência atrás das pás da turbina para atrair mais fluxo de ar para a 
turbina. A estrutura efetivamente aumentou a velocidade da turbina e, portanto, levou a um 
aumento de 20 a 25% da potência gerada. 

Considerando a necessidade de se ampliar o entendimento sobre como ocorre o 
escoamento do ar nestes dispositivos e o aumento na geração de potência destes sistemas, 
muitos trabalhos tem sido desenvolvidos utilizando ferramentas de CFD (Computational 
Fluid Dynamics). O uso de CFD, permite estudar e interpretar fenômenos de um escoamento 
de fluido e prever resultados, possibilitando a otimização de sistemas com significativa 
redução do número de experimentos e de custos para a construção de protótipos. 

Diante disso, com este trabalho tem-se, o objetivo de analisar por meio de 
ferramentas de CFD, o efeito de diferentes geometrias de difusores aumentados sobre o 


comportamento do escoamento internamente e ao redor destes dispositivos. 


21 COMPUTATIONAL PROCEDURE 


As etapas para o desenvolvimento de estudos utilizando ferramentas de CFD 
incluiram: 
1) Definição da geometria a ser analisada e do domínio computacional (limite físico 
onde as equações básicas da mecânica dos fluidos serão resolvidas). 
Uma representação genérica da geometria a ser avaliada está representada na 
Figura 1, porém foram avaliadas diferentes geometrias variando-se ângulos de divergência 
do bocal (a), e variando-se as relações de comprimento e diâmetro de entrada (L/D). Na 


Tabela 1 estão apresentados os parâmetros utilizados para a criação das geometrias. 
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Table 1. Parâmetros utilizados para a criação das geometrias. 





2) Geração de malhas, onde o volume ocupado pelo fluido é dividido em células 
discretas. 


3) Definição das condições de contorno, onde as condições limites foram definidas. 
Isso envolveu a especificação do comportamento do escoamento, considerado 
como incompressível e em regime permanente e das propriedades do fluido nos 
limites do problema. 


4) Simulação computacional, etapa em que as equações foram resolvidas 
iterativamente usando um modelo de turbulência. 


5) Pós processamento, etapa onde os dados numéricos e gráficos (resultados) de 
gradientes de velocidade e de pressão em todo o campo de escoamento foram 
analisados. 

As etapas deste estudo foram desenvolvidas utilizando o pacote Ansys Fluent 17.0 
versão acadêmica. 

Para a solução das equações o escoamento foi considerado isentrópico com 
velocidade de corrente livre de 4,0 m/s a montante e pressão manométrica igual a zero a 
jusante do escoamento. A massa específica do ar foi considerada de 1,2 kg/m, a viscosidade 
dinâmica de 1,79 x10º kg/m.s e velocidade zero nas superfícies interna e externa do bico 
respeitando a condição de parede ou seja a condição de não escorregamennto. 
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O domínio do modelo foi construído de acordo com o mostrado na Fig. 2, ou seja, 
baseado em um plano secional 2D longitudinal com o eixo de simetria no centro do bocal 
divergente. 


16,5 M 





Figura 2. Domain of the model. 


A discretização do domínio foi baseada em uma malha estruturada com elementos 
quadrilaterais. Para a malha padrão, foi gerado um total de 4 malhas com número de 
elementos diferentes e, em seguida, foi realizada uma análise de convergência da malha 
em função do gradiente de velocidade no centro do bocal, na região central da simetria 


conforme mostrado na Figura 3. 





Figura 3. Domínio (a) e malha (b) gerados para o modelo computacional 
Na etapa de simulação, foi utilizado o modelo de turbulência k-w SST com 5% de 


intensidade de turbulência e 0,45 de escala de turbulência. O algoritmo SIMPLE foi utilizado 


para tratar o acoplamento velocidade e pressão devido ao conhecido bom desempenho em 
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uma grande gama de problemas. Para a discretização espacial das equações de transporte 
utilizou-se o esquema Upwind de segunda ordem. O critério de convergência de 10% foi 


aplicado. 


31 RESULTADOS E DISCUSSÃO 


Os gradientes de velocidade obtidos por simulação computacional do escoamento 


em bocal divergente com ângulo de 5º estão apresentados na Figura 4. 





Figura 4. Gradientes de velocidade para os modelos com ângulo de divergência de 5º. 


Para os modelos de bocais com ângulo de divergência de 5º observa-se um aumento 
na velocidade do fluido na região interna do bocal divergente em relação a velocidade de 
corrente livre (4,0 m/s) para todas as relações L/D. O modelo que proporcionou maior 
aumento da velocidade do fluido dentro do bocal de ângulo de 5º foi o Difusor 2 cuja 
relação L/D é igual a 2,5, seguido do Difusor 1 com relação L/d igual a 3. A relação entre a 
velocidade máxima do fluido dentro do bocal e a velocidade de corrente livre para o Difusor 
2 foi de 1,8, indicando que este modelo promove um aumento na velocidade de corrente 
livre de 80%. 

De acordo com os resultados obtidos por simulação, para relação L/D inferior a 2,5, 
a diminuição do comprimento do bocal, causa uma redução no aumento de velocidade, ou 
seja, quanto menor o comprimento do bocal, menor o aumento da velocidade na região 


interna do mesmo. 
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Figura 5. Gradientes de velocidade para os modelos com ângulo de divergência de 10º. 
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Para os modelos com ângulo de divergência em 10º, verificou-se que existe uma 
variação de acordo com a relação L/D, na posição em x onde ocorrem os aumentos de 
velocidade dentro do bocal, conforme mostrado na Figura 5. 

Possivelmente, isso ocorre devido à presença de reentrada de fluido (efeitos de 
vorticidade) a jusante do bocal. Isso se deve ao fato de que a recirculação é mais intensa 
na parte posterior onde ocorre separação do escoamento. 

Os modelos D7 e D8 foram os que apresentaram maior aumento na velocidade do 
vento. O melhor caso para este ângulo foi o modelo D8, onde sua relação de comprimento 
L/D também é 2,5 assim como observado para os modelos com ângulo de difusão de 
5º. Para este modelo de difusor, a relação entre a velocidade máxima do fluido dentro 
do bocal e a velocidade de corrente livre foi de 1,45. Os demais casos apresentaram 
pequenas regiões com velocidade aumentada e, portanto, podem contribuir bem menos 
para o aumento da potência eólica a ser gerada por uma turbina eólica instalada nestes 
dispositivos. 







PELDL 
—— 
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Figura 6. Gradientes de velocidade para os modelos com ângulo de divergência de 15. 
Os gradientes de velocidade obtidos para os bocais com ângulo de 15º estão 


apresentados na Figura 6. Observa-se um ganho em todos os difusores com ângulo de 
difusão de 15º, porém os difusores D13 (L/D = 3), D16 (L/D = 1,5) e D18 (L/D = 0,5) se 
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destacaram apresentando maior aumento na velocidade de escoamento. 

Nesta categoria de ângulo de divergência, o modelo D16 se destacou por possuir 
uma velocidade máxima de 8,3 m/s, porém isso ocorreu na região próxima a jusante do 
bico. O caso D13 trouxe bons resultados quando comparado com o também, porém, possui 
uma relação comprimento superior de L/D. A redução de L/D dos bocais D13 para D16 
proporcionaram algumas interferências na magnitude da velocidade, com diferença na 
localização das velocidades máximas dentro do bocal. Este fato pode estar associado ao 
aumento da vorticidade resultante da mudança na geometria a jusante do bocal. 

Fazendo uma análise de todos os modelos simulados, os modelos D2, D8, D16 e 
D13 são os melhores casos de todos os modelos analisados. Comparando esses modelos, 
o destaque para os resultados dos casos D2, D13 e D16, por possuírem uma maior área 
de uso com velocidade aumentada para aproveitamento eólico. Em um caso isolado, o 
modelo D18 apresenta uma pequena região onde possui a velocidade máxima, porém, 
suas velocidades para a área total atingem uma média aproximada em 7,1 m/s ao longo 
do raio do bocal, o que pode ser um bom resultado considerando sua menor relação 
de comprimento (L/D = 0,5) o que é uma vantagem quando comparado com os outros 
modelos pois pode favorecer a redução significativa de material de fabricação e também 
reduzir custos associados a ancoragem destas estruturas. Considerando a densidade do 
ar igual a 1,2 kg/m3 para a área circular com o diâmetro unitário, a potência eólica para 
uma turbina idealizada sem o bocal divergente instalado (calculado conforme Equação 1) e 
considerando ainda a velocidade do fluxo livre de 4 m/s seria de 30,2 W. A Tabela 2, mostra 


as potências máximas teóricas para os modelos D2, D8, D13, Di6 e D18. 


Potencia o 


e» sé ss 





Tabela 2. Resultados para potência eólica teórica dos difusores que promoveram maior 
aumento de velocidade do vento. 


De acordo com os valores obtidos neste trabalho, os resultados para o grupo de 
divergência de 5º aproximam-se dos analisados por (Ohya et al., 2008a) e (Matsushima 
et al., 2006a), que apresentam os melhores resultados para o bocal difusor próximo ao 
ângulo de 5º e 6º. No trabalho do (Chaker et al., 2012), são apresentadas discussão sobre 
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os resultados do bocal com o ângulo de divergência de 5º a 15º, neste trabalho os autores 
relatam que para ângulos maiores que 5º a aplicação de um flange na jusante do bocal 
tornou-se atraente. Com o flange, o melhor resultado encontrado por esses autores são 12º 
de ângulo de divergência. Portanto, o uso deste flange afeta na região de deslocamento 
onde ocorre o efeito de recirculação, para evitar que o fluido retorne após a saída do bocal, 
isso promove o ganho de velocidade e garante que seus grandes gradientes de velocidade 
possam ser deslocados para a entrada do bico, melhorando o aproveitamento da área da 
turbina do rotor. 

No presente trabalho, todos os modelos avaliados mostraram a separação do 
escoamento do fluido na região interna do difusor próximo à parede. Este efeito é uma 
função da mudança abrupta da borda de ataque e seu ângulo de divergência. Os modelos 
de grupo do ângulo 10º e 15º possuem uma grande região de divisão de fluxo e isso 
interfere na área útil para a geração de energia, portanto, um estudo deste efeito também 
nas superfícies desse tipo de bocal buscando a suavização na entrada do difusor terá de 
ser considerado. 

Mansour e Meskinkhoda (2014) em estudo utilizando a simulação computacional 
para predizer o comportamento do fluxo do vento em bocal com o ângulo difusor de 4º e 
relação L/D de 1,5 obtiveram resultados da razão de velocidade em torno de 1,3 também 
sem flange na saída do bocal. Para a relação angular de 5º e L/D de 1,5, a relação de 
velocidade observada no presente estudo foi próxima a 1,52. Esses autores verificaram 
que é possível aumentar a razão de velocidade de 1,3 para aproximadamente para 1,7 com 
a mesma relação de comprimento usando um tamanho adequado de flange na saída do 
bocal. Um aumento significativo na velocidade do vento foi obtido com a geometria do caso 
D16, onde a taxa da velocidade na região interna do bocal chegou a 2,07. Para este caso 
a relação entre a potência teórica em uma turbina sem difusor e uma turbina com difusor é 
de aproximadamente 8,86. 

Abe et al., (2005), e Ohya et al., (2008) em análises numéricas e experimentais, 
observaram que a potência para turbinas com o difusor é quase 4 vezes maior para a razão 
da velocidade do 1,5 e 1,6. O aumento de potência observado com a utilização da turbina 
eólica no bico deve-se, entre outros fatores, ao incremento da vazão mássica ocasionado 
pela redução da pressão na saída da vazão. A diferença observada pode ser explicada no 
primeiro estágio como a variação da relação de velocidade. A diferença nessas relações de 
velocidade será observada quando a simulação ocorrer com a turbina eólica. Além disso, 
estima-se que para resultados experimentais as razões de velocidade sejam menores do 


que casos numéricos apresentadas neste trabalho. 


3.1 Efeito do ângulo do difusor sobre o escoamento 


Para avaliar o efeito do ângulo do difusor sobre a velocidade do escoamento do 
vento dentro do bocal, os casos foram agrupados mantendo-se constante as relações de 
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comprimento, ou seja, comparando-se o caso D1 com o D7 e D13, D2 com o D8 e D14, D3 
com o D9 e D15, D4 com o Di0O e D16, D5S com o D11 e D17, DO como D1i2 e D18. 

De acordo com os resultados mostrados, conforme esperado observa-se que O 
ângulo de abertura do difusor tem efeito sobre os gradientes de velocidade. Para ângulo 
de 5º é possível observar gradientes de velocidade bem definidos com as máximas 
velocidades ocorrendo na região do centro do bocal em praticamente todas as relações 
de comprimento. São perfis de velocidade em desenvolvimento, ou seja, com variações 
observadas também na direção x. Comportamento semelhante é observado para ângulo de 
10 e relações de comprimento de 3 e 2,5 e para ângulo de 15 com relação de comprimento 
de 3. Para os ângulos de abertura de 10 e 15º, observa-se uma maior espessura de 
escoamento separado na parte interna do bocal e consequentemente menor área de seção 
de escoamento do fluido quando comparado com os casos de ângulo de 5º, minimizando o 
aumento de velocidade em grande parte do bocal. 

Apesar disso em algumas regiões para os casos com ângulo de 15º e relações de 
comprimento de 3 e de 1,5 observou-se aumento expressivo da velocidade do fluido em 
relação à velocidade de corrente livre. 


3.2 Efeito do comprimento do bocal sobre a velocidade do escoamento do 
vento 

Para avaliar o efeito do comprimento do bocal sobre a velocidade do escoamento, 
os gradientes de velocidade foram organizados mantendo-se o ângulo do difusor constante 
variando-se a relação de comprimento do bocal. 

A diminuição da relação de comprimento do bocal causa efeito sobre o escoamento, 
independente do ângulo de abertura, porém quando o ângulo de abertura é de 10 e de 15. 
este efeito é somado com o efeito da pressão externa em torno do bocal. Com a diminuição 
do comprimento do bocal, aumenta-se o efeito dos gradientes de pressão externa sobre 
o escoamento interno ao bocal. Quanto menor a relação de comprimento, menor é a 
região com máximas velocidades, provavelmente devido ao efeito de recirculação do fluido 
na saída do bocal. Para minimizar os efeitos do escoamento externo sobre o interno, a 
colocação de um flange na saída do bocal tem sido testada com resultados interessantes 
(Abe et al., 2005a), (Ohya et al., 2008a), porem o uso deste tipo de dispositivo não foi objeto 
de avaliação neste trabalho. 


41 CONCLUSÃO 


Para ângulo de 5º é possível observar gradientes de velocidade bem definidos com 
as máximas velocidades ocorrendo na região do centro do bocal em praticamente todas as 
relações de comprimento. Para os ângulos de abertura de 10 e 15º, observa-se uma maior 
espessura de escoamento separado na parte interna do bocal e consequentemente menor 
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àrea de seção de escoamento do fluido quando comparado com os casos de ângulo de 5º, 
minimizando o aumento de velocidade em grande parte do bocal. 

A diminuição da relação de comprimento do bocal causa efeito sobre o escoamento, 
independente do ângulo de abertura, porém quando o ângulo de abertura é de 10 e de 15. 
este efeito é somado com o efeito da pressão externa em torno do bocal. Quanto menor 
a relação de comprimento, menor é a região com máximas velocidades, provavelmente 
devido ao efeito de recirculação do fluido na saída do bocal. 

Os modelos D2, D8, D16 e D13 foram os que promoveram maior aumento de 
velocidade do vento dentro do bocal, porém com desvantagem em relação ao modelo D18 
devido a possuírem maior relação de comprimento L/D o que na prática pode contribuir 
para um maior custo de aplicação destes modelos. O modelo D16 se destacou por 
apresentar uma razão de potência de 8,86, e uma relação entre a velocidade máxima e a 
velocidade de corrente livre igual a 2,07. Para futuros estudos, pretende-se ainda trabalhar 
com melhoramento das geometrias, suavizando as curvas do bocal criando-se um bordo de 
ataque que permita o alinhamento do escoamento para minimizar os efeitos de separação 
do fluido, encurtando-se também o comprimento do bocal utilizando espécies de flanges 
na saída do dispositivo para melhorar aproveitamento de área de seção de escoamento e 
evitar a recirculação do ar para dentro do bocal. 
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RESUMO: Rochas, solos e águas possuem 
radionuclídeos de ocorrência natural. Dentre 
estes, destaca-se o “8U e 2Z2Th e seus 
produtos de decaimento. A distribuição de tais 
radioelementos em termos de concentração de 
atividade depende do tipo mineral e da origem. 
Dessa forma, todo processamento técnico de 
minerais resulta na liberação de radionuclídeos 
de meia-vida longa e curta, principalmente de 
radônio e sua progênie. Assim, é importante 
monitorar esse gás e seus produtos em minas 
subterrâneas a fim de avaliar o risco radiológico 
para os trabalhadores  ocupacionalmente 
expostos. Inserido nesta preocupação que é 
global, medidas da concentração de radônio e 
da sua progênie, cálculo do fator de equilibrio 
e estimativa de dose foram realizadas para 5 
minas subterrâneas do Brasil. A concentração de 
radônio e da sua progênie foram determinadas 
utilizando os monitores contínuos AlphaGUARD 
PQ2000PRO (Saphymo) e DOSEman Pro 
(Sarad), respectivamente. O estado de equilíbrio 
entre o radônio e sua progênie foi avaliado. 
Baseado nestes dados, estimou-se a dose 
efetiva anual para os mineiros de acordo com 
a metodologia recomendada pela UNSCEAR, 
2000. Assim, a Concentração de Radônio 
variou de 60 a 4964 Bq.m* e a Concentração 
Equivalente de Equilíbrio se encontra faixa de 12 
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a 1174 Bq.m*. O valor médio dos fatores de equilíbrio obtido para as 8 minas foi de 0,4 
(0,2 — 0,7). Com isso, os valores de dose efetiva para os trabalhadores das minas medeiam o 
seguinte intervalo: 0,2 a 21 mSv.a!, com valor médio de 9,8 mSv.a”. Portanto, os resultados 
mostram a relevância de avaliar continuamente e permanentemente o comportamento do 
radônio e da sua progênie e de adotar medidas de segurança contra a radiação natural no 
ambiente de minas subterrâneas. 

PALAVRAS-CHAVE: Radônio, Fator de Equilíbrio e Dose. 


EQUILIBRIUM FACTORS AND DOSES IN BRAZILIAN UNDERGROUND MINES 


ABSTRACT: Rocks, soils and waters have naturally occurring radionuclides. Among these, 
the 28U and 22Th and its decay products stand out. The distribution of such radioelements 
in terms of activity concentration depends on the mineral type and origin. Thus, all technical 
processing of minerals results in the release of radionuclides of long and short halíf-life, mainly 
radon and its progeny. Thus, it is important to monitor this gas and its products in underground 
mines in order to assess the radiological risk for occupationally exposed workers. In this global 
concern, measurements of radon concentration and its progeny, calculation of the equilibrium 
factor and dose estimation were performed for 5 underground mines in Brazil. The radon 
concentration and its progeny were determined using the continuous monitors AlphaGUARD 
PQ2000PRO (Saphymo) and DOSEman Pro (Sarad), respectively. The state of equilibrium 
between radon and its progeny was evaluated. Based on these data, the annual effective dose 
was estimated for miners according to the methodology recommended by UNSCEAR, 2000. 
Thus, the Radon Concentration ranged from 60 to 4964 Bq.m* and the Equivalent Equilibrium 
Concentration ranged from 12 to 1174 Bq.m?. The mean value of the equilibrium factors 
obtained for the 8 mines was 0.4 (0.2 — 0.7). As a course, the effective dose values for mine 
workers measure the following range: 0.2 to 21 mSv.a”, with an average value of 9.8 mSv.a”. 
Therefore, the results show the relevance of continuously and permanenily evaluating the 
behavior of radon and its progeny and of adopting safety measures against natural radiation 
in the environment of underground mines. 

KEYWORDS: Radon, Equilibrium Factor, Dose. 


11 INTRODUÇÃO 


Os radionuclídeos de ocorrência natural, urânio e tório seus produtos de decaimento, 
estão presentes em concentrações variadas nos materiais geológicos, especialmente 
em rochas e solos. Os níveis específicos desses radioelementos estão relacionados ao 
tipo e origem das rochas (UNSCEAR, 2000). Sabe-se que são litófilos e se concentram 
preferencialmente nas rochas ígneas ácidas, entre as quais se encontram o granito, e não 
nas básicas e ultrabásicas (BONOTTO E SILVEIRA, 2006). As rochas sedimentares, em 
geral, contêm níveis baixos de radioatividade natural. Entretanto, existem algumas exceções 
como os arenitos e as rochas fosfáticas que com frequência apresentam concentrações 
relativamente altas de radionuclídeos naturais (UNSCEAR, 2000). 

Durante a abertura de minas subterrâneas, radionuclídeos de meia-vida longa e 
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curta provenientes das séries naturais do urânio e tório são liberados e concentram-se 
em tais ambientes tornando poluidores potenciais (OTHMAN et al.,1992). Dentre esses, 
destaca-se o radônio e sua progênie que são considerados a principal fonte de exposição 
a radiação em minas subterrâneas (ICRP, 1986; BALDIK et al., 2006). 

No que tange à sua formação, o radônio é produzido pelo decaimento do rádio 
presente no corpo mineral e por recuo emana da matriz sólida dos materiais para os 
interstícios dos grãos preenchidos com ar ou água ou para os planos de fraturas das rochas 
(SANTOS, 2008). Logo após ser emanado, o gás tende a migrar-se por difusão molecular ou 
fluxo convectivo no sentido da superfície, onde a porção que não decaiu no percurso exala 
para o interior das galerias. Outro mecanismo de transporte do radônio para os ambientes 
das minas é a circulação de águas cuja concentração de radônio depende dos tipos de 
rochas, do volume relativo de água e da temperatura. A importância de cada processo 
de transferência relaciona-se com as fraturas tectônicas e geológicas da formação e do 
comportamento hidrogeológico do aquífero nas operações de mineração (ICRP, 1986). 

O radônio, depois de exalado, migra ao longo das correntes de ventilação 
enquanto forma seus produtos de decaimento de meia-vida curta: “Po, “Pb, “Bi e 
*4Po. O radônio e sua progênie no interior de minas subterrâneas é uma preocupação 
internacional (OTHMAN et al., 1992; BALDIK et al., 2006; FATHABADI, 2006; CILE, 2010; 
EVANGELISTA, 2002; PAGE E SMITH, 1992; ANJOS et al, 2010; THINOVA E ROVENSKA, 
2011). A dose correspondente ao radônio é atribuída à sua progênie que quando inalada 
deposita-se no pulmão, especialmente no trato respiratório superior, e irradia o tecido 
pulmonar ao decair. Estudos epidemiológicos revelam uma forte correlação entre câncer 
de pulmão e exposição ao radônio. Assim, a Organização Mundial da Saúde classifica o 
radônio como um carcinógeno humano (ICRP, 2010). Diante disso, torna-se imprescindível 
avaliar a distribuição e a origem do radônio em minas subterrâneas visando mitigação e o 
atendimento aos padrões de proteção radiológica. 

Diante disso, torna-se imprescindível avaliar a dose devida aos produtos de 
decaimento do radônio, a qual os mineiros estão expostos. Para estimar a dose, em geral, 
mede-se a concentração de radônio e assume-se um Fator de Equilíbrio (F) entre o radônio 
e sua progênie usualmente 0,4 ou 0,5. Entretanto, na realidade, as concentrações de 
radônio e sua progênie variam significantemente com o tempo e o espaço (YU et al). Alguns 
fatores contribuem para a variabilidade da concentração desses radionuclídeos, são eles: a 
formação geológica, o teor de urânio, tório e rádio, a taxa de exalação de radônio devido às 
diferenças de pressão, o grau de diluição, a porosidade, a umidade, o tipo de mina, a idade 
do ar dentro da mina e as condições de trabalho, particularmente o grau de ventilação 
(ICRP, 1986; BALDIK et al, 2006; FATHABADI, 2006; CILE, 2010). Contudo, a adoção de 
um valor para F pode não refletir as condições reais de uma mina (YU et al., 2008). Esse 
problema pode ser resolvido através de determinações experimentais de tal fator. 

Neste contexto, o presente trabalho objetiva determinar o fator de equilíbrio 
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a fim de avaliar a dose em 8 minas subterrâneas no Brasil. As minas incluem extração 
de algamatolito, carvão, esmeralda, fluorita, turmalina, schelita e zinco. Este trabalho é 
resultado do “Projeto Radônio” realizado pela Comissão Nacional de Energia Nuclear e os 
resultados obtidos contribuirão para as discussões necessárias à elaboração das normas 


ou posições regulatórias. 


21 METODOLOGIA 


2.1 Concentração de Radônio no Ar 


As medições da concentração de radônio no ar foram executadas com o detector 
AlphaGUARD PQ2000PRO (Saphymo GmbH). O AlphaGUARD opera segundo os 
princípios de uma câmara de ionização pulsada. A técnica consistiu em instalar o 
AlphaGUARD no ambiente interior das minas no modo difusão por aproximadamente 2 
dias de acordo com o protocolo para medições de curto prazo (UNSCEAR, 1993). Neste 
trabalho, o equipamento foi programado para intervalos de amostragem de 10 minutos 
ou eventualmente de 60 minutos, atuando como um detector passivo contínuo. Os 
pontos amostrados consistiram de frentes de trabalho, galerias e locais onde os mineiros 
permaneciam a maior parte do tempo e o equipamento encontrava-se seguro. Além da 
concentração de radônio, informações de data, hora, temperatura, umidade relativa e 
pressão barométrica foram também registradas. O limite de detecção do AlphaGUARD é 2 
Bq.m*. Os dados foram obtidos e tratados com o software DataExpert (Saphymo GmbH, 
2010). 


2.2 Progênie do Radônio no Ar 


As análises por espectrometria alfa, para determinação da progênie do radônio, 
foram realizadas utilizando o detector DOSEman Pro (SARAD). Grandezas e unidades 
especiais são usadas para caracterizar a concentração da progênie do radônio no ar como: 
Concentração de Energia Alfa Potencial (CEAP) e Concentração Equivalente de Equilíbrio 
(CEE). A CEAP é definida como sendo a soma das energias alfa potenciais emitidas pela 
progênie do radônio por unidade de volume de ar. Essa grandeza pode ser expressa em 
termos da CEE, que corresponde à concentração de radônio em equilíbrio radioativo com 
sua progênie de meia-vida curta, que liberaria a mesma CEAP que a mistura real não 
equilibrada (ICRP, 1993). Neste equipamento, amostras de particulados de aerossóis 
são coletadas, por sucção de ar, em um filtro, que é analisado continuamente por um 
detector de silício com respeito à energia de decaimento alfa. Utilizaram-se intervalos de 
amostragem de 60 minutos. 


2.3 Fator de Equilíbrio 


Para a determinação do fator de equilíbrio, o AlphaGUARD foi instalado juntamente 
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com o DOSEman nas minas subterrâneas estudadas. Dessa forma o fator de equilíbrio 
(F) foi calculado através da seguinte equação (UNSCEAR, 2000): F=CEE /C. Onde C, 
é a concentração de radônio no ar. O F caracteriza o desequilíbrio entre radônio e sua 
progênie em termos da energia alfa potencial. Este fator tem sido amplamente estudado 
em condições distintas para avaliar a dose que as pessoas estão expostas para uma dada 
concentração de radônio (CIGNA, 2005). Normalmente, quando não se tem informação ou 
resultados experimentais, recomendam-se a adoção de um fator de equilíbrio médio de 0,4 
(UNSCEAR, 2000; ICRP, 1993). 


2.4 Dose efetiva anual 


A dose efetiva anual (H), por sua vez, foi estimada de acordo com a metodologia 
descrita no relatório UNSCEAR 2000, na qual: H = C.FTk. Onde T é o tempo de 
permanência em locais de trabalho, T= 2000 h.a (ICRP, 1993), e ké o fator de conversão, 
k= 9 nSv (Bq.h.m*)* (UNSCEAR, 2000) 


31 RESULTADOS 


Os valores médios da Concentração de Radônio (C, ), da Concentração Equivalente 
de Equilíbrio (CEE), do Fator de Equilíbrio (F) e da Dose Efetiva Anual devido à progênie do 
radônio (H) para as 8 minas subterrâneas brasileiras são apresentadas na Tabela 1. Como 
pode ser observado na Tabela 1, obteve-se valores médios C e CEE na faixa de 60 a 4964 
Bq.mº* e de 12 a 1174 Bqg.m*º, respectivamente. Os dados encontrados mostram que as 
concentrações de radônio e da sua progênie variam consideralvelmente de mina para mina. 
Isso decorre da diferença de atividade específica dos radionuclídeos naturais precursores 
do radônio nas rochas e das diferentes condições de atividades operacionais das minas. 
Neste caso, destaca-se como principal fator a ventilação que dispõe de capacidade variada 
para cada mina estudada (FATHABADI, 2006). Baseado nesses dados, verifica-se que 
em algumas minas a C, excedeu o nível de referência de 1000 Bq.m* recomendado pela 
Comissão Internacional de Proteção Radiológica (ICRP, 2010) e pela Comissão Nacional 
da Energia Nuclear (CNEN, 2005) para locais de trabalho. Em três minas, concentrações 
médias de radônio elevadas estão presentes, são elas: 3889 Bq.m* na Mina A na condição 
“parada”, 4964 Bq.m* na Mina D e 1142 Bqg.m* na mina E. 

O F calculado para as 5 minas variou de 0,2 — 0,7 (Tabela 1). Isso demonstra que 
uma ampla faixa de valores de F ocorre nas minas, o que indica que adotar um valor médio 
para tal fator a fim de avaliar dose para casos individuais pode proporcionar grandes erros 
como mencionado em (FATHABADI, 2006). Nota-se que valores de equilíbrio diferentes 
foram encontrados para a Mina A na condição “em operação” e na condição “parada”, os 
quais são: 0,3 e 0,5 respectivamente. Esse fato pode ser explicado visto que atividades 
recorrentes de mineração como detonação, por exemplo, geram grandes quantidades de 
aerossóis e de liberação de radônio no interior das galerias. Além disso, C, diferentes 
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também foram observadas para tal mina na condição “parada” e na condição “em 
operação”. Estas discrepâncias estão claramente influenciadas também pela ventilação 
ativa na mina em operação. 

No que se refere à dose efetiva para os trabalhadores das minas, os valores 
medeiam o seguinte intervalo: 0,2 a 21 mSv.a! de acordo com a Tabela 1. O valor de 
referência de 1000 Bq.m* recomendado tanto pela CNEN (CNEN, 2005), como pela ICRP 
(ICRP, 2010) para locais de trabalho corresponde a uma dose efetiva anual de 7 mSy, 
considerando um fator de equilíbrio de 0,4 e 2000 horas de trabalho anual (metodologia 
UNSCEAR, 2000). Com base nisso, verifica-se na Tabela 1 que em duas minas (mina A na 
condição “parada” e mina D), a dose efetiva referente ao radônio está acima desse valor 
de referência. Além disso, a dose efetiva devido ao radônio para os mineiros dessas minas 
excede o limite de 20 mSv.a! aplicado para indivíduo ocupacionalmente exposto no Brasil 
(CNEN, 2014). Isso impõe a necessidade de controle de doses ou de aperfeiçoamento dos 
sistemas de ventilação. 


C,. (Bq.m*) CEE (Bq.m*) 
Condição da mina [média [média 
(max-min)J (max-min)P 


Em Operação 3889 (528-4128) 1174 (431317) 
rá 2992088) | 5T7 2a 
| MinaB | Em Operação 949 (584-1264) 259 (1293-444) 


mina D | EmoOporação | 406s (Is02-10860) | ri4o(2062765) [02 [ ar | 
MinaE | Panda | smoçmcoso) | ao(oasm) [02[ 4 | 
minar | EmOperação | sertossim | sastioar, [04[ 
minaG | EmoOporação | Goo(ra-S00) | Bioma meo) [07[ 4 + 
ima | EmOporação | cogo-5) | 16-56) [02| 02. 


a medidas contínuas. 





> Presença de atividades no interior da mina; sistema de ventilação desligado. 


Tabela 1: Concentração de Radônio, Concentração Equivalente de Equilíbrio, Fator de 
Equelíbrio e Dose Efetiva Anual devido à progênie do radônio para 8 minas subterrâneas no 
Brasil. 


41 CONCLUSÃO 


De acordo com os resultados mostrados e discutidos acima, somente em três minas 
subterrâneas (Mina A na condição “parada”, Mina D e Mina E), as concentrações de radônio 
obtidas excederam os limites recomendados pela ICRP. Diante disso, ações remediadoras 
devem ser adotadas em tais minas como, por exemplo, ventilação mecânica mais intensa 


observando o princípio de otimização de doses estabelecido pela CNEN (CNEN, 2014) e/ 
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ou escalas de trabalhos adequadas com o objetivo de reduzir o tempo de exposição dos 
trabalhadores. Os fatores de equilíbrio obtidos variaram amplamente de mina para mina. 
O valor médio de 0,4 foi determinado, com um mínimo de 0,2 e um máximo de 0,7. Dessa 
forma, ressalta-se que a utilização de um fator médio a fim de determinar dose não é 
recomendada. Embora o valor médio encontrado esteja coerente com o adotado pela ICRP 
e pela UNSCEAR. A dose para os trabalhadores das minas brasileiras abrangeu o intervalo 
de 0,2 a 21 mSv.a!. Consequentemente, na Mina A parada e na Mina D constatou-se que 
os valores de dose ultrapassaram o limite obtido com a metodologia da UNSCEAR e o 
limite para indivíduos ocupacionalmente expostos no Brasil, o que reforça e esclarece a 
necessidade de mitigação. Contudo, existem diversas minas subterrâneas no Brasil que 
ainda serão investigadas. Trabalhos futuros serão realizados para melhor compreensão do 
comportamento do radônio no ambiente subterrâneo. É importante caracterizar as fontes 
de radônio de cada mina, entender a circulação do radônio e de sua progênie dentro desses 


ambientes e avaliar o risco correspondente. 
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ABSTRACT: Vertebroplasty and Kyphoplasty 
have been studied for several years as 
percutaneous procedures for treatment of bone 
fractures, osteoporosis and other abnormalities in 
the skeletal system. Currently, these procedures 
have already been established as effective 
minimally invasive surgical procedures very 
promising in orthopedics and traumatology. 
Those procedures are guided by fluoroscopy 
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VERTEBRAE 


using radiopaque substances such as barium 
sulfate (BaSO,). The improvements of the 
radiological contrast at the image is still 
necessary. In present study, we addressed the 
benefits of Hydroxyapatite (HAp) as radiopaque 
element in the bone cement and the possibility of 
removing barium sulfate to reduce the toxicity of 
the material. 

KEYWORDS: Bone cement, Hydroxyapatite, 
Radiological contrast. 


DISTRIBUIÇÃO FÍSICA E CONTRASTE 
RADIOLÓGICO DE CIMENTOS 
IMPLANTADOS NAS VÉRTEBRAS IN 
VITRO 


RESUMO: La Vertebroplastia e Cifoplastia foram 
estudadas por vários anos como procedimentos 
percutâneos para tratamentos de fraturas 
ósseas, osteoporose e outras anomalias 
no sistema esquelético. Atualmente, esses 
procedimentos já tem sido estabelecidos como 
procedimentos cirúrgicos minimamente invasivos 
muito promissórios na ortopedia e traumatologia. 
Aqueles procedimentos estão guiados por meio 
da fluoroscopia usando substâncias radiopacas 
tais como sulfato de Bário BaSO,. As melhoras 
do contraste radiológico na imagem é ainda 
necessário. No presente estudo, encaminhamos 
os benefícios do Hidroxiapatita (HAp) como 
o elemento radiopaco no cimento ósseo e a 
possibilidade da remoção do sulfato de Bário 
para reduzir a toxicidade do material. 
PALAVRAS-CHAVE: Cimento 
Hidroxiapatita, Contraste radiológico. 
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11 INTRODUCTION 


A procedure was developed and first described in France in 1987 known as 
Vertebroplasty used as a percutaneous technique by means of an acrylic injected with a 
syringe into the fractured vertebra body. In the surgical context to this procedure, itis a 
minimally invasive technique for ambulatory nature. The technique was initially implemented 
for the treatment of vertebral hemangiomas and later was applied for decompression of 
fractures of osteolytic metastases and osteoporosis. Another percutaneous technique 
was introduced with use of PMMA (polymethylmethacrylate with chemical nomenclature 
[CH,C(CH,)(CO,CH.)] ) as bone cement to join the fractured bone paris. In addition, it has 
been inserted into a balloon that it is inflated in order to restore the anatomical structure 
of the vertebra. Such percutaneous technique is known as Kyphoplasty. One of the major 
problems of the percutaneous techniques is the possible extravasations of cement outside 
the bone region where it aims to be located. One solution to this undesirable event is 
the decompression of the vertebra by re-intervention by surgery (CHEN; LUO; ZHANG; 
NALAJALA et al., 2013; GERSZTEN, 2007; HENDRICKSON; SHEHATA; KIRCHNER, 
1976). These percutaneous procedures are complementary to therapeutic treatments 
within tumor control (KANEKO; SEHGAL; SKINNER; AL-GHAZI et al., 2010). In most cases 
patients have advanced metastatic disease stage and the radiation therapy is palliative. 
Some studies have proposed IMRT (Intensity Modulated Radiotherapy) or IGRT (Image- 
Guided Radiotherapy) as therapeutic treatments with less neurological impact on the 
marrow. Although demyelinating phenomenon of nerve tissues may occur since such 
structures are highly radiosensitive; it presents in lesser extension with IMRT and IGRT 
(EMAMI; LYMAN; BROWN; COLA et al., 1991; RYU; FANG YIN; ROCK; ZHU et al., 2008; 
RYU; JIN; JIN; ROCK et aíl., 2007; YAMADA; BILSKY; LOVELOCK; VENKATRAMAN et al., 
2008). At 2006, other options have been proposed in the field of radiation therapy based on 
the idealization of radioactive bone cement by our research group (MONTANO; CAMPOS: 
LEMOS; YOSHIDA et al., 2020). Such radioactive cement was addressed to treat bone 
metastases dealing with the elimination of the clonogenic properties of the cancer cells 
together with the reinforcement of the bone structure. At 2009, it has been followed by 
the use of Quadramet ('“Sm-EDTMP) administered percutaneously in the affected region 
by means of protocol Khyphoplasty (Kyphon) (ASHAMALLA; CARDOSO; MACEDON; 
GUIRGUIS et al., 2009). 

There are a number of biophosphonates considered bioceramics which have already 
been studied. Among those that can highlight is the Hydroxyapatite (HAp with chemical 
nomenclature [Ca,(PO,) (OH)]), contains 69% of its composition equivalent to naturally 
bone and the 20% of collagen matrices fibers as connective tissue. Therefore, since HAp 
is part of the nature of the bone, it has already been proposed to be mixed with the PMMA 
reaching optimal adhesion's response of osteoblasts in the biomaterial (DOROZHKIN, 
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2009; 2010; HIRSCH; ROSENSTEIN; MEDICH; MARTEL et al., 2009). 

Considering the cement implant, the improvement of the radiological contrast at the 
image is still necessary. In present study, we addressed the benefits of HAp as radiopaque 
element in the bone cement and the possibility of removing barium sulfate to reduce the 
toxicity of the material. Radiological images from X-ray and Ultrasonography are evaluated 


in cement implants in vertebrae. 


21 METHODS 


2.1 In vitro anatomical sample preparation 


The separation of a section of a pig vertebral column, in a special cut provided by 
a meat market, was performed in vitro embed in an equivalent muscle tissue, maintaining 
anatomically its structure without considering distinct anthropometric characteristics of the 
model. The structure was immediately cooled to -18 C to reduce the effect of decomposition. 


2.2 Synthesis of the composite 


The composite was prepared in cold based on PMMA, HAp, barium sulfate so that 
the dilution held excess water as dissolvent. The solution in deionized water was done 
primarily to modulate the phases of the acrylic in order to increase the polymerization time 
by mean of stirring of the mixture. Both PMMA as the instrumentis in the mixture were also 
cooled due to polymerizing effect that is proportional to the temperature increase. The HAp 
was synthesized by the sol-gel method according to (LEGEROS; LIN; ROHANIZADEH; 
MIJARES et al., 2003; MONTANO; CAMPOS; LEMOS; YOSHIDA et al., 2020) after to 
mix the components the solution was left for 24 hours in a closed beaker to force the 
precipitate's product, the nucleation and formation of colloids. The compounds used for the 
synthesis of HAp were 3.937 g of Ca(NO,),*H,O, 0.69 mL of H.PO,, 2 mL of CH (OH) and 
water as solvent in excess. Subsequently the sample was heated in an oven ramped from 
room temperature to 100 “ºC for 22 hours as follows: ffom room temperature to 80º Cat 
a rate of 0.306 “C.min”, then an isotherm 360 minutes, subsequently to 100 “ºC a ramp at 
the rate of 0.333 “C-min! and finally an isotherm 720 minutes. The next day, the sample 
was heated from room temperature to 720 “Cata rate of 6 “C-min”, then 60 min isotherm. 
The HAp powder is macerated and mixed in different proportions with the PMMA (Metil Etil 
Metacrilate (C,H,O,)) to produce 0.5 g of bone cementin a binary system (1-x)PMMA-xHAp, 
where [x] is the concentration of HAp in the mixture as shown in Table 1. The concentrations 


(1) and (2) were presented in this preparation only. 
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Nº Concentration [x] PMMA (C,H,0,) HAp(Ca,(PO,),(0H)) 
0.00000 0.50000 0.00000 
2 0.50000 0.08310 0.41690 


Table 1. Concentrations [x] bone cement to produce 0.5 g of the binary system (1- x)PMMA- 
xHAp. 


2.3 Cement injection 


The cement was injected into vertebral models designed in the laboratory drilling 
the vertebral body with a bent 45º to the midsagittal plane simulating drilling needle 


Vertebroplasty kit. 


2.4 Radiological contrast 


The X-ray machine used was the BR 100, a transportable model with performance 
100 mA and 90 kV. The BR 100 is composed with an X-ray tube sealed and linked to 
a telescope that is connected to the apparatus column than have an arm which allows 
depending on the movement necessary adjustment of the distance between tube and film 
that is studied radiographically. It also has a command table with the control elements 
necessary for desired voltage and intensity. According to the radiological techniques, 
exposure time of 0.50 s, current of 70 mA, voltage of 60 kV and DFF of 85 cm were applied 
in thoracic vertebrae. After the development processing, films were digitized to study the 
radiological contrast in each of the vertebrae that was implanted the bone cement. Contrast 


was evaluated in different poinis. 


2.5 Ultrasonography images 


The GE Healthcare ultrasound equipment LOGIQe REF. 5199704 made in China 
in 2008 was used for image guided in the Vertebroplasty simulating procedure with the pig 
vertebral column. One 12L-RS probe was used in B-cine mode with a signal of 12 MHz and 
with a maximum depth of 7 cm for the configuration abdominal type examination in addition 
before to the application of the contact gel for pig column model in vitro. In the second image 
signal was 10 MHz and the depth of field of ultrasound was defined to 4 cm. 


31 RESULTS 


The images of the phantom taken with radiology team have good contrast even 
though the paris are embedded in a tissue simulating equivalent soft tissue. According to Fig. 
1(a), an anterior-posterior frontal shot was made, in which the X-ray beam had to overcome 
apart phantom also the styrofoam lid container 2 cm that contains. The polystyrene and the 


surrounding soft equivalent tissue have a lower attenuation coefficient than to bone tissues. 
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In Fig. 1 (b), a side-sagittal image was shooting. The x- rays had to overcome two lids that 
are part of styrofoam container with 2 cm every one, maintaining good contrast. In Figi(c) 
and Fig. 1(d), images provided by ultrasound equipment, simulating the lumbar puncture 
protocol, are presented. A good contrast of the spinous processes, transverse and upper 
and lower joints, is observed. Just as the needle is guided by the spinous processes in the 
lumbar punciure for patients in gestation, it is possible to do the same in Vertebroplasty, with 
the difficulty of invading the field temporarily while the needle is inserted into the body of the 


vertebra. However, to control the flow of cement during insertion fluoroscopy is necessary. 





(a) (D) (c) 


Figure 1. X-ray and ultrasound images in vitro ffom vertebrae model is presented. In (a) and 
(b), frontal and sagittal radiological images are observed, respectively. In (c) and (d), there are 
ultrasound images; and, (d), a dynamic anomaly shows the insertion of the needle. 





Figure 2. Radiological images of four vertebral models. In (a) the vertebrae has molten lead. (b) 
PMMA injected. (c) PMMA + 10% BaSO, were injected. (d) x = 0.5 of the bone cement. 


The following images as shown in Fig. 2 were taken from four vertebral bodies. In 
Fig. 2(a), one of them is observed with an inlay of a small mass of lead in the linear shape 
that will be used as reference point to evaluate the radiological contrast of the other three 
images. In Fig. 2(b), a vertebral body with an abnormality is observed corresponding to 
the perforation made to the vertebral body in which PMMA was injected. In Fig. 2(c), it 
was injected PMMA with 10% BaSO, and no abnormalities indicating that the cement was 
largely diffused in bone piece. In Fig. 2(d), polymerizing bone cement (1-x) PMMA-xHAp 
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with a concentration [x], x = 0.5 HAp in the mix was injected. As in Fig. 2(c) the cement is 
well spread in the bone piece and with better contrasts which identifiable in the area of the 
perforation of the vertebral body. An interesting aspect was that, despite BaSO, suppress in 
the bone cement, this still shows a good radiological contrast. 

The physical distribution of this biomaterial was observed through X-ray images after 
application of the cementis in liquid phase. It has been spread 5 up to 20 mm into the 
vertebrae from the inject point. It demonstrated that such biomaterial can be incorporated 
adequately and a large amount of the cold cement can be well distributed in the vertebral 
structure. The spatial distribution of cement in regions Ill and IV of each vertebra is optimal 
intervened in the spongy tissue. Radiological response of the cement was demonstrated on 
the radiological films. 


41 CONCLUSIONS 


The spatial distribution of this material was observed through radiological images 
obtained after cement application showing the possibility of incorporating adequate amount 
of cement mass in the bone structure. It guaranteed the radiopaque nature of the bone 
cement. After cross cutting in some of the vertebrae, it was observed that the compound 
was solidified after being injected. Also, the spreading of the material in bone tissue was 
larger in the spongy bone tissue. A pair of puncture on the lumbar vertebrae can reduce the 
unwanted effect of extravasation due to the in situ decompression. 

Since these compounds are used for the purpose of be absorbed and reabsorbed 
bone tissues, it is very important to reduce toxicity suppressing as much as possible the 
substances which are not completely biocompatible with organic tissues, as radiopaque 
compounds used in medical imaging, without affecting significantly their action. All biological 
effects that may be induced by the excess of HAp must be considered in order to calculate 
the proportions of HAp in the cement. In addition, the desirable porosity for migration of 
bone cells is required. Also, porosity is essential to spread the cement and facilitate the join 
of fractured pieces, as the case is required. 

Percutaneous column procedures guided by fluoroscopy have been already 
proposed. One of the major problems in radioactive cemenis is directly related to vertebrae 
dose because dose must be planned with high precision to avoid negative responses 
subsequent to radiotherapeutic treatmenis. In the Ultrasound image the commercial teams 
gradually have improved their services and today are already possible to reconstruct live 
images in 3D which is known as the 4D mode. Thus, fluoroscopic image remains needed to 
avoid extravasation of the cement outside the vertebral body. 
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CAPÍTULO 16 


PROJETO DE UM PADRÃO UNIVERSAL DE BAIXO 
CUSTO PARA CALIBRAÇÃO DE EQUIPAMENTOS E 
INSPEÇÃO DE SOLDAGEM 


Data de aceite: 04/01/2021 welding inspection processes, it is necessary 
Data de submissão: 30/09/2020 that the equipment is properly calibrated. The 
objective of this work is to present the design of 
a multipurpose part to assist welding inspectors. 
The project unifies the main defecis in welding 
with the main functionalites of the welding 
gauges, also counting on a gauge of cuts and 
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http:/attes.cnpq.br/3888719768059144 electrodes and wires, in addition to enabling the 
visual characterization of materials. 
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11 INTRODUÇÃO 


Os processos de soldagem são 


RESUMO: Para garantir a confiabilidade amplamente utilizados na indústria tendo em 


dos processos de inspeção de soldagem, é vista que oferecem união permanente entre 
necessário que os equipamentos estejam elementos de diversas formas, através de 
devidamente calibrados. O objetivo deste diferentes fontes de energia que possibilitam 
trabalho é apresentar o projeto de uma peça a variação de parâmetros nos processos de 


multiuso para auxiliar inspetores de solda. O 
projeto unifica os principais defeitos em soldagem 
com as principais funcionalidades dos calibres 
de solda, contando também com um calibre 


acordo com a necessidade de cada projeto. 
É possível que ocorram defeitos nas juntas 


soldadas, causados na maioria das vezes 


de cortes e gabaritos para eletrodos e arames, pela escolha incorreta dos parâmetros de 
além de possibilitar a caracterização visual de soldagem (CHIAVERINI, 1986; MODENESI, 
materiais. MARQUES, 2000; KIMINAMI et al., 2013). A fim 


PALAVRAS-CHAVE: Soldagem, Calibração, 
Calibre Universal, Inspeção de Solda. 


de evitar que estes defeitos causem ruptura ou 
inutilização das peças ou estruturas produzidas, 


é importante que haja inspeção de solda. 
DESIGN OF A LOW COST UNIVERSAL A inspeção tem a finalidade de localizar 
STANDARD FOR EQUIPMENT 


CALIBRATION AND WELDING 
INSPECTION qualidade nos processos de fabricação, e 


os defeitos e garantir que seja mantida a 


ABSTRACT: To guarantee the reliability of the | pode ser realizada visualmente com o auxílio 
de instrumentos como calibres de solda, 
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espelhos e câmeras, e/ou com a utilização de equipamentos eletrônicos, como aparelhos 
de ultrassom, radiografia, e até mesmo geradores de campo elétrico em ensaios com 
partículas magnéticas (DE SOUZA et al., 2007). 

Nas técnicas de inspeção visual, o calibre de solda é muito utilizado pois auxilia o 
inspetor em diversas medições, tais como aberturas de chanfros, comprimento e reforço 
de cordões, espessuras de chapas, ângulos de bisel, profundidade de chanfros, altura de 
garganta e perna de solda, entre outros (INSIZE, 2013; CARBOGRAFITE, 2018). Existem 
diversos modelos de calibres disponíveis no mercado, alguns com função única como o 
medidor de perna e garganta em superfícies ortogonais (Figura 1a). Outros equipamentos 
conhecidos como calibres universais, que unificam o maior número de funções como 
medição de ângulos, espessuras, garganta e perna, profundidade de aberturas e 
comprimento de cordões (Figura 1b). 

Quanto aos equipamentos utilizados nas demais técnicas de inspeção, é necessário 
que estejam bem calibrados, pois o funcionamento correto destes equipamentos garante 
que os defeitos sejam encontrados, sem que haja erro em leituras, o que poderia apontar 
defeitos inexistentes ou deixar de apontar defeitos existentes. Aconselha-se que a 
calibração destes instrumentos seja realizada a cada seis meses por empresa confiável e 
de acordo com a norma ISO 9000:2000 (GIMENES; URTADO, 2011). 
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Figura 1. (a) - calibre de solda com função única (INSIZE, 2013) e (b) - Calibre de solda 
universal (Carbografite, 2018) 


O objetivo deste trabalho é apresentar o projeto de uma peça multiuso a ser utilizada 
por inspetores de solda tanto em inspeções visuais quanto na calibração de equipamentos 
de inspeção, tendo em vista que a peça contempla um calibre de solda universal e uma 
chapa com os principais defeitos em cordões de solda, causados propositalmente, para 
que os equipamentos de inspeção os localizem, comprovando seu funcionamento. Embora 
existam calibres universais no mercado, não há atualmente uma peça que contemple as 
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principais funcionalidades de calibres em conjunto com a possibilidade de calibração de 
equipamentos de inspeção. 

A peça apresenta ainda funcionalidades como medidores de diâmetros de eletrodos 
(TIG e ER) e arames de solda (MIG/MAG e SAS), além de informações pertinentes à área 


de soldagem. 


2 | MATERIAIS E MÉTODOS 


Algumas etapas foram seguidas para idealização do projeto. Primeiro, houve a 
definição do design da peça, sendo esta com três espessuras diferentes, contendo um tubo 
removível e um calibre universal. Posteriormente foi definido que os materiais utilizados 
seriam aço 1020 na chapa e aço inoxidável nas partes removíveis (tubo e calibre). Logo 
após, as dimensões foram estabelecidas, bem como quais defeitos seriam simulados. 
Houve também a cotação dos materiais para fabricação e uma estimativa do custo de mão 
de obra. Embora esteja pronto para ser fabricado, não houve a confecção de um protótipo 
do projeto em razão da pandemia mundial (COVID-19). 

Quanto ao design, a peça terá três espessuras diferentes para possibilitar a 
utilização de cordões multipasse nas partes de maior espessura, e ao mesmo tempo trazer 
economia de material com a redução da espessura em partes nas quais não foi necessária 
a utilização de multipasse. O calibre universal foi posicionado na parte de menor espessura 
assim como o tubo. A chapa conterá os 10 defeitos mais comuns em soldagem, para que o 
maior número de equipamentos de inspeção possam ser calibrados através destes defeitos. 
Foram utilizados diferentes cortes nas laterais com diferentes qualidades, servindo também 
como calibre de cortes. 

Sob a chapa foram gravadas à laser algumas informações pertinentes à área 
de soldagem, e a chapa também contará com amostras de materiais com diferentes 
colorações afim de permitir a identificação visual destes materiais, proporcionando que o 
inspetor reconheça os metais de base para solda. O material utilizado na chapa foi o aço 
1020, muito comum na indústria, assim como todos os aços com baixo teor de carbono, 
o que possibilita a geração dos defeitos em um material muito próximo aos quais surgirão 
os defeitos reais. Além disso, os aços de baixo carbono também apresentam ótimo custo 
benefício. 

Para o calibre e o tubo, partes removíveis que possivelmente terão mais contato 
com as mãos dos inspetores, deverá ser utilizado aço inoxidável, material predominante 
em todos os calibres encontrados no mercado, principalmente pela característica anti- 
oxidação. Quanto à oxidação da chapa, esta deverá ser evitada através da galvanização 


por zinco, por exemplo, ou outros métodos como pintura ou revestimento. 
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31 RESULTADOS E DISCUSSÃO 


Após o levantamento das ideias sobre quais seriam as dimensões e os componentes 
do projeto, definiu-se que as espessuras da chapa foram de 19,05 mm (0,75 in) 12,70 mm 
(0,50 in) e 6,35 mm (0,25 in). O calibre com largura de 70 mm e comprimento de 100 mm, 
contará com uma régua com escala de 100 mm, um medidor de ângulos com escala de 
0º a 180º, um medidor para eletrodos revestidos com diâmetros de 1,60 mm a 6,00 mm, 
um medidor de arames para solda MIG/MAG com diâmetros de 0,80 mm a 1,60 mm, um 
medidor de profundidade com escala de 15 mm e uma lâmina deslizante para medição 
de garganta, perna e espessuras com escala de 50 mm. As dimensões do calibre são 
demonstradas na Figura 2a (vista frontal) e na Figura 2b (vista posterior). 


HINO (b) 
ESCALA É 





Figura 2. (a) - vista frontal do calibre de solda e (b) - Calibre de solda universal 
(CARBOGRAFITE, 2018) 


A chapa também apresentará informações gravadas à laser, sendo estas 
informações: Diagrama Ferro-Carbono (Figura 3a), tabela com os principais tipos de 
chanfros e juntas nas quais são utilizados (Figura 3b), identificação por cor dos eletrodos 
TIG (Tabela 1), códigos e nomes dos defeitos produzidos nos cordões conforme a norma 
ISO 6520-1 (ISO, 1998) (Tabela 2), número de referência dos principais processos (Tabela 
3) e exemplo de identificação do código da nomenclatura dos eletrodos revestidos (Figura 
30). 

Ainda na chapa, na parte de maior espessura haverá a inclusão de alguns materiais 
utilizados na indústria com diferentes colorações para possibilitar a caracterização visual, 
sendo estes: ouro, cobre, alumínio, níquel, chumbo, bronze, latão e ferro. A chapa terá 300 


mm de largura e comprimento de 480 mm, tal como demonstrado na Figura 4. 
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[Eletrodo | a 


— | Indicam resistência a tração X 
1000 psi 





a e Refere-se à posição de soldagem 
——— (1 = Todas as posições, 2 = horizontal 
sd —o e plana, 3 = plana, 4 = plana, 
XX (X) Y E o sobrecarga, horizontal, vertical 
descendente) 


” | Indica o grau de utilização do eletrodo. 
Ex.: o tipo de corrente e o tipo de 
revestimento 








- | Composição química do depósito de 
soldagem 





Figura 3. (a) - diagrama Ferro-Carbono (CALLISTER, 2002), (b) - tipos de juntas e estilos de 
chanfros (MODENESI, 2008) e (c) - nomenclatura de eletrodos revestidos (FORTES, 2005) 


Designação Composição (Impurezas Totais < 0,10%) 
pi e 


o EWTh2 | ThO,:1,70-2,20% | Vermelha | 


E E 


Tabela 1. Identificação dos eletrodos TIG por cor (Adaptado de FONSECA, 2004) 








Tabela 2. Códigos de identificação dos defeitos em soldagem, norma ISO 6520-1 (ISO, 1998) 
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Soldagem por arco voltaico 35 | Arco metálico e gás ativo (MAG) 


a Arco metálico e gás ativo com eletrodo 
101 Arco metáli 
0 co metálico 6 Favestido 
1 Aicasjolisies manual q Arco metálico e gás inerte com eletrodo 
revestido 


12 [Acosbmero uai [GésinenodoMngetino iG) 
[15 [Acomeistcoogásprottor [15 [Sodaagama 


Tabela 3. Número de referência dos principais processos (arco elétrico) (Adaptado de 
FISCHER et al., 2011) 





480,00 
100,00 250,00 130,00 


300,00 





Figura 4. Vista superior da peça 


As laterais foram utilizadas para demonstrar as superfícies de quatro diferentes tipos 
de cortes, sendo estes cortes à plasma, oxicorte, à laser e por jato d'água. Nas laterais 
com 450 mm, foram utilizados os cortes à plasma e oxicorte com diferenças de qualidade 
entre as diferentes espessuras, na maior espessura o corte foi de maior qualidade, na 
espessura média foi de qualidade regular e na menor espessura a qualidade foi inferior, 
como exemplifica a Figura 5. 
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Figura 5. Exemplo de cortes laterais. 


Nas laterais com 350 mm foram utilizados os cortes à laser e por jato d'água, além 
da gravação à laser de escalas padrão em mm e em polegadas. Na Figura 6 é possível 
notar o exemplo de escala em mm. 





Figura 6. Escala graduada lateral em mm 


A montagem final da peça pode ser vista nas Figuras 7a, 7b e 7c. 
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(b) 





(c) 


Figura 7. (a) - vista lateral direita, (b) - vista lateral esquerda e (c) - montagem final da peça 
(vista em perspectiva isométrica) 


O projeto apresenta como principal vantagem o fato de ser multiuso, pois pode ser 
utilizado para diversas finalidades em soldagem e conta também com um calibre dos cortes 
mais utilizados além de amostras de materiais para identificação visual. O calibre universal 
anexado à peça possibilita que o inspetor utilize um único instrumento em inspeções 
visuais, reduzindo em volume e quantidade os instrumentos para transporte nas inspeções. 
Os defeitos presentes na peça facilitam a calibração de instrumentos de inspeção de solda, 
uma vez que estes podem ser testados através dos defeitos. 

A utilização de diferentes tipos de corte nas laterais com diferentes qualidades 
demonstra a aparência superficial resultante de cada processo, além de demonstrar 
também as qualidades dos processos mais utilizados (plasma e oxicorte), tornando possível 
que o inspetor verifique a qualidade do corte das peças por comparação. As amostras de 
materiais com diferentes colorações fazem com que tanto o inspetor quanto o soldador 
possam reconhecer o material soldado de forma rápida e fácil, o que pode evitar que sejam 
utilizados parâmetros incorretos durante o processo de soldagem. 

Outra vantagem do projeto é a segurança. A peça não apresenta riscos graves 
quanto à segurança de uso, os únicos riscos previsíveis são a possibilidade da peça se 
chocar contra o corpo do operador e o risco ergonômico, caso o operador não esteja atento 
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as regras de ergonomia ao erguer a peça do chão e carregá-la por um longo trajeto, por 
exemplo. 

Uma limitação do projeto é a dificuldade na realização de ensaios por líquido 
penetrante, o que descarta a verificação deste método de inspeção, pois o ensaio por líquido 
penetrante deixa resquícios que invalidariam outros testes com penetrante no mesmo local, 
logo, deveria ser utilizada uma nova superfície a cada ensaio, o que inviabiliza este tipo 
de ensaio na peça. A oxidação do aço 1020 também pode ser vista como limitação, porém 
é possível que este problema seja resolvido com a galvanização ou utilização de outros 
processos que impeçam a oxidação, como pintura ou revestimento. 

O custo de fabricação do projeto foi calculado com base na densidade dos materiais 
a serem utilizados. As densidades utilizadas nos cálculos foram de 8000 kg/mº para o aço 
inoxidável e 7870 kg/m? para o aço 1020, o peso aproximado da peça foi de 14,10 kg, 
desconsiderando o peso dos cordões de solda e de parafusos e chapas de fixação. 

O custo de fabricação da peça, incluindo os materiais e os serviços de gravação à 
laser, cortes e soldagem, seria de aproximadamente R$ 1500,00, com base em pesquisa 
realizada na cidade de Uberlândia-MG em maio de 2020. Este valor não engloba os custos 
de projeto realizado por engenheiro, bem como o custo para desenho em software ou 
lucros pretendidos com a possível comercialização. 

O projeto poderá impactar a indústria no sentido de que irá contribuir para a melhoria 
na qualidade das peças e estruturas produzidas, uma vez que auxiliará na calibração 
dos instrumentos de inspeção tornando-os mais precisos e confiáveis na localização de 
defeitos, possibilitando que estes sejam corrigidos, o que causa diminuição nos acidentes 
com a ruptura de estruturas, impactando também na sociedade. A peça também facilitaria 
o trabalho dos inspetores por oferecer em um único equipamento diversas funcionalidades, 
diminuindo o número de instrumentos isolados e oferecendo além da possibilidade de 
calibração dos instrumentos de diversos ensaios de inspeção, um calibre de solda universal 
para inspeções visuais. 

A calibração oferecida pela peça seria útil em todos os processos de soldagem nos 
quais ocorrem os defeitos ou descontinuidades citados (porosidade, trinca à quente, trinca 
a frio, cratera, falta de fusão e penetração, inclusões, mordeduras e respingos), pois esta 
funcionalidade é atingida através do apontamento dos defeitos existentes pelo aparelho 
de inspeção, podendo ser utilizado em qualquer material, não se limitando aos cordões 
realizados em aços com baixo teor de carbono. 


41 CONCLUSÃO 


Foi apresentado o projeto de uma peça multiuso que poderá ser utilizada por 
inspetores de solda em inspeções visuais e na calibração de equipamentos de inspeção. 
Através da análise do projeto é possível concluir que a peça apresenta os principais 
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defeitos em soldagem para calibração de instrumentos de inspeção, e que o calibre de 
solda incluso contempla as principais funcionalidades dos calibres de solda existentes no 
mercado, podendo portanto ser considerado um calibre universal. 

O projeto é multifuncional, devido à variedade de instrumentos passíveis de 
calibração, ao auxílio em inspeção visual, à possibilidade de identificação de materiais 
de base, à capacidade de identificação de diâmetros de eletrodos e arames de solda e à 
viabilidade em reconhecer a qualidade dos cortes mais utilizados na indústria, funcionando 


também como calibre de cortes. 
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